diff options
Diffstat (limited to 'trunk')
39 files changed, 0 insertions, 5708 deletions
diff --git a/trunk/chemCal/DESCRIPTION b/trunk/chemCal/DESCRIPTION deleted file mode 100644 index a4c8e82..0000000 --- a/trunk/chemCal/DESCRIPTION +++ /dev/null @@ -1,19 +0,0 @@ -Package: chemCal -Version: 0.2-02 -Date: 2012-09-11 -Title: Calibration functions for analytical chemistry -Author: Johannes Ranke <jranke@uni-bremen.de> -Maintainer: Johannes Ranke <jranke@uni-bremen.de> -Depends: R (>= 2.4.0) -Suggests: MASS -Description: chemCal provides simple functions for plotting linear - calibration functions and estimating standard errors for measurements - according to the Handbook of Chemometrics and Qualimetrics: Part A - by Massart et al. There are also functions estimating the limit - of detection (LOD) and limit of quantification (LOQ). - The functions work on model objects from - optionally weighted - linear - regression (lm) or robust linear regression (rlm from the MASS package). -License: GPL version 2 or newer -URL: http://www.r-project.org, - http://www.uft.uni-bremen.de/chemie/ranke/?page=chemCal, - http://kriemhild.uft.uni-bremen.de/viewcvs/?root=chemCal diff --git a/trunk/chemCal/INDEX b/trunk/chemCal/INDEX deleted file mode 100644 index f880eb6..0000000 --- a/trunk/chemCal/INDEX +++ /dev/null @@ -1,9 +0,0 @@ -calplot Plot calibration graphs from univariate linear - models -chemCal-package Calibration functions for analytical chemistry -din32645 Calibration data from DIN 32645 -inverse.predict Predict x from y for a linear calibration -lod Estimate a limit of detection (LOD) -loq Estimate a limit of quantification (LOQ) -massart97ex3 Calibration data from Massart et al. (1997), - example 3 diff --git a/trunk/chemCal/R/calfunctions.R b/trunk/chemCal/R/calfunctions.R deleted file mode 100644 index 6ce29f7..0000000 --- a/trunk/chemCal/R/calfunctions.R +++ /dev/null @@ -1,2 +0,0 @@ -powfunc <- function(x,a,b) a * x^b -ipowfunc <- function(y,a,b) (y/a)^1/b diff --git a/trunk/chemCal/R/calplot.R b/trunk/chemCal/R/calplot.R deleted file mode 100644 index 6aed9c0..0000000 --- a/trunk/chemCal/R/calplot.R +++ /dev/null @@ -1,80 +0,0 @@ -calplot <- function(object, - xlim = c("auto", "auto"), ylim = c("auto", "auto"), - xlab = "Concentration", ylab = "Response", alpha = 0.05, - varfunc = NULL) -{ - UseMethod("calplot") -} - -calplot.default <- function(object, - xlim = c("auto","auto"), ylim = c("auto","auto"), - xlab = "Concentration", ylab = "Response", - alpha=0.05, varfunc = NULL) -{ - stop("Calibration plots only implemented for univariate lm objects.") -} - -calplot.lm <- function(object, - xlim = c("auto","auto"), ylim = c("auto","auto"), - xlab = "Concentration", ylab = "Response", alpha=0.05, - varfunc = NULL) -{ - if (length(object$coef) > 2) - stop("More than one independent variable in your model - not implemented") - - if (alpha <= 0 | alpha >= 1) - stop("Alpha should be between 0 and 1 (exclusive)") - - m <- object - level <- 1 - alpha - y <- m$model[[1]] - x <- m$model[[2]] - if (xlim[1] == "auto") xlim[1] <- 0 - if (xlim[2] == "auto") xlim[2] <- max(x) - xlim <- as.numeric(xlim) - newdata <- list( - x = seq(from = xlim[[1]], to = xlim[[2]], length=250)) - names(newdata) <- names(m$model)[[2]] - if (is.null(varfunc)) { - varfunc <- if (length(m$weights)) { - function(variable) mean(m$weights) - } else function(variable) rep(1,250) - } - pred.lim <- predict(m, newdata, interval = "prediction", - level=level, weights.newdata = varfunc(m)) - conf.lim <- predict(m, newdata, interval = "confidence", - level=level) - yrange.auto <- range(c(0,pred.lim)) - if (ylim[1] == "auto") ylim[1] <- yrange.auto[1] - if (ylim[2] == "auto") ylim[2] <- yrange.auto[2] - plot(1, - type = "n", - xlab = xlab, - ylab = ylab, - xlim = as.numeric(xlim), - ylim = as.numeric(ylim) - ) - points(x,y, pch = 21, bg = "yellow") - matlines(newdata[[1]], pred.lim, lty = c(1, 4, 4), - col = c("black", "red", "red")) - if (length(object$weights) > 0) { - legend(min(x), - max(pred.lim, na.rm = TRUE), - legend = c("Fitted Line", "Confidence Bands"), - lty = c(1, 3), - lwd = 2, - col = c("black", "green4"), - horiz = FALSE, cex = 0.9, bg = "gray95") - } else { - matlines(newdata[[1]], conf.lim, lty = c(1, 3, 3), - col = c("black", "green4", "green4")) - legend(min(x), - max(pred.lim, na.rm = TRUE), - legend = c("Fitted Line", "Confidence Bands", - "Prediction Bands"), - lty = c(1, 3, 4), - lwd = 2, - col = c("black", "green4", "red"), - horiz = FALSE, cex = 0.9, bg = "gray95") - } -} diff --git a/trunk/chemCal/R/inverse.predict.lm.R b/trunk/chemCal/R/inverse.predict.lm.R deleted file mode 100644 index d57275c..0000000 --- a/trunk/chemCal/R/inverse.predict.lm.R +++ /dev/null @@ -1,115 +0,0 @@ -# This is an implementation of Equation (8.28) in the Handbook of Chemometrics -# and Qualimetrics, Part A, Massart et al (1997), page 200, validated with -# Example 8 on the same page, extended as specified in the package vignette - -inverse.predict <- function(object, newdata, ..., - ws = "auto", alpha = 0.05, var.s = "auto") -{ - UseMethod("inverse.predict") -} - -inverse.predict.default <- function(object, newdata, ..., - ws = "auto", alpha = 0.05, var.s = "auto") -{ - stop("Inverse prediction only implemented for univariate lm and nls objects.") -} - -inverse.predict.lm <- function(object, newdata, ..., - ws = "auto", alpha = 0.05, var.s = "auto") -{ - yname = names(object$model)[[1]] - xname = names(object$model)[[2]] - if (ws == "auto") { - ws <- ifelse(length(object$weights) > 0, mean(object$weights), 1) - } - if (length(object$weights) > 0) { - wx <- split(object$weights,object$model[[xname]]) - w <- sapply(wx,mean) - } else { - w <- rep(1,length(split(object$model[[yname]],object$model[[xname]]))) - } - .inverse.predict(object = object, newdata = newdata, - ws = ws, alpha = alpha, var.s = var.s, w = w, xname = xname, yname = yname) -} - -inverse.predict.rlm <- function(object, newdata, ..., - ws = "auto", alpha = 0.05, var.s = "auto") -{ - yname = names(object$model)[[1]] - xname = names(object$model)[[2]] - if (ws == "auto") { - ws <- mean(object$w) - } - wx <- split(object$weights,object$model[[xname]]) - w <- sapply(wx,mean) - .inverse.predict(object = object, newdata = newdata, - ws = ws, alpha = alpha, var.s = var.s, w = w, xname = xname, yname = yname) -} - -inverse.predict.nls <- function(object, newdata, ..., - ws = "auto", alpha = 0.05, var.s = "auto") -{ - yname = names(object$model)[[1]] - xname = names(object$model)[[2]] - if (ws == "auto") { - ws <- ifelse(length(object$weights) > 0, mean(object$weights), 1) - } - if (length(object$weights) > 0) { - wx <- split(object$weights,object$model[[xname]]) - w <- sapply(wx,mean) - } else { - w <- rep(1,length(split(object$model[[yname]],object$model[[xname]]))) - } - if (length(object$coef) > 2) - stop("More than one independent variable in your model - not implemented") -} - -.inverse.predict <- function(object, newdata, ws, alpha, var.s, w, xname, yname) -{ - if (length(object$coef) > 2) - stop("More than one independent variable in your model - not implemented") - - if (alpha <= 0 | alpha >= 1) - stop("Alpha should be between 0 and 1 (exclusive)") - - ybars <- mean(newdata) - m <- length(newdata) - - yx <- split(object$model[[yname]], object$model[[xname]]) - n <- length(yx) - df <- n - length(object$coef) - x <- as.numeric(names(yx)) - ybar <- sapply(yx,mean) - yhatx <- split(object$fitted.values, object$model[[xname]]) - yhat <- sapply(yhatx, mean) - se <- sqrt(sum(w * (ybar - yhat)^2)/df) - - if (var.s == "auto") { - var.s <- se^2/ws - } - - b1 <- object$coef[[xname]] - - ybarw <- sum(w * ybar)/sum(w) - -# This is the adapted form of equation 8.28 (see package vignette) - sxhats <- 1/b1 * sqrt( - (var.s / m) + - se^2 * (1/sum(w) + - (ybars - ybarw)^2 * sum(w) / - (b1^2 * (sum(w) * sum(w * x^2) - sum(w * x)^2))) - ) - - if (names(object$coef)[1] == "(Intercept)") { - b0 <- object$coef[["(Intercept)"]] - } else { - b0 <- 0 - } - - xs <- (ybars - b0) / b1 - t <- qt(1-0.5*alpha, n - 2) - conf <- t * sxhats - result <- list("Prediction"=xs,"Standard Error"=sxhats, - "Confidence"=conf, "Confidence Limits"=c(xs - conf, xs + conf)) - return(result) -} diff --git a/trunk/chemCal/R/lod.R b/trunk/chemCal/R/lod.R deleted file mode 100644 index f5bbb7d..0000000 --- a/trunk/chemCal/R/lod.R +++ /dev/null @@ -1,53 +0,0 @@ -lod <- function(object, ..., alpha = 0.05, beta = 0.05, method = "default") -{ - UseMethod("lod") -} - -lod.default <- function(object, ..., alpha = 0.05, beta = 0.05, method = "default") -{ - stop("lod is only implemented for univariate lm objects.") -} - -lod.lm <- function(object, ..., alpha = 0.05, beta = 0.05, method = "default") -{ - if (length(object$weights) > 0) { - stop(paste( - "\nThe detemination of a lod from calibration models obtained by", - "weighted linear regression requires confidence intervals for", - "predicted y values taking into account weights for the x values", - "from which the predictions are to be generated.", - "This is not supported by the internally used predict.lm method.", - sep = "\n" - )) - } - xname <- names(object$model)[[2]] - yname <- names(object$model)[[1]] - newdata <- data.frame(0) - names(newdata) <- xname - y0 <- predict(object, newdata, interval = "prediction", - level = 1 - 2 * alpha) - yc <- y0[[1,"upr"]] - if (method == "din") { - y0.d <- predict(object, newdata, interval = "prediction", - level = 1 - 2 * beta) - deltay <- y0.d[[1, "upr"]] - y0.d[[1, "fit"]] - lod.y <- yc + deltay - lod.x <- inverse.predict(object, lod.y)$Prediction - } else { - f <- function(x) { - newdata <- data.frame(x) - names(newdata) <- xname - pi.y <- predict(object, newdata, interval = "prediction", - level = 1 - 2 * beta) - yd <- pi.y[[1,"lwr"]] - (yd - yc)^2 - } - lod.x <- optimize(f,interval=c(0,max(object$model[[xname]])))$minimum - newdata <- data.frame(x = lod.x) - names(newdata) <- xname - lod.y <- predict(object, newdata) - } - lod <- list(lod.x, lod.y) - names(lod) <- c(xname, yname) - return(lod) -} diff --git a/trunk/chemCal/R/loq.R b/trunk/chemCal/R/loq.R deleted file mode 100644 index 5776096..0000000 --- a/trunk/chemCal/R/loq.R +++ /dev/null @@ -1,40 +0,0 @@ -loq <- function(object, ..., alpha = 0.05, k = 3, n = 1, w.loq = "auto", - var.loq = "auto") -{ - UseMethod("loq") -} - -loq.default <- function(object, ..., alpha = 0.05, k = 3, n = 1, w.loq = "auto", - var.loq = "auto") -{ - stop("loq is only implemented for univariate lm objects.") -} - -loq.lm <- function(object, ..., alpha = 0.05, k = 3, n = 1, w.loq = "auto", - var.loq = "auto") -{ - if (length(object$weights) > 0 && var.loq == "auto" && w.loq == "auto") { - stop(paste("If you are using a model from weighted regression,", - "you need to specify a reasonable approximation for the", - "weight (w.loq) or the variance (var.loq) at the", - "limit of quantification")) - } - xname <- names(object$model)[[2]] - yname <- names(object$model)[[1]] - f <- function(x) { - newdata <- data.frame(x = x) - names(newdata) <- xname - y <- predict(object, newdata) - p <- inverse.predict(object, rep(y, n), ws = w.loq, - var.s = var.loq, alpha = alpha) - (p[["Prediction"]] - k * p[["Confidence"]])^2 - } - tmp <- optimize(f,interval=c(0,max(object$model[[2]]))) - loq.x <- tmp$minimum - newdata <- data.frame(x = loq.x) - names(newdata) <- xname - loq.y <- predict(object, newdata) - loq <- list(loq.x, loq.y) - names(loq) <- c(xname, yname) - return(loq) -} diff --git a/trunk/chemCal/TODO b/trunk/chemCal/TODO deleted file mode 100644 index 1fc44a7..0000000 --- a/trunk/chemCal/TODO +++ /dev/null @@ -1,5 +0,0 @@ -- lod(): At the moment, it is not possible to calculate an lod for the - case of more than one replicates (m is fixed to 1). The formula - for the prediction of y from mean(x) in Massart et al, p. 433 - could be used to generalize this. -- Write methods for nonlinear calibration functions diff --git a/trunk/chemCal/data/din32645.rda b/trunk/chemCal/data/din32645.rda Binary files differdeleted file mode 100644 index 4e2913a..0000000 --- a/trunk/chemCal/data/din32645.rda +++ /dev/null diff --git a/trunk/chemCal/data/massart97ex1.rda b/trunk/chemCal/data/massart97ex1.rda Binary files differdeleted file mode 100644 index 1a6fd4b..0000000 --- a/trunk/chemCal/data/massart97ex1.rda +++ /dev/null diff --git a/trunk/chemCal/data/massart97ex3.rda b/trunk/chemCal/data/massart97ex3.rda Binary files differdeleted file mode 100644 index 7f60f3a..0000000 --- a/trunk/chemCal/data/massart97ex3.rda +++ /dev/null diff --git a/trunk/chemCal/man/calplot.lm.Rd b/trunk/chemCal/man/calplot.lm.Rd deleted file mode 100644 index b60a032..0000000 --- a/trunk/chemCal/man/calplot.lm.Rd +++ /dev/null @@ -1,69 +0,0 @@ -\name{calplot} -\alias{calplot} -\alias{calplot.default} -\alias{calplot.lm} -\title{Plot calibration graphs from univariate linear models} -\description{ - Produce graphics of calibration data, the fitted model as well - as confidence, and, for unweighted regression, prediction bands. -} -\usage{ - calplot(object, xlim = c("auto", "auto"), ylim = c("auto", "auto"), - xlab = "Concentration", ylab = "Response", alpha=0.05, varfunc = NULL) -} -\arguments{ - \item{object}{ - A univariate model object of class \code{\link{lm}} or - \code{\link[MASS:rlm]{rlm}} - with model formula \code{y ~ x} or \code{y ~ x - 1}. - } - \item{xlim}{ - The limits of the plot on the x axis. - } - \item{ylim}{ - The limits of the plot on the y axis. - } - \item{xlab}{ - The label of the x axis. - } - \item{ylab}{ - The label of the y axis. - } - \item{alpha}{ - The error tolerance level for the confidence and prediction bands. Note that this - includes both tails of the Gaussian distribution, unlike the alpha and beta parameters - used in \code{\link{lod}} (see note below). - } - \item{varfunc}{ - The variance function for generating the weights in the model. - Currently, this argument is ignored (see note below). - } -} -\value{ - A plot of the calibration data, of your fitted model as well as lines showing - the confidence limits. Prediction limits are only shown for models from - unweighted regression. -} -\note{ - Prediction bands for models from weighted linear regression require weights - for the data, for which responses should be predicted. Prediction intervals - using weights e.g. from a variance function are currently not supported by - the internally used function \code{\link{predict.lm}}, therefore, - \code{calplot} does not draw prediction bands for such models. - - It is possible to compare the \code{\link{calplot}} prediction bands with the - \code{\link{lod}} values if the \code{lod()} alpha and beta parameters are - half the value of the \code{calplot()} alpha parameter. - -} -\examples{ -data(massart97ex3) -m <- lm(y ~ x, data = massart97ex3) -calplot(m) -} -\author{ - Johannes Ranke - \email{jranke@uni-bremen.de} - \url{http://www.uft.uni-bremen.de/chemie/ranke} -} -\keyword{regression} diff --git a/trunk/chemCal/man/chemCal-package.Rd b/trunk/chemCal/man/chemCal-package.Rd deleted file mode 100644 index fb09cb1..0000000 --- a/trunk/chemCal/man/chemCal-package.Rd +++ /dev/null @@ -1,16 +0,0 @@ -\name{chemCal-package} -\alias{chemCal-package} -\docType{package} -\title{ - Calibration functions for analytical chemistry -} -\description{ - See \url{../DESCRIPTION} -} -\details{ - There is a package vignette located in \url{../doc/chemCal.pdf}. -} -\author{ - Author and Maintainer: Johannes Ranke <jranke@uni-bremen.de> -} -\keyword{manip} diff --git a/trunk/chemCal/man/din32645.Rd b/trunk/chemCal/man/din32645.Rd deleted file mode 100644 index cacbf07..0000000 --- a/trunk/chemCal/man/din32645.Rd +++ /dev/null @@ -1,61 +0,0 @@ -\name{din32645} -\docType{data} -\alias{din32645} -\title{Calibration data from DIN 32645} -\description{ - Sample dataset to test the package. -} -\usage{data(din32645)} -\format{ - A dataframe containing 10 rows of x and y values. -} -\examples{ -data(din32645) -m <- lm(y ~ x, data = din32645) -calplot(m) - -## Prediction of x with confidence interval -(prediction <- inverse.predict(m, 3500, alpha = 0.01)) - -# This should give 0.07434 according to test data from Dintest, which -# was collected from Procontrol 3.1 (isomehr GmbH) in this case -round(prediction$Confidence,5) - -## Critical value: -(crit <- lod(m, alpha = 0.01, beta = 0.5)) - -# According to DIN 32645, we should get 0.07 for the critical value -# (decision limit, "Nachweisgrenze") -round(crit$x, 2) -# and according to Dintest test data, we should get 0.0698 from -round(crit$x, 4) - -## Limit of detection (smallest detectable value given alpha and beta) -# In German, the smallest detectable value is the "Erfassungsgrenze", and we -# should get 0.14 according to DIN, which we achieve by using the method -# described in it: -lod.din <- lod(m, alpha = 0.01, beta = 0.01, method = "din") -round(lod.din$x, 2) - -## Limit of quantification -# This accords to the test data coming with the test data from Dintest again, -# except for the last digits of the value cited for Procontrol 3.1 (0.2121) -(loq <- loq(m, alpha = 0.01)) -round(loq$x,4) - -# A similar value is obtained using the approximation -# LQ = 3.04 * LC (Currie 1999, p. 120) -3.04 * lod(m,alpha = 0.01, beta = 0.5)$x -} -\references{ - DIN 32645 (equivalent to ISO 11843), Beuth Verlag, Berlin, 1994 - - Dintest. Plugin for MS Excel for evaluations of calibration data. Written - by Georg Schmitt, University of Heidelberg. - \url{http://www.rzuser.uni-heidelberg.de/~df6/download/dintest.htm} - - Currie, L. A. (1997) Nomenclature in evaluation of analytical methods including - detection and quantification capabilities (IUPAC Recommendations 1995). - Analytica Chimica Acta 391, 105 - 126. -} -\keyword{datasets} diff --git a/trunk/chemCal/man/inverse.predict.Rd b/trunk/chemCal/man/inverse.predict.Rd deleted file mode 100644 index 347d670..0000000 --- a/trunk/chemCal/man/inverse.predict.Rd +++ /dev/null @@ -1,69 +0,0 @@ -\name{inverse.predict} -\alias{inverse.predict} -\alias{inverse.predict.lm} -\alias{inverse.predict.rlm} -\alias{inverse.predict.default} -\title{Predict x from y for a linear calibration} -\usage{inverse.predict(object, newdata, \dots, - ws, alpha=0.05, var.s = "auto") -} -\arguments{ - \item{object}{ - A univariate model object of class \code{\link{lm}} or - \code{\link[MASS:rlm]{rlm}} - with model formula \code{y ~ x} or \code{y ~ x - 1}. - } - \item{newdata}{ - A vector of observed y values for one sample. - } - \item{\dots}{ - Placeholder for further arguments that might be needed by - future implementations. - } - \item{ws}{ - The weight attributed to the sample. This argument is obligatory - if \code{object} has weights. - } - \item{alpha}{ - The error tolerance level for the confidence interval to be reported. - } - \item{var.s}{ - The estimated variance of the sample measurements. The default is to take - the residual standard error from the calibration and to adjust it - using \code{ws}, if applicable. This means that \code{var.s} - overrides \code{ws}. - } -} -\value{ - A list containing the predicted x value, its standard error and a - confidence interval. -} -\description{ - This function predicts x values using a univariate linear model that has been - generated for the purpose of calibrating a measurement method. Prediction - intervals are given at the specified confidence level. - The calculation method was taken from Massart et al. (1997). In particular, - Equations 8.26 and 8.28 were combined in order to yield a general treatment - of inverse prediction for univariate linear models, taking into account - weights that have been used to create the linear model, and at the same - time providing the possibility to specify a precision in sample measurements - differing from the precision in standard samples used for the calibration. - This is elaborated in the package vignette. -} -\note{ - The function was validated with examples 7 and 8 from Massart et al. (1997). -} -\references{ - Massart, L.M, Vandenginste, B.G.M., Buydens, L.M.C., De Jong, S., Lewi, P.J., - Smeyers-Verbeke, J. (1997) Handbook of Chemometrics and Qualimetrics: Part A, - p. 200 -} -\examples{ -# This is example 7 from Chapter 8 in Massart et al. (1997) -data(massart97ex1) -m <- lm(y ~ x, data = massart97ex1) -inverse.predict(m, 15) # 6.1 +- 4.9 -inverse.predict(m, 90) # 43.9 +- 4.9 -inverse.predict(m, rep(90,5)) # 43.9 +- 3.2 -} -\keyword{manip} diff --git a/trunk/chemCal/man/ipowfunc.Rd b/trunk/chemCal/man/ipowfunc.Rd deleted file mode 100644 index e09e590..0000000 --- a/trunk/chemCal/man/ipowfunc.Rd +++ /dev/null @@ -1,33 +0,0 @@ -\name{ipowfunc} -\alias{ipowfunc} -\title{Power function} -\description{ - Inverse of the arithmetic power function \code{\link{powfunc}} used for - modelling univariate nonlinear calibration data. } -\usage{ - ipowfunc(x,a,b) -} -\arguments{ - \item{x}{ - Independent variable} - \item{a}{ - Coefficient} - \item{b}{ - Exponent} -} -\value{ - The result of evaluating the function - \deqn{f(x) = \frac{y}{a}^\frac{1}{b}}{f(x) = y/a^1/b} - which is the inverse of the function defined by \code{\link{powfunc}} -} -\author{ - Johannes Ranke - \email{jranke@uni-bremen.de} - \url{http://www.uft.uni-bremen.de/chemie/ranke} -} -\seealso{ - The original function \code{\link{powfunc}}. -} -\keyword{models} -\keyword{regression} -\keyword{nonlinear} diff --git a/trunk/chemCal/man/lod.Rd b/trunk/chemCal/man/lod.Rd deleted file mode 100644 index e468e1d..0000000 --- a/trunk/chemCal/man/lod.Rd +++ /dev/null @@ -1,83 +0,0 @@ -\name{lod} -\alias{lod} -\alias{lod.lm} -\alias{lod.rlm} -\alias{lod.default} -\title{Estimate a limit of detection (LOD)} -\usage{ - lod(object, \dots, alpha = 0.05, beta = 0.05, method = "default") -} -\arguments{ - \item{object}{ - A univariate model object of class \code{\link{lm}} or - \code{\link[MASS:rlm]{rlm}} - with model formula \code{y ~ x} or \code{y ~ x - 1}, - optionally from a weighted regression. - } - \item{\dots}{ - Placeholder for further arguments that might be needed by - future implementations. - } - \item{alpha}{ - The error tolerance for the decision limit (critical value). - } - \item{beta}{ - The error tolerance beta for the detection limit. - } - \item{method}{ - The default method uses a prediction interval at the LOD - for the estimation of the LOD, which obviously requires - iteration. This is described for example in Massart, p. 432 ff. - The \dQuote{din} method uses the prediction interval at - x = 0 as an approximation. - } -} -\value{ - A list containig the corresponding x and y values of the estimated limit of - detection of a model used for calibration. -} -\description{ - The decision limit (German: Nachweisgrenze) is defined as the signal or - analyte concentration that is significantly different from the blank signal - with a first order error alpha (one-sided significance test). - The detection limit, or more precise, the minimum detectable value - (German: Erfassungsgrenze), is then defined as the signal or analyte - concentration where the probability that the signal is not detected although - the analyte is present (type II or false negative error), is beta (also a - one-sided significance test). -} -\note{ - - The default values for alpha and beta are the ones recommended by IUPAC. - - The estimation of the LOD in terms of the analyte amount/concentration - xD from the LOD in the signal domain SD is done by simply inverting the - calibration function (i.e. assuming a known calibration function). - - The calculation of a LOD from weighted calibration models requires - a weights argument for the internally used \code{\link{predict.lm}} - function, which is currently not supported in R. -} -\references{ - Massart, L.M, Vandenginste, B.G.M., Buydens, L.M.C., De Jong, S., Lewi, P.J., - Smeyers-Verbeke, J. (1997) Handbook of Chemometrics and Qualimetrics: Part A, - Chapter 13.7.8 - - J. Inczedy, T. Lengyel, and A.M. Ure (2002) International Union of Pure and - Applied Chemistry Compendium of Analytical Nomenclature: Definitive Rules. - Web edition. - - Currie, L. A. (1997) Nomenclature in evaluation of analytical methods including - detection and quantification capabilities (IUPAC Recommendations 1995). - Analytica Chimica Acta 391, 105 - 126. -} -\examples{ -data(din32645) -m <- lm(y ~ x, data = din32645) -lod(m) - -# The critical value (decision limit, German Nachweisgrenze) can be obtained -# by using beta = 0.5: -lod(m, alpha = 0.01, beta = 0.5) -} -\seealso{ - Examples for \code{\link{din32645}} -} -\keyword{manip} diff --git a/trunk/chemCal/man/loq.Rd b/trunk/chemCal/man/loq.Rd deleted file mode 100644 index 7541e77..0000000 --- a/trunk/chemCal/man/loq.Rd +++ /dev/null @@ -1,77 +0,0 @@ -\name{loq} -\alias{loq} -\alias{loq.lm} -\alias{loq.rlm} -\alias{loq.default} -\title{Estimate a limit of quantification (LOQ)} -\usage{ - loq(object, \dots, alpha = 0.05, k = 3, n = 1, w.loq = "auto", - var.loq = "auto") -} -\arguments{ - \item{object}{ - A univariate model object of class \code{\link{lm}} or - \code{\link[MASS:rlm]{rlm}} - with model formula \code{y ~ x} or \code{y ~ x - 1}, - optionally from a weighted regression. If weights are specified - in the model, either \code{w.loq} or \code{var.loq} have to - be specified. - } - \item{alpha}{ - The error tolerance for the prediction of x values in the calculation. - } - \item{\dots}{ - Placeholder for further arguments that might be needed by - future implementations. - } - \item{k}{ - The inverse of the maximum relative error tolerated at the - desired LOQ. - } - \item{n}{ - The number of replicate measurements for which the LOQ should be - specified. - } - \item{w.loq}{ - The weight that should be attributed to the LOQ. Defaults - to one for unweighted regression, and to the mean of the weights - for weighted regression. See \code{\link{massart97ex3}} for - an example how to take advantage of knowledge about the - variance function. - } - \item{var.loq}{ - The approximate variance at the LOQ. The default value is - calculated from the model. - } -} -\value{ - The estimated limit of quantification for a model used for calibration. -} -\description{ - The limit of quantification is the x value, where the relative error - of the quantification given the calibration model reaches a prespecified - value 1/k. Thus, it is the solution of the equation - \deqn{L = k c(L)}{L = k * c(L)} - where c(L) is half of the length of the confidence interval at the limit L - (DIN 32645, equivalent to ISO 11843). c(L) is internally estimated by - \code{\link{inverse.predict}}, and L is obtained by iteration. -} -\note{ - - IUPAC recommends to base the LOQ on the standard deviation of the signal - where x = 0. - - The calculation of a LOQ based on weighted regression is non-standard - and therefore not tested. Feedback is welcome. -} -\examples{ -data(massart97ex3) -attach(massart97ex3) -m <- lm(y ~ x) -loq(m) - -# We can get better by using replicate measurements -loq(m, n = 3) -} -\seealso{ - Examples for \code{\link{din32645}} -} -\keyword{manip} diff --git a/trunk/chemCal/man/massart97ex1.Rd b/trunk/chemCal/man/massart97ex1.Rd deleted file mode 100644 index 44e1b85..0000000 --- a/trunk/chemCal/man/massart97ex1.Rd +++ /dev/null @@ -1,17 +0,0 @@ -\name{massart97ex1} -\docType{data} -\alias{massart97ex1} -\title{Calibration data from Massart et al. (1997), example 1} -\description{ - Sample dataset from p. 175 to test the package. -} -\usage{data(massart97ex1)} -\format{ - A dataframe containing 6 observations of x and y data. -} -\source{ - Massart, L.M, Vandenginste, B.G.M., Buydens, L.M.C., De Jong, S., Lewi, P.J., - Smeyers-Verbeke, J. (1997) Handbook of Chemometrics and Qualimetrics: Part A, - Chapter 8. -} -\keyword{datasets} diff --git a/trunk/chemCal/man/massart97ex3.Rd b/trunk/chemCal/man/massart97ex3.Rd deleted file mode 100644 index efdcf02..0000000 --- a/trunk/chemCal/man/massart97ex3.Rd +++ /dev/null @@ -1,51 +0,0 @@ -\name{massart97ex3} -\docType{data} -\alias{massart97ex3} -\title{Calibration data from Massart et al. (1997), example 3} -\description{ - Sample dataset from p. 188 to test the package. -} -\usage{data(massart97ex3)} -\format{ - A dataframe containing 6 levels of x values with 5 - observations of y for each level. -} -\examples{ -data(massart97ex3) -attach(massart97ex3) -yx <- split(y, x) -ybar <- sapply(yx, mean) -s <- round(sapply(yx, sd), digits = 2) -w <- round(1 / (s^2), digits = 3) -weights <- w[factor(x)] -m <- lm(y ~ x, w = weights) -calplot(m) - -# The following concords with the book p. 200 -inverse.predict(m, 15, ws = 1.67) # 5.9 +- 2.5 -inverse.predict(m, 90, ws = 0.145) # 44.1 +- 7.9 - -# The LOD is only calculated for models from unweighted regression -# with this version of chemCal -m0 <- lm(y ~ x) -lod(m0) - -# Limit of quantification from unweighted regression -loq(m0) - -# For calculating the limit of quantification from a model from weighted -# regression, we need to supply weights, internally used for inverse.predict -# If we are not using a variance function, we can use the weight from -# the above example as a first approximation (x = 15 is close to our -# loq approx 14 from above). -loq(m, w.loq = 1.67) -# The weight for the loq should therefore be derived at x = 7.3 instead -# of 15, but the graphical procedure of Massart (p. 201) to derive the -# variances on which the weights are based is quite inaccurate anyway. -} -\source{ - Massart, L.M, Vandenginste, B.G.M., Buydens, L.M.C., De Jong, S., Lewi, P.J., - Smeyers-Verbeke, J. (1997) Handbook of Chemometrics and Qualimetrics: Part A, - Chapter 8. -} -\keyword{datasets} diff --git a/trunk/chemCal/man/powfunc.Rd b/trunk/chemCal/man/powfunc.Rd deleted file mode 100644 index 73fe3b0..0000000 --- a/trunk/chemCal/man/powfunc.Rd +++ /dev/null @@ -1,32 +0,0 @@ -\name{powfunc} -\alias{powfunc} -\title{Power function} -\description{ - Arithmetic power function for modelling univariate nonlinear calibration data. -} -\usage{ - powfunc(x,a,b) -} -\arguments{ - \item{x}{ - Independent variable} - \item{a}{ - Coefficient} - \item{b}{ - Exponent} -} -\value{ - The result of evaluating the function - \deqn{f(x) = a x^b}{f(x) = a * x^b} -} -\author{ - Johannes Ranke - \email{jranke@uni-bremen.de} - \url{http://www.uft.uni-bremen.de/chemie/ranke} -} -\seealso{ - The inverse of this function \code{\link{ipowfunc}}. -} -\keyword{models} -\keyword{regression} -\keyword{nonlinear} diff --git a/trunk/chemCal/tests/din32645.R b/trunk/chemCal/tests/din32645.R deleted file mode 100644 index e5ffed7..0000000 --- a/trunk/chemCal/tests/din32645.R +++ /dev/null @@ -1,7 +0,0 @@ -require(chemCal) -data(din32645) -m <- lm(y ~ x, data = din32645) -inverse.predict(m, 3500, alpha = 0.01) -lod <- lod(m, alpha = 0.01, beta = 0.5) -lod(m, alpha = 0.01, beta = 0.01) -loq <- loq(m, alpha = 0.01) diff --git a/trunk/chemCal/tests/din32645.Rout.save b/trunk/chemCal/tests/din32645.Rout.save deleted file mode 100644 index c5ed5a7..0000000 --- a/trunk/chemCal/tests/din32645.Rout.save +++ /dev/null @@ -1,45 +0,0 @@ - -R : Copyright 2006, The R Foundation for Statistical Computing -Version 2.3.1 (2006-06-01) -ISBN 3-900051-07-0 - -R is free software and comes with ABSOLUTELY NO WARRANTY. -You are welcome to redistribute it under certain conditions. -Type 'license()' or 'licence()' for distribution details. - -R is a collaborative project with many contributors. -Type 'contributors()' for more information and -'citation()' on how to cite R or R packages in publications. - -Type 'demo()' for some demos, 'help()' for on-line help, or -'help.start()' for an HTML browser interface to help. -Type 'q()' to quit R. - -> require(chemCal) -Loading required package: chemCal -[1] TRUE -> data(din32645) -> m <- lm(y ~ x, data = din32645) -> inverse.predict(m, 3500, alpha = 0.01) -$Prediction -[1] 0.1054792 - -$`Standard Error` -[1] 0.02215619 - -$Confidence -[1] 0.07434261 - -$`Confidence Limits` -[1] 0.03113656 0.17982178 - -> lod <- lod(m, alpha = 0.01, beta = 0.5) -> lod(m, alpha = 0.01, beta = 0.01) -$x -[1] 0.132904 - -$y -[1] 3764.977 - -> loq <- loq(m, alpha = 0.01) -> diff --git a/trunk/chemCal/tests/massart97.R b/trunk/chemCal/tests/massart97.R deleted file mode 100644 index 00f837f..0000000 --- a/trunk/chemCal/tests/massart97.R +++ /dev/null @@ -1,25 +0,0 @@ -require(chemCal) -data(massart97ex1) -m <- lm(y ~ x, data = massart97ex1) -inverse.predict(m, 15) # 6.1 +- 4.9 -inverse.predict(m, 90) # 43.9 +- 4.9 -inverse.predict(m, rep(90,5)) # 43.9 +- 3.2 - -data(massart97ex3) -attach(massart97ex3) -yx <- split(y, x) -ybar <- sapply(yx, mean) -s <- round(sapply(yx, sd), digits = 2) -w <- round(1 / (s^2), digits = 3) -weights <- w[factor(x)] -m <- lm(y ~ x, w = weights) -#calplot(m) - -inverse.predict(m, 15, ws = 1.67) # 5.9 +- 2.5 -inverse.predict(m, 90, ws = 0.145) # 44.1 +- 7.9 - -m0 <- lm(y ~ x) -lod(m0) - -loq(m0) -loq(m, w.loq = 1.67) diff --git a/trunk/chemCal/tests/massart97.Rout.save b/trunk/chemCal/tests/massart97.Rout.save deleted file mode 100644 index cb113d0..0000000 --- a/trunk/chemCal/tests/massart97.Rout.save +++ /dev/null @@ -1,123 +0,0 @@ - -R : Copyright 2006, The R Foundation for Statistical Computing -Version 2.3.1 (2006-06-01) -ISBN 3-900051-07-0 - -R is free software and comes with ABSOLUTELY NO WARRANTY. -You are welcome to redistribute it under certain conditions. -Type 'license()' or 'licence()' for distribution details. - -R is a collaborative project with many contributors. -Type 'contributors()' for more information and -'citation()' on how to cite R or R packages in publications. - -Type 'demo()' for some demos, 'help()' for on-line help, or -'help.start()' for an HTML browser interface to help. -Type 'q()' to quit R. - -> require(chemCal) -Loading required package: chemCal -[1] TRUE -> data(massart97ex1) -> m <- lm(y ~ x, data = massart97ex1) -> inverse.predict(m, 15) # 6.1 +- 4.9 -$Prediction -[1] 6.09381 - -$`Standard Error` -[1] 1.767278 - -$Confidence -[1] 4.906751 - -$`Confidence Limits` -[1] 1.187059 11.000561 - -> inverse.predict(m, 90) # 43.9 +- 4.9 -$Prediction -[1] 43.93983 - -$`Standard Error` -[1] 1.767747 - -$Confidence -[1] 4.908053 - -$`Confidence Limits` -[1] 39.03178 48.84788 - -> inverse.predict(m, rep(90,5)) # 43.9 +- 3.2 -$Prediction -[1] 43.93983 - -$`Standard Error` -[1] 1.141204 - -$Confidence -[1] 3.168489 - -$`Confidence Limits` -[1] 40.77134 47.10832 - -> -> data(massart97ex3) -> attach(massart97ex3) -> yx <- split(y, x) -> ybar <- sapply(yx, mean) -> s <- round(sapply(yx, sd), digits = 2) -> w <- round(1 / (s^2), digits = 3) -> weights <- w[factor(x)] -> m <- lm(y ~ x, w = weights) -> #calplot(m) -> -> inverse.predict(m, 15, ws = 1.67) # 5.9 +- 2.5 -$Prediction -[1] 5.865367 - -$`Standard Error` -[1] 0.892611 - -$Confidence -[1] 2.478285 - -$`Confidence Limits` -[1] 3.387082 8.343652 - -> inverse.predict(m, 90, ws = 0.145) # 44.1 +- 7.9 -$Prediction -[1] 44.06025 - -$`Standard Error` -[1] 2.829162 - -$Confidence -[1] 7.855012 - -$`Confidence Limits` -[1] 36.20523 51.91526 - -> -> m0 <- lm(y ~ x) -> lod(m0) -$x -[1] 5.406637 - -$y -[1] 13.63822 - -> -> loq(m0) -$x -[1] 13.97767 - -$y -[1] 30.62355 - -> loq(m, w.loq = 1.67) -$x -[1] 7.346231 - -$y -[1] 17.90784 - -> diff --git a/trunk/chemCal/vignettes/Makefile b/trunk/chemCal/vignettes/Makefile deleted file mode 100644 index 8eca69e..0000000 --- a/trunk/chemCal/vignettes/Makefile +++ /dev/null @@ -1,33 +0,0 @@ -# Makefile for Sweave documents containing both Latex and R code -# Author: Johannes Ranke <jranke@uni-bremen.de> -# Last Change: 2006 Jun 23 -# based on the Makefile of Nicholas Lewin-Koh -# in turn based on work of Rouben Rostmaian -# SVN: $Id: Makefile.rnoweb 62 2006-05-24 08:30:59Z ranke $ - -RNWFILES = $(wildcard *.Rnw) -TARGETS = $(patsubst %.Rnw,%.tex,$(RNWFILE)) $(patsubst %.Rnw,%.pdf,$(RNWFILES)) - -RERUN = "(There were undefined references|Rerun to get (cross-references|the bars) right|Table widths have changed. Rerun LaTeX.|Linenumber reference failed)" -RERUNBIB = "No file.*\.bbl|Citation.*undefined" - -%.tex: %.Rnw - echo 'Sweave("$<")' | R --no-save --no-restore - -%.pdf: %.tex - pdflatex $< - egrep -c $(RERUNBIB) $*.log && (bibtex $*;pdflatex $<); true - egrep $(RERUN) $*.log && (pdflatex $<) ; true - egrep $(RERUN) $*.log && (pdflatex $<) ; true - -all: all-recursive $(TARGETS) - -clean: clean-recursive - rm -f *.aux *.log *.bbl *.blg *.brf *.cb *.ind *.idx *.ilg \ - *.inx *.ps *.dvi *.toc *.out *.lot *~ *.lof *.ttt *.fff - -all-recursive: - for dir in $(wildcard *); do if [ -d $$dir ] && [ -f $$dir/Makefile ]; then cd $$dir; $(MAKE) all; cd ..; fi; done - -clean-recursive: - for dir in $(wildcard *); do if [ -d $$dir ] && [ -f $$dir/Makefile ]; then cd $$dir; $(MAKE) clean; cd ..; fi; done diff --git a/trunk/chemCal/vignettes/Rplots.pdf b/trunk/chemCal/vignettes/Rplots.pdf Binary files differdeleted file mode 100644 index 4029222..0000000 --- a/trunk/chemCal/vignettes/Rplots.pdf +++ /dev/null diff --git a/trunk/chemCal/vignettes/Rplots.ps b/trunk/chemCal/vignettes/Rplots.ps deleted file mode 100644 index c5e68ec..0000000 --- a/trunk/chemCal/vignettes/Rplots.ps +++ /dev/null @@ -1,1970 +0,0 @@ -%!PS-Adobe-3.0 -%%DocumentNeededResources: font Helvetica -%%+ font Helvetica-Bold -%%+ font Helvetica-Oblique -%%+ font Helvetica-BoldOblique -%%+ font Symbol -%%DocumentMedia: a4 595 841 0 () () -%%Title: R Graphics Output -%%Creator: R Software -%%Pages: (atend) -%%Orientation: Landscape -%%BoundingBox: 18 18 577 824 -%%EndComments -%%BeginProlog -/bp { gs 595.00 0 translate 90 rotate gs } def -% begin .ps.prolog -/gs { gsave } def -/gr { grestore } def -/ep { showpage gr gr } def -/m { moveto } def -/l { rlineto } def -/np { newpath } def -/cp { closepath } def -/f { fill } def -/o { stroke } def -/c { newpath 0 360 arc } def -/r { 4 2 roll moveto 1 copy 3 -1 roll exch 0 exch rlineto 0 rlineto -1 mul 0 exch rlineto closepath } def -/p1 { stroke } def -/p2 { gsave bg setrgbcolor fill grestore newpath } def -/p3 { gsave bg setrgbcolor fill grestore stroke } def -/t { 6 -2 roll moveto gsave rotate - ps mul neg 0 2 1 roll rmoveto - 1 index stringwidth pop - mul neg 0 rmoveto show grestore } def -/cl { grestore gsave newpath 3 index 3 index moveto 1 index - 4 -1 roll lineto exch 1 index lineto lineto - closepath clip newpath } def -/rgb { setrgbcolor } def -/s { scalefont setfont } def -% end .ps.prolog -%%IncludeResource: font Helvetica -/Helvetica findfont -dup length dict begin - {1 index /FID ne {def} {pop pop} ifelse} forall - /Encoding ISOLatin1Encoding def - currentdict - end -/Font1 exch definefont pop -%%IncludeResource: font Helvetica-Bold -/Helvetica-Bold findfont -dup length dict begin - {1 index /FID ne {def} {pop pop} ifelse} forall - /Encoding ISOLatin1Encoding def - currentdict - end -/Font2 exch definefont pop -%%IncludeResource: font Helvetica-Oblique -/Helvetica-Oblique findfont -dup length dict begin - {1 index /FID ne {def} {pop pop} ifelse} forall - /Encoding ISOLatin1Encoding def - currentdict - end -/Font3 exch definefont pop -%%IncludeResource: font Helvetica-BoldOblique -/Helvetica-BoldOblique findfont -dup length dict begin - {1 index /FID ne {def} {pop pop} ifelse} forall - /Encoding ISOLatin1Encoding def - currentdict - end -/Font4 exch definefont pop -%%IncludeResource: font Symbol -/Symbol findfont -dup length dict begin - {1 index /FID ne {def} {pop pop} ifelse} forall - currentdict - end -/Font5 exch definefont pop -%%EndProlog -%%Page: 1 1 -bp -18.00 18.00 823.89 577.28 cl -0 0 0 rgb -0.75 setlinewidth -[] 0 setdash -1 setlinecap -1 setlinejoin -10.00 setmiterlimit -np -103.58 91.44 m -663.53 0 l -o -np -103.58 91.44 m -0 -7.20 l -o -np -236.29 91.44 m -0 -7.20 l -o -np -368.99 91.44 m -0 -7.20 l -o -np -501.70 91.44 m -0 -7.20 l -o -np -634.40 91.44 m -0 -7.20 l -o -np -767.11 91.44 m -0 -7.20 l -o -/ps 12 def /Font1 findfont 12 s -103.58 65.52 (0) .5 0 0 t -236.29 65.52 (10) .5 0 0 t -368.99 65.52 (20) .5 0 0 t -501.70 65.52 (30) .5 0 0 t -634.40 65.52 (40) .5 0 0 t -767.11 65.52 (50) .5 0 0 t -np -77.04 119.83 m -0 352.62 l -o -np -77.04 119.83 m --7.20 0 l -o -np -77.04 190.35 m --7.20 0 l -o -np -77.04 260.88 m --7.20 0 l -o -np -77.04 331.40 m --7.20 0 l -o -np -77.04 401.93 m --7.20 0 l -o -np -77.04 472.45 m --7.20 0 l -o -59.76 119.83 (0) .5 0 90 t -59.76 190.35 (20) .5 0 90 t -59.76 260.88 (40) .5 0 90 t -59.76 331.40 (60) .5 0 90 t -59.76 401.93 (80) .5 0 90 t -59.76 472.45 (100) .5 0 90 t -np -77.04 91.44 m -716.61 0 l -0 426.80 l --716.61 0 l -0 -426.80 l -o -18.00 18.00 823.89 577.28 cl -/ps 12 def /Font1 findfont 12 s -0 0 0 rgb -435.34 36.72 (Concentration) .5 0 0 t -30.96 304.84 (Response) .5 0 90 t -77.04 91.44 793.65 518.24 cl -/bg { 1 1 0 } def -0 0 0 rgb -0.75 setlinewidth -[] 0 setdash -1 setlinecap -1 setlinejoin -10.00 setmiterlimit -103.58 133.93 2.70 c p3 -236.29 197.41 2.70 c p3 -368.99 274.98 2.70 c p3 -501.70 331.40 2.70 c p3 -634.40 384.30 2.70 c p3 -767.11 486.56 2.70 c p3 -103.58 130.41 2.70 c p3 -236.29 190.35 2.70 c p3 -368.99 282.03 2.70 c p3 -501.70 341.98 2.70 c p3 -634.40 405.45 2.70 c p3 -767.11 504.19 2.70 c p3 -103.58 133.93 2.70 c p3 -236.29 193.88 2.70 c p3 -368.99 278.51 2.70 c p3 -501.70 331.40 2.70 c p3 -634.40 398.40 2.70 c p3 -767.11 497.14 2.70 c p3 -103.58 137.46 2.70 c p3 -236.29 197.41 2.70 c p3 -368.99 274.98 2.70 c p3 -501.70 341.98 2.70 c p3 -634.40 394.87 2.70 c p3 -767.11 475.98 2.70 c p3 -103.58 133.93 2.70 c p3 -236.29 193.88 2.70 c p3 -368.99 274.98 2.70 c p3 -501.70 341.98 2.70 c p3 -634.40 391.35 2.70 c p3 -767.11 490.08 2.70 c p3 -np -103.58 130.14 m -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -370.06 270.46 lineto -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -636.53 410.78 lineto -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -o -1 0 0 rgb -0.75 setlinewidth -[ 0.00 4.00 3.00 4.00] 0 setdash -np -103.58 107.25 m -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.42 l -2.66 1.41 l -2.67 1.42 l -2.66 1.41 l -2.67 1.42 l -2.66 1.41 l -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.41 l -2.67 1.42 l -2.66 1.41 l -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.41 l -2.67 1.42 l -2.66 1.41 l -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.41 l -2.66 1.41 l -2.67 1.42 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.42 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.40 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.41 l -2.67 1.40 l -2.66 1.41 l -2.67 1.41 l -2.66 1.40 l -2.66 1.41 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.41 l -370.06 248.29 lineto -2.66 1.41 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.39 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.40 l -2.67 1.39 l -2.66 1.40 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.40 l -636.53 388.36 lineto -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.40 l -2.67 1.39 l -2.66 1.40 l -2.67 1.39 l -2.66 1.40 l -2.67 1.39 l -2.66 1.40 l -2.67 1.39 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.39 l -2.66 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.39 l -2.66 1.40 l -2.67 1.39 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -2.67 1.39 l -2.66 1.40 l -2.67 1.39 l -o -np -103.58 153.03 m -2.67 1.39 l -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -2.66 1.40 l -2.67 1.39 l -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.39 l -2.66 1.40 l -2.67 1.39 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.39 l -2.66 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.39 l -2.66 1.40 l -2.67 1.39 l -2.66 1.40 l -2.67 1.39 l -2.66 1.40 l -2.67 1.39 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.40 l -2.67 1.39 l -2.66 1.40 l -2.67 1.39 l -2.66 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.40 l -2.67 1.39 l -2.66 1.40 l -2.67 1.39 l -2.66 1.40 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.40 l -2.67 1.39 l -2.66 1.40 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.39 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -370.06 292.63 lineto -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.41 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.41 l -2.67 1.41 l -2.66 1.40 l -2.66 1.41 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.41 l -2.67 1.40 l -2.66 1.41 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.40 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -636.53 433.20 lineto -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.42 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.41 l -2.67 1.42 l -2.66 1.41 l -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.41 l -2.67 1.42 l -2.66 1.41 l -2.67 1.42 l -2.66 1.41 l -2.67 1.41 l -2.66 1.42 l -2.66 1.41 l -2.67 1.42 l -2.66 1.41 l -2.67 1.42 l -2.66 1.41 l -2.67 1.42 l -o -0 0 0 rgb -0.75 setlinewidth -[] 0 setdash -np -103.58 130.14 m -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -370.06 270.46 lineto -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -636.53 410.78 lineto -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -o -0 0.5451 0 rgb -0.75 setlinewidth -[ 0.00 4.00] 0 setdash -np -103.58 123.09 m -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.45 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.66 1.43 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.43 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.43 l -2.66 1.44 l -2.67 1.44 l -2.66 1.43 l -2.67 1.44 l -2.66 1.43 l -2.67 1.44 l -2.66 1.43 l -2.67 1.44 l -2.66 1.43 l -2.67 1.44 l -2.66 1.43 l -2.66 1.44 l -2.67 1.43 l -2.66 1.44 l -2.67 1.43 l -2.66 1.43 l -2.67 1.44 l -2.66 1.43 l -2.67 1.43 l -2.66 1.43 l -2.67 1.44 l -2.66 1.43 l -2.67 1.43 l -2.66 1.43 l -2.67 1.43 l -2.66 1.43 l -2.67 1.43 l -2.66 1.43 l -2.67 1.43 l -2.66 1.43 l -2.67 1.43 l -2.66 1.43 l -2.66 1.43 l -2.67 1.43 l -2.66 1.43 l -2.67 1.42 l -2.66 1.43 l -2.67 1.43 l -2.66 1.42 l -2.67 1.43 l -2.66 1.43 l -2.67 1.42 l -2.66 1.43 l -2.67 1.42 l -2.66 1.43 l -2.67 1.42 l -2.66 1.42 l -2.67 1.43 l -2.66 1.42 l -2.67 1.42 l -2.66 1.42 l -2.67 1.43 l -2.66 1.42 l -2.67 1.42 l -2.66 1.42 l -2.66 1.42 l -2.67 1.42 l -2.66 1.41 l -2.67 1.42 l -2.66 1.42 l -2.67 1.42 l -2.66 1.41 l -370.06 266.32 lineto -2.66 1.42 l -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.66 1.40 l -2.67 1.39 l -2.66 1.40 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.39 l -2.66 1.40 l -2.67 1.39 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -2.67 1.38 l -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -2.66 1.38 l -2.67 1.39 l -2.66 1.38 l -2.67 1.39 l -2.66 1.38 l -2.67 1.39 l -2.66 1.38 l -2.67 1.38 l -2.66 1.39 l -2.67 1.38 l -2.66 1.38 l -2.67 1.38 l -2.66 1.38 l -2.67 1.39 l -2.66 1.38 l -2.67 1.38 l -2.66 1.38 l -2.67 1.38 l -2.66 1.37 l -2.67 1.38 l -2.66 1.38 l -2.66 1.38 l -2.67 1.38 l -2.66 1.38 l -2.67 1.37 l -2.66 1.38 l -2.67 1.38 l -2.66 1.37 l -2.67 1.38 l -2.66 1.37 l -2.67 1.38 l -2.66 1.37 l -2.67 1.38 l -2.66 1.37 l -2.67 1.38 l -2.66 1.37 l -2.67 1.38 l -2.66 1.37 l -2.67 1.37 l -2.66 1.38 l -2.67 1.37 l -2.66 1.37 l -636.53 405.46 lineto -2.67 1.38 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.36 l -2.66 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.36 l -2.66 1.37 l -2.67 1.37 l -2.66 1.36 l -2.67 1.37 l -2.66 1.37 l -2.67 1.36 l -2.66 1.37 l -2.67 1.37 l -2.66 1.36 l -2.67 1.37 l -2.66 1.37 l -2.67 1.36 l -2.66 1.37 l -2.67 1.36 l -2.66 1.37 l -2.67 1.36 l -2.66 1.37 l -2.66 1.36 l -2.67 1.37 l -2.66 1.36 l -2.67 1.37 l -2.66 1.36 l -2.67 1.37 l -o -np -103.58 137.19 m -2.67 1.36 l -2.66 1.37 l -2.67 1.36 l -2.66 1.37 l -2.66 1.36 l -2.67 1.37 l -2.66 1.36 l -2.67 1.37 l -2.66 1.36 l -2.67 1.37 l -2.66 1.36 l -2.67 1.37 l -2.66 1.37 l -2.67 1.36 l -2.66 1.37 l -2.67 1.36 l -2.66 1.37 l -2.67 1.37 l -2.66 1.36 l -2.67 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.36 l -2.67 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.36 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.36 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.67 1.38 l -2.66 1.37 l -2.67 1.37 l -2.66 1.37 l -2.66 1.37 l -2.67 1.37 l -2.66 1.38 l -2.67 1.37 l -2.66 1.37 l -2.67 1.38 l -2.66 1.37 l -2.67 1.37 l -2.66 1.38 l -2.67 1.37 l -2.66 1.38 l -2.67 1.37 l -2.66 1.38 l -2.67 1.37 l -2.66 1.38 l -2.67 1.37 l -2.66 1.38 l -2.67 1.38 l -2.66 1.37 l -2.67 1.38 l -2.66 1.38 l -2.66 1.38 l -2.67 1.37 l -2.66 1.38 l -2.67 1.38 l -2.66 1.38 l -2.67 1.38 l -2.66 1.38 l -2.67 1.38 l -2.66 1.38 l -2.67 1.38 l -2.66 1.38 l -2.67 1.38 l -2.66 1.39 l -2.67 1.38 l -2.66 1.38 l -2.67 1.39 l -2.66 1.38 l -2.67 1.38 l -2.66 1.39 l -2.67 1.38 l -2.66 1.39 l -2.67 1.38 l -2.66 1.39 l -2.66 1.39 l -2.67 1.38 l -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -370.06 274.60 lineto -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -2.67 1.39 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.39 l -2.67 1.39 l -2.66 1.40 l -2.67 1.40 l -2.66 1.39 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.39 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.40 l -2.67 1.40 l -2.66 1.40 l -2.67 1.41 l -2.66 1.40 l -2.67 1.40 l -2.66 1.41 l -2.67 1.40 l -2.66 1.41 l -2.67 1.41 l -2.66 1.40 l -2.67 1.41 l -2.66 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.41 l -2.67 1.41 l -2.66 1.42 l -2.67 1.41 l -2.66 1.41 l -2.67 1.42 l -2.66 1.41 l -2.67 1.42 l -2.66 1.41 l -2.67 1.42 l -2.66 1.42 l -2.67 1.41 l -2.66 1.42 l -2.67 1.42 l -2.66 1.42 l -2.67 1.42 l -2.66 1.42 l -2.66 1.42 l -2.67 1.42 l -2.66 1.42 l -2.67 1.42 l -2.66 1.42 l -2.67 1.43 l -2.66 1.42 l -2.67 1.42 l -2.66 1.43 l -2.67 1.42 l -2.66 1.42 l -2.67 1.43 l -2.66 1.42 l -2.67 1.43 l -2.66 1.42 l -2.67 1.43 l -2.66 1.43 l -2.67 1.42 l -2.66 1.43 l -2.67 1.43 l -2.66 1.43 l -2.66 1.43 l -2.67 1.42 l -2.66 1.43 l -2.67 1.43 l -2.66 1.43 l -2.67 1.43 l -2.66 1.43 l -2.67 1.43 l -2.66 1.43 l -2.67 1.43 l -2.66 1.44 l -2.67 1.43 l -2.66 1.43 l -2.67 1.43 l -2.66 1.43 l -2.67 1.44 l -2.66 1.43 l -2.67 1.43 l -2.66 1.44 l -2.67 1.43 l -2.66 1.43 l -636.53 416.10 lineto -2.67 1.43 l -2.66 1.44 l -2.67 1.43 l -2.66 1.44 l -2.67 1.43 l -2.66 1.44 l -2.67 1.43 l -2.66 1.44 l -2.67 1.43 l -2.66 1.44 l -2.67 1.44 l -2.66 1.43 l -2.67 1.44 l -2.66 1.44 l -2.67 1.43 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.43 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.43 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.45 l -2.67 1.44 l -2.66 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.44 l -2.66 1.44 l -2.67 1.45 l -o -/bg { 0.9490 0.9490 0.9490 } def -0 0 0 rgb -0.75 setlinewidth -[] 0 setdash -103.58 502.43 133.62 -51.84 r p3 -1.50 setlinewidth -[] 0 setdash -np -113.30 489.47 m -19.44 0 l -o -0 0.5451 0 rgb -1.50 setlinewidth -[ 0.00 6.00] 0 setdash -np -113.30 476.51 m -19.44 0 l -o -1 0 0 rgb -1.50 setlinewidth -[ 0.00 6.00 4.50 6.00] 0 setdash -np -113.30 463.55 m -19.44 0 l -o -/ps 11 def /Font1 findfont 11 s -0 0 0 rgb -142.46 485.52 (Fitted Line) 0 0 0 t -142.46 472.56 (Confidence Bands) 0 0 0 t -142.46 459.60 (Prediction Bands) 0 0 0 t -ep -%%Page: 2 2 -bp -77.04 91.44 793.65 518.24 cl -18.00 18.00 823.89 577.28 cl -0 0 0 rgb -0.75 setlinewidth -[] 0 setdash -1 setlinecap -1 setlinejoin -10.00 setmiterlimit -np -84.00 91.44 m -669.65 0 l -o -np -84.00 91.44 m -0 -7.20 l -o -np -217.93 91.44 m -0 -7.20 l -o -np -351.86 91.44 m -0 -7.20 l -o -np -485.79 91.44 m -0 -7.20 l -o -np -619.72 91.44 m -0 -7.20 l -o -np -753.65 91.44 m -0 -7.20 l -o -/ps 12 def /Font1 findfont 12 s -84.00 65.52 (0) .5 0 0 t -217.93 65.52 (20) .5 0 0 t -351.86 65.52 (40) .5 0 0 t -485.79 65.52 (60) .5 0 0 t -619.72 65.52 (80) .5 0 0 t -753.65 65.52 (100) .5 0 0 t -np -77.04 107.25 m -0 378.00 l -o -np -77.04 107.25 m --7.20 0 l -o -np -77.04 233.25 m --7.20 0 l -o -np -77.04 359.25 m --7.20 0 l -o -np -77.04 485.25 m --7.20 0 l -o -59.76 107.25 (0.0) .5 0 90 t -59.76 233.25 (0.5) .5 0 90 t -59.76 359.25 (1.0) .5 0 90 t -59.76 485.25 (1.5) .5 0 90 t -np -77.04 91.44 m -716.61 0 l -0 426.80 l --716.61 0 l -0 -426.80 l -o -18.00 18.00 823.89 577.28 cl -/ps 12 def /Font1 findfont 12 s -0 0 0 rgb -435.34 36.72 (Fitted values) .5 0 0 t -0.75 setlinewidth -[] 0 setdash -1 setlinecap -1 setlinejoin -10.00 setmiterlimit -np -30.17 238.93 m --0.87 1.55 l -4.54 1.56 l --11.58 2.06 l -0 126.09 l -o -np -33.84 245.91 m --9.07 0 l -o -33.61 246.47 (S) 0 0 90 t -33.61 254.47 (t) 0 0 90 t -33.61 257.81 (a) 0 0 90 t -33.61 264.48 (n) 0 0 90 t -33.61 271.15 (d) 0 0 90 t -33.61 277.82 (a) 0 0 90 t -33.61 284.50 (r) 0 0 90 t -33.61 288.49 (d) 0 0 90 t -33.61 295.16 (i) 0 0 90 t -33.61 297.83 (z) 0 0 90 t -33.61 303.83 (e) 0 0 90 t -33.61 310.50 (d) 0 0 90 t -33.61 317.17 ( ) 0 0 90 t -33.61 320.51 (r) 0 0 90 t -33.61 324.50 (e) 0 0 90 t -33.61 331.18 (s) 0 0 90 t -33.61 337.18 (i) 0 0 90 t -33.61 339.84 (d) 0 0 90 t -33.61 346.51 (u) 0 0 90 t -33.61 353.18 (a) 0 0 90 t -33.61 359.86 (l) 0 0 90 t -33.61 362.52 (s) 0 0 90 t -np -33.84 369.08 m --9.07 0 l -o -77.04 91.44 793.65 518.24 cl -0 0 0 rgb -0.75 setlinewidth -[] 0 setdash -1 setlinecap -1 setlinejoin -10.00 setmiterlimit -103.58 262.03 2.70 c p1 -236.29 234.09 2.70 c p1 -368.99 283.13 2.70 c p1 -501.70 332.99 2.70 c p1 -634.40 502.43 2.70 c p1 -767.11 317.74 2.70 c p1 -103.58 148.43 2.70 c p1 -236.29 351.20 2.70 c p1 -368.99 378.99 2.70 c p1 -501.70 223.02 2.70 c p1 -634.40 268.15 2.70 c p1 -767.11 501.72 2.70 c p1 -103.58 262.03 2.70 c p1 -236.29 301.67 2.70 c p1 -368.99 336.13 2.70 c p1 -501.70 332.99 2.70 c p1 -634.40 370.53 2.70 c p1 -767.11 440.55 2.70 c p1 -103.58 322.23 2.70 c p1 -236.29 234.09 2.70 c p1 -368.99 283.13 2.70 c p1 -501.70 223.02 2.70 c p1 -634.40 408.96 2.70 c p1 -767.11 257.16 2.70 c p1 -103.58 262.03 2.70 c p1 -236.29 301.67 2.70 c p1 -368.99 283.13 2.70 c p1 -501.70 223.02 2.70 c p1 -634.40 443.02 2.70 c p1 -767.11 365.26 2.70 c p1 -1 0 0 rgb -np -103.58 255.58 m -0 0 l -0 0 l -0 0 l -0 0 l -132.71 27.87 l -0 0 l -0 0 l -0 0 l -0 0 l -132.70 9.01 l -0 0 l -0 0 l -0 0 l -0 0 l -132.71 24.97 l -0 0 l -0 0 l -0 0 l -0 0 l -132.70 36.59 l -0 0 l -0 0 l -0 0 l -0 0 l -132.71 43.75 l -0 0 l -0 0 l -0 0 l -0 0 l -o -18.00 18.00 823.89 577.28 cl -/ps 12 def /Font1 findfont 12 s -0 0 0 rgb -435.34 22.32 (lm\(y ~ x\)) .5 0 0 t -18.00 18.00 823.89 577.28 cl -/ps 12 def /Font1 findfont 12 s -0 0 0 rgb -435.34 524.72 (Scale-Location) .5 0 0 t -18.00 18.00 823.89 577.28 cl -/ps 9 def /Font1 findfont 9 s -0 0 0 rgb -638.00 500.27 (5) 0 0 0 t -763.51 499.57 (12) 1 0 0 t -638.00 440.86 (29) 0 0 0 t -ep -%%Trailer -%%Pages: 2 -%%EOF diff --git a/trunk/chemCal/vignettes/chemCal-001.eps b/trunk/chemCal/vignettes/chemCal-001.eps deleted file mode 100644 index acc5f4b..0000000 --- a/trunk/chemCal/vignettes/chemCal-001.eps +++ /dev/null @@ -1,1762 +0,0 @@ -%!PS-Adobe-3.0 -%%DocumentNeededResources: font Helvetica -%%+ font Helvetica-Bold -%%+ font Helvetica-Oblique -%%+ font Helvetica-BoldOblique -%%+ font Symbol -%%DocumentMedia: special 432 432 0 () () -%%Title: R Graphics Output -%%Creator: R Software -%%Pages: (atend) -%%BoundingBox: 0 0 432 432 -%%EndComments -%%BeginProlog -/bp { gs gs } def -% begin .ps.prolog -/gs { gsave } def -/gr { grestore } def -/ep { showpage gr gr } def -/m { moveto } def -/l { rlineto } def -/np { newpath } def -/cp { closepath } def -/f { fill } def -/o { stroke } def -/c { newpath 0 360 arc } def -/r { 4 2 roll moveto 1 copy 3 -1 roll exch 0 exch rlineto 0 rlineto -1 mul 0 exch rlineto closepath } def -/p1 { stroke } def -/p2 { gsave bg setrgbcolor fill grestore newpath } def -/p3 { gsave bg setrgbcolor fill grestore stroke } def -/t { 6 -2 roll moveto gsave rotate - ps mul neg 0 2 1 roll rmoveto - 1 index stringwidth pop - mul neg 0 rmoveto show grestore } def -/cl { grestore gsave newpath 3 index 3 index moveto 1 index - 4 -1 roll lineto exch 1 index lineto lineto - closepath clip newpath } def -/rgb { setrgbcolor } def -/s { scalefont setfont } def -% end .ps.prolog -%%IncludeResource: font Helvetica -/Helvetica findfont -dup length dict begin - {1 index /FID ne {def} {pop pop} ifelse} forall - /Encoding ISOLatin1Encoding def - currentdict - end -/Font1 exch definefont pop -%%IncludeResource: font Helvetica-Bold -/Helvetica-Bold findfont -dup length dict begin - {1 index /FID ne {def} {pop pop} ifelse} forall - /Encoding ISOLatin1Encoding def - currentdict - end -/Font2 exch definefont pop -%%IncludeResource: font Helvetica-Oblique -/Helvetica-Oblique findfont -dup length dict begin - {1 index /FID ne {def} {pop pop} ifelse} forall - /Encoding ISOLatin1Encoding def - currentdict - end -/Font3 exch definefont pop -%%IncludeResource: font Helvetica-BoldOblique -/Helvetica-BoldOblique findfont -dup length dict begin - {1 index /FID ne {def} {pop pop} ifelse} forall - /Encoding ISOLatin1Encoding def - currentdict - end -/Font4 exch definefont pop -%%IncludeResource: font Symbol -/Symbol findfont -dup length dict begin - {1 index /FID ne {def} {pop pop} ifelse} forall - currentdict - end -/Font5 exch definefont pop -%%EndProlog -%%Page: 1 1 -bp -0.00 0.00 432.00 432.00 cl -0 0 0 rgb -0.75 setlinewidth -[] 0 setdash -1 setlinecap -1 setlinejoin -10.00 setmiterlimit -np -71.73 73.44 m -317.34 0 l -o -np -71.73 73.44 m -0 -7.20 l -o -np -135.20 73.44 m -0 -7.20 l -o -np -198.67 73.44 m -0 -7.20 l -o -np -262.13 73.44 m -0 -7.20 l -o -np -325.60 73.44 m -0 -7.20 l -o -np -389.07 73.44 m -0 -7.20 l -o -/ps 12 def /Font1 findfont 12 s -71.73 47.52 (0) .5 0 0 t -135.20 47.52 (10) .5 0 0 t -198.67 47.52 (20) .5 0 0 t -262.13 47.52 (30) .5 0 0 t -325.60 47.52 (40) .5 0 0 t -389.07 47.52 (50) .5 0 0 t -np -59.04 93.36 m -0 247.47 l -o -np -59.04 93.36 m --7.20 0 l -o -np -59.04 142.86 m --7.20 0 l -o -np -59.04 192.35 m --7.20 0 l -o -np -59.04 241.84 m --7.20 0 l -o -np -59.04 291.34 m --7.20 0 l -o -np -59.04 340.83 m --7.20 0 l -o -41.76 93.36 (0) .5 0 90 t -41.76 142.86 (20) .5 0 90 t -41.76 192.35 (40) .5 0 90 t -41.76 241.84 (60) .5 0 90 t -41.76 291.34 (80) .5 0 90 t -41.76 340.83 (100) .5 0 90 t -np -59.04 73.44 m -342.72 0 l -0 299.52 l --342.72 0 l -0 -299.52 l -o -0.00 0.00 432.00 432.00 cl -/ps 12 def /Font1 findfont 12 s -0 0 0 rgb -230.40 18.72 (Concentration) .5 0 0 t -12.96 223.20 (Response) .5 0 90 t -59.04 73.44 401.76 372.96 cl -/bg { 1 1 0 } def -0 0 0 rgb -0.75 setlinewidth -[] 0 setdash -1 setlinecap -1 setlinejoin -10.00 setmiterlimit -71.73 103.26 2.70 c p3 -135.20 147.81 2.70 c p3 -198.67 202.25 2.70 c p3 -262.13 241.84 2.70 c p3 -325.60 278.96 2.70 c p3 -389.07 350.73 2.70 c p3 -71.73 100.79 2.70 c p3 -135.20 142.86 2.70 c p3 -198.67 207.20 2.70 c p3 -262.13 249.27 2.70 c p3 -325.60 293.81 2.70 c p3 -389.07 363.10 2.70 c p3 -71.73 103.26 2.70 c p3 -135.20 145.33 2.70 c p3 -198.67 204.72 2.70 c p3 -262.13 241.84 2.70 c p3 -325.60 288.86 2.70 c p3 -389.07 358.15 2.70 c p3 -71.73 105.74 2.70 c p3 -135.20 147.81 2.70 c p3 -198.67 202.25 2.70 c p3 -262.13 249.27 2.70 c p3 -325.60 286.39 2.70 c p3 -389.07 343.30 2.70 c p3 -71.73 103.26 2.70 c p3 -135.20 145.33 2.70 c p3 -198.67 202.25 2.70 c p3 -262.13 249.27 2.70 c p3 -325.60 283.91 2.70 c p3 -389.07 353.20 2.70 c p3 -np -71.73 100.60 m -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -199.18 199.07 lineto -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -326.62 297.55 lineto -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -o -1 0 0 rgb -0.75 setlinewidth -[ 0.00 4.00 3.00 4.00] 0 setdash -np -71.73 84.53 m -1.28 1.00 l -1.27 0.99 l -1.28 0.99 l -1.27 1.00 l -1.28 0.99 l -1.27 0.99 l -1.27 0.99 l -1.28 1.00 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 1.00 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.27 0.99 l -1.28 1.00 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 1.00 l -1.28 0.99 l -1.27 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 1.00 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 1.00 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -199.18 183.52 lineto -1.27 0.98 l -1.28 0.99 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -326.62 281.82 lineto -1.27 0.98 l -1.28 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.97 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.97 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.97 l -1.27 0.98 l -1.28 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.98 l -1.28 0.97 l -1.27 0.98 l -1.28 0.98 l -o -np -71.73 116.66 m -1.28 0.98 l -1.27 0.98 l -1.28 0.97 l -1.27 0.98 l -1.28 0.98 l -1.27 0.97 l -1.27 0.98 l -1.28 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.98 l -1.28 0.97 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.97 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.97 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -199.18 214.63 lineto -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.99 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -326.62 313.28 lineto -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 1.00 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 1.00 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 1.00 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 1.00 l -1.27 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 1.00 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 1.00 l -1.27 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 1.00 l -1.28 0.99 l -1.27 0.99 l -1.28 1.00 l -o -0 0 0 rgb -0.75 setlinewidth -[] 0 setdash -np -71.73 100.60 m -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -199.18 199.07 lineto -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -326.62 297.55 lineto -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -o -0 0.5451 0 rgb -0.75 setlinewidth -[ 0.00 4.00] 0 setdash -np -71.73 95.65 m -1.28 1.01 l -1.27 1.01 l -1.28 1.02 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.02 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.27 1.00 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.00 l -1.28 1.01 l -1.27 1.01 l -1.27 1.01 l -1.28 1.00 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.00 l -1.27 1.01 l -1.27 1.01 l -1.28 1.00 l -1.27 1.01 l -1.28 1.00 l -1.27 1.01 l -1.28 1.01 l -1.27 1.00 l -1.28 1.01 l -1.27 1.00 l -1.27 1.01 l -1.28 1.00 l -1.27 1.01 l -1.28 1.00 l -1.27 1.00 l -1.28 1.01 l -1.27 1.00 l -1.28 1.00 l -1.27 1.01 l -1.27 1.00 l -1.28 1.00 l -1.27 1.01 l -1.28 1.00 l -1.27 1.00 l -1.28 1.00 l -1.27 1.00 l -1.28 1.01 l -1.27 1.00 l -1.27 1.00 l -1.28 1.00 l -1.27 1.00 l -1.28 1.00 l -1.27 1.00 l -1.28 1.00 l -1.27 1.00 l -1.28 1.00 l -1.27 1.00 l -1.27 0.99 l -1.28 1.00 l -1.27 1.00 l -1.28 1.00 l -1.27 0.99 l -1.28 1.00 l -1.27 1.00 l -1.28 0.99 l -1.27 1.00 l -1.27 1.00 l -1.28 0.99 l -1.27 1.00 l -1.28 0.99 l -1.27 1.00 l -199.18 196.17 lineto -1.27 0.99 l -1.28 1.00 l -1.27 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 1.00 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.97 l -1.27 0.98 l -1.28 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.97 l -1.27 0.97 l -1.28 0.98 l -1.27 0.97 l -1.28 0.97 l -1.27 0.98 l -1.28 0.97 l -1.27 0.97 l -1.28 0.97 l -1.27 0.97 l -1.27 0.97 l -1.28 0.97 l -1.27 0.97 l -1.28 0.97 l -1.27 0.97 l -1.28 0.97 l -1.27 0.97 l -1.28 0.97 l -1.27 0.97 l -1.27 0.97 l -1.28 0.96 l -1.27 0.97 l -1.28 0.97 l -1.27 0.97 l -1.28 0.97 l -1.27 0.96 l -1.28 0.97 l -1.27 0.97 l -1.27 0.96 l -1.28 0.97 l -1.27 0.96 l -1.28 0.97 l -1.27 0.96 l -1.28 0.97 l -1.27 0.96 l -1.28 0.97 l -1.27 0.96 l -1.27 0.97 l -1.28 0.96 l -1.27 0.97 l -1.28 0.96 l -1.27 0.96 l -1.28 0.97 l -1.27 0.96 l -1.28 0.96 l -326.62 293.82 lineto -1.27 0.96 l -1.28 0.96 l -1.27 0.97 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.27 0.97 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.97 l -1.27 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.95 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.95 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.27 0.96 l -1.28 0.95 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -o -np -71.73 105.54 m -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.95 l -1.27 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.95 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.95 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.27 0.96 l -1.28 0.97 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.97 l -1.27 0.96 l -1.28 0.96 l -1.27 0.96 l -1.28 0.96 l -1.27 0.97 l -1.28 0.96 l -1.27 0.96 l -1.27 0.97 l -1.28 0.96 l -1.27 0.96 l -1.28 0.97 l -1.27 0.96 l -1.28 0.96 l -1.27 0.97 l -1.28 0.96 l -1.27 0.97 l -1.27 0.96 l -1.28 0.97 l -1.27 0.96 l -1.28 0.97 l -1.27 0.96 l -1.28 0.97 l -1.27 0.96 l -1.28 0.97 l -1.27 0.96 l -1.27 0.97 l -1.28 0.97 l -1.27 0.96 l -1.28 0.97 l -1.27 0.97 l -1.28 0.97 l -1.27 0.97 l -1.28 0.96 l -1.27 0.97 l -1.27 0.97 l -1.28 0.97 l -1.27 0.97 l -1.28 0.97 l -1.27 0.97 l -1.28 0.97 l -1.27 0.97 l -1.28 0.97 l -1.27 0.97 l -1.27 0.97 l -1.28 0.97 l -1.27 0.97 l -1.28 0.97 l -1.27 0.98 l -1.28 0.97 l -1.27 0.97 l -1.28 0.98 l -1.27 0.97 l -1.27 0.97 l -1.28 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.97 l -199.18 201.98 lineto -1.27 0.97 l -1.28 0.98 l -1.27 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.97 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.27 0.98 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.27 0.99 l -1.28 0.98 l -1.27 0.99 l -1.28 0.99 l -1.27 0.98 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 0.99 l -1.28 0.99 l -1.27 1.00 l -1.28 0.99 l -1.27 0.99 l -1.27 0.99 l -1.28 1.00 l -1.27 0.99 l -1.28 0.99 l -1.27 1.00 l -1.28 0.99 l -1.27 1.00 l -1.28 0.99 l -1.27 1.00 l -1.27 1.00 l -1.28 0.99 l -1.27 1.00 l -1.28 1.00 l -1.27 0.99 l -1.28 1.00 l -1.27 1.00 l -1.28 1.00 l -1.27 0.99 l -1.27 1.00 l -1.28 1.00 l -1.27 1.00 l -1.28 1.00 l -1.27 1.00 l -1.28 1.00 l -1.27 1.00 l -1.28 1.00 l -1.27 1.00 l -1.27 1.00 l -1.28 1.01 l -1.27 1.00 l -1.28 1.00 l -1.27 1.00 l -1.28 1.00 l -1.27 1.01 l -1.28 1.00 l -1.27 1.00 l -1.27 1.01 l -1.28 1.00 l -1.27 1.00 l -1.28 1.01 l -1.27 1.00 l -1.28 1.00 l -1.27 1.01 l -1.28 1.00 l -1.27 1.01 l -1.27 1.00 l -1.28 1.01 l -1.27 1.00 l -1.28 1.01 l -1.27 1.01 l -1.28 1.00 l -1.27 1.01 l -1.28 1.00 l -326.62 301.28 lineto -1.27 1.01 l -1.28 1.00 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.00 l -1.27 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.00 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.00 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.02 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.27 1.01 l -1.28 1.01 l -1.27 1.01 l -1.28 1.02 l -1.27 1.01 l -1.28 1.01 l -o -/bg { 0.9490 0.9490 0.9490 } def -0 0 0 rgb -0.75 setlinewidth -[] 0 setdash -71.73 361.87 133.62 -51.84 r p3 -1.50 setlinewidth -[] 0 setdash -np -81.45 348.91 m -19.44 0 l -o -0 0.5451 0 rgb -1.50 setlinewidth -[ 0.00 6.00] 0 setdash -np -81.45 335.95 m -19.44 0 l -o -1 0 0 rgb -1.50 setlinewidth -[ 0.00 6.00 4.50 6.00] 0 setdash -np -81.45 322.99 m -19.44 0 l -o -/ps 11 def /Font1 findfont 11 s -0 0 0 rgb -110.61 344.96 (Fitted Line) 0 0 0 t -110.61 332.00 (Confidence Bands) 0 0 0 t -110.61 319.04 (Prediction Bands) 0 0 0 t -ep -%%Trailer -%%Pages: 1 -%%EOF diff --git a/trunk/chemCal/vignettes/chemCal-001.pdf b/trunk/chemCal/vignettes/chemCal-001.pdf Binary files differdeleted file mode 100644 index a8c2629..0000000 --- a/trunk/chemCal/vignettes/chemCal-001.pdf +++ /dev/null diff --git a/trunk/chemCal/vignettes/chemCal-002.eps b/trunk/chemCal/vignettes/chemCal-002.eps deleted file mode 100644 index c0ee749..0000000 --- a/trunk/chemCal/vignettes/chemCal-002.eps +++ /dev/null @@ -1,288 +0,0 @@ -%!PS-Adobe-3.0 -%%DocumentNeededResources: font Helvetica -%%+ font Helvetica-Bold -%%+ font Helvetica-Oblique -%%+ font Helvetica-BoldOblique -%%+ font Symbol -%%DocumentMedia: special 432 432 0 () () -%%Title: R Graphics Output -%%Creator: R Software -%%Pages: (atend) -%%BoundingBox: 0 0 432 432 -%%EndComments -%%BeginProlog -/bp { gs gs } def -% begin .ps.prolog -/gs { gsave } def -/gr { grestore } def -/ep { showpage gr gr } def -/m { moveto } def -/l { rlineto } def -/np { newpath } def -/cp { closepath } def -/f { fill } def -/o { stroke } def -/c { newpath 0 360 arc } def -/r { 4 2 roll moveto 1 copy 3 -1 roll exch 0 exch rlineto 0 rlineto -1 mul 0 exch rlineto closepath } def -/p1 { stroke } def -/p2 { gsave bg setrgbcolor fill grestore newpath } def -/p3 { gsave bg setrgbcolor fill grestore stroke } def -/t { 6 -2 roll moveto gsave rotate - ps mul neg 0 2 1 roll rmoveto - 1 index stringwidth pop - mul neg 0 rmoveto show grestore } def -/cl { grestore gsave newpath 3 index 3 index moveto 1 index - 4 -1 roll lineto exch 1 index lineto lineto - closepath clip newpath } def -/rgb { setrgbcolor } def -/s { scalefont setfont } def -% end .ps.prolog -%%IncludeResource: font Helvetica -/Helvetica findfont -dup length dict begin - {1 index /FID ne {def} {pop pop} ifelse} forall - /Encoding ISOLatin1Encoding def - currentdict - end -/Font1 exch definefont pop -%%IncludeResource: font Helvetica-Bold -/Helvetica-Bold findfont -dup length dict begin - {1 index /FID ne {def} {pop pop} ifelse} forall - /Encoding ISOLatin1Encoding def - currentdict - end -/Font2 exch definefont pop -%%IncludeResource: font Helvetica-Oblique -/Helvetica-Oblique findfont -dup length dict begin - {1 index /FID ne {def} {pop pop} ifelse} forall - /Encoding ISOLatin1Encoding def - currentdict - end -/Font3 exch definefont pop -%%IncludeResource: font Helvetica-BoldOblique -/Helvetica-BoldOblique findfont -dup length dict begin - {1 index /FID ne {def} {pop pop} ifelse} forall - /Encoding ISOLatin1Encoding def - currentdict - end -/Font4 exch definefont pop -%%IncludeResource: font Symbol -/Symbol findfont -dup length dict begin - {1 index /FID ne {def} {pop pop} ifelse} forall - currentdict - end -/Font5 exch definefont pop -%%EndProlog -%%Page: 1 1 -bp -0.00 0.00 432.00 432.00 cl -0 0 0 rgb -0.75 setlinewidth -[] 0 setdash -1 setlinecap -1 setlinejoin -10.00 setmiterlimit -np -62.37 73.44 m -320.26 0 l -o -np -62.37 73.44 m -0 -7.20 l -o -np -126.42 73.44 m -0 -7.20 l -o -np -190.47 73.44 m -0 -7.20 l -o -np -254.53 73.44 m -0 -7.20 l -o -np -318.58 73.44 m -0 -7.20 l -o -np -382.63 73.44 m -0 -7.20 l -o -/ps 12 def /Font1 findfont 12 s -62.37 47.52 (0) .5 0 0 t -126.42 47.52 (20) .5 0 0 t -190.47 47.52 (40) .5 0 0 t -254.53 47.52 (60) .5 0 0 t -318.58 47.52 (80) .5 0 0 t -382.63 47.52 (100) .5 0 0 t -np -59.04 84.53 m -0 265.28 l -o -np -59.04 84.53 m --7.20 0 l -o -np -59.04 172.96 m --7.20 0 l -o -np -59.04 261.38 m --7.20 0 l -o -np -59.04 349.81 m --7.20 0 l -o -41.76 84.53 (0.0) .5 0 90 t -41.76 172.96 (0.5) .5 0 90 t -41.76 261.38 (1.0) .5 0 90 t -41.76 349.81 (1.5) .5 0 90 t -np -59.04 73.44 m -342.72 0 l -0 299.52 l --342.72 0 l -0 -299.52 l -o -0.00 0.00 432.00 432.00 cl -/ps 12 def /Font1 findfont 12 s -0 0 0 rgb -230.40 18.72 (Fitted values) .5 0 0 t -0.75 setlinewidth -[] 0 setdash -1 setlinecap -1 setlinejoin -10.00 setmiterlimit -np -12.17 157.30 m --0.87 1.55 l -4.54 1.55 l --11.58 2.07 l -0 126.08 l -o -np -15.84 164.28 m --9.07 0 l -o -15.61 164.83 (S) 0 0 90 t -15.61 172.84 (t) 0 0 90 t -15.61 176.17 (a) 0 0 90 t -15.61 182.84 (n) 0 0 90 t -15.61 189.52 (d) 0 0 90 t -15.61 196.19 (a) 0 0 90 t -15.61 202.86 (r) 0 0 90 t -15.61 206.86 (d) 0 0 90 t -15.61 213.53 (i) 0 0 90 t -15.61 216.19 (z) 0 0 90 t -15.61 222.19 (e) 0 0 90 t -15.61 228.86 (d) 0 0 90 t -15.61 235.54 ( ) 0 0 90 t -15.61 238.87 (r) 0 0 90 t -15.61 242.87 (e) 0 0 90 t -15.61 249.54 (s) 0 0 90 t -15.61 255.54 (i) 0 0 90 t -15.61 258.20 (d) 0 0 90 t -15.61 264.88 (u) 0 0 90 t -15.61 271.55 (a) 0 0 90 t -15.61 278.22 (l) 0 0 90 t -15.61 280.88 (s) 0 0 90 t -np -15.84 287.44 m --9.07 0 l -o -59.04 73.44 401.76 372.96 cl -0 0 0 rgb -0.75 setlinewidth -[] 0 setdash -1 setlinecap -1 setlinejoin -10.00 setmiterlimit -71.73 193.16 2.70 c p1 -135.20 173.55 2.70 c p1 -198.67 207.97 2.70 c p1 -262.13 242.95 2.70 c p1 -325.60 361.87 2.70 c p1 -389.07 232.26 2.70 c p1 -71.73 113.44 2.70 c p1 -135.20 255.74 2.70 c p1 -198.67 275.24 2.70 c p1 -262.13 165.78 2.70 c p1 -325.60 197.45 2.70 c p1 -389.07 361.37 2.70 c p1 -71.73 193.16 2.70 c p1 -135.20 220.98 2.70 c p1 -198.67 245.16 2.70 c p1 -262.13 242.95 2.70 c p1 -325.60 269.30 2.70 c p1 -389.07 318.44 2.70 c p1 -71.73 235.40 2.70 c p1 -135.20 173.55 2.70 c p1 -198.67 207.97 2.70 c p1 -262.13 165.78 2.70 c p1 -325.60 296.27 2.70 c p1 -389.07 189.74 2.70 c p1 -71.73 193.16 2.70 c p1 -135.20 220.98 2.70 c p1 -198.67 207.97 2.70 c p1 -262.13 165.78 2.70 c p1 -325.60 320.17 2.70 c p1 -389.07 265.60 2.70 c p1 -1 0 0 rgb -np -71.73 188.63 m -0 0 l -0 0 l -0 0 l -0 0 l -63.47 19.56 l -0 0 l -0 0 l -0 0 l -0 0 l -63.47 6.32 l -0 0 l -0 0 l -0 0 l -0 0 l -63.46 17.52 l -0 0 l -0 0 l -0 0 l -0 0 l -63.47 25.69 l -0 0 l -0 0 l -0 0 l -0 0 l -63.47 30.70 l -0 0 l -0 0 l -0 0 l -0 0 l -o -0.00 0.00 432.00 432.00 cl -/ps 12 def /Font1 findfont 12 s -0 0 0 rgb -230.40 4.32 (lm\(y ~ x\)) .5 0 0 t -0.00 0.00 432.00 432.00 cl -/ps 12 def /Font1 findfont 12 s -0 0 0 rgb -230.40 379.44 (Scale-Location) .5 0 0 t -0.00 0.00 432.00 432.00 cl -/ps 9 def /Font1 findfont 9 s -0 0 0 rgb -329.20 359.71 (5) 0 0 0 t -385.47 359.22 (12) 1 0 0 t -329.20 318.02 (29) 0 0 0 t -ep -%%Trailer -%%Pages: 1 -%%EOF diff --git a/trunk/chemCal/vignettes/chemCal-002.pdf b/trunk/chemCal/vignettes/chemCal-002.pdf Binary files differdeleted file mode 100644 index 86bf0ac..0000000 --- a/trunk/chemCal/vignettes/chemCal-002.pdf +++ /dev/null diff --git a/trunk/chemCal/vignettes/chemCal.Rnw b/trunk/chemCal/vignettes/chemCal.Rnw deleted file mode 100644 index 3c46540..0000000 --- a/trunk/chemCal/vignettes/chemCal.Rnw +++ /dev/null @@ -1,161 +0,0 @@ -\documentclass[a4paper]{article} -%\VignetteIndexEntry{Short manual for the chemCal package} -\usepackage{hyperref} - -\title{Basic calibration functions for analytical chemistry} -\author{Johannes Ranke} - -\begin{document} -\maketitle - -The \texttt{chemCal} package was first designed in the course of a lecture and lab -course on "analytics of organic trace contaminants" at the University of Bremen -from October to December 2004. In the fall 2005, an email exchange with -Ron Wehrens led to the belief that it would be desirable to implement the -inverse prediction method given in \cite{massart97} since it also covers the -case of weighted regression. Studies of the IUPAC orange book and of DIN 32645 -as well as publications by Currie and the Analytical Method Committee of the -Royal Society of Chemistry and a nice paper by Castillo and Castells provided -further understanding of the matter. - -At the moment, the package consists of four functions, working on univariate -linear models of class \texttt{lm} or \texttt{rlm}, plus to datasets for -validation. - -A \href{http://bugs.r-project.org/bugzilla3/show_bug.cgi?id=8877}{bug -report (PR\#8877)} and the following e-mail exchange on the r-devel mailing list about -prediction intervals from weighted regression entailed some further studies -on this subject. However, I did not encounter any proof or explanation of the -formula cited below yet, so I can't really confirm that Massart's method is correct. - -When calibrating an analytical method, the first task is to generate a suitable -model. If we want to use the \texttt{chemCal} functions, we will have to -restrict ourselves to univariate, possibly weighted, linear regression so far. - -Once such a model has been created, the calibration can be graphically -shown by using the \texttt{calplot} function: - -<<echo=TRUE,fig=TRUE>>= -library(chemCal) -data(massart97ex3) -m0 <- lm(y ~ x, data = massart97ex3) -calplot(m0) -@ - -As we can see, the scatter increases with increasing x. This is also -illustrated by one of the diagnostic plots for linear models -provided by R: - -<<echo=TRUE,fig=TRUE>>= -plot(m0,which=3) -@ - -Therefore, in Example 8 in \cite{massart97} weighted regression -is proposed which can be reproduced by - -<<>>= -attach(massart97ex3) -yx <- split(y, x) -ybar <- sapply(yx, mean) -s <- round(sapply(yx, sd), digits = 2) -w <- round(1 / (s^2), digits = 3) -weights <- w[factor(x)] -m <- lm(y ~ x, w = weights) -@ - -If we now want to predict a new x value from measured y values, -we use the \texttt{inverse.predict} function: - -<<>>= -inverse.predict(m, 15, ws=1.67) -inverse.predict(m, 90, ws = 0.145) -@ - -The weight \texttt{ws} assigned to the measured y value has to be -given by the user in the case of weighted regression, or alternatively, -the approximate variance \texttt{var.s} at this location. - -\section*{Theory for \texttt{inverse.predict}} -Equation 8.28 in \cite{massart97} gives a general equation for predicting the -standard error $s_{\hat{x_s}}$ for an x value predicted from measurements of y -according to the linear calibration function $ y = b_0 + b_1 \cdot x$: - -\begin{equation} -s_{\hat{x_s}} = \frac{s_e}{b_1} \sqrt{\frac{1}{w_s m} + \frac{1}{\sum{w_i}} + - \frac{(\bar{y_s} - \bar{y_w})^2 \sum{w_i}} - {{b_1}^2 \left( \sum{w_i} \sum{w_i {x_i}^2} - - {\left( \sum{ w_i x_i } \right)}^2 \right) }} -\end{equation} - -with - -\begin{equation} -s_e = \sqrt{ \frac{\sum w_i (y_i - \hat{y_i})^2}{n - 2}} -\end{equation} - -where $w_i$ is the weight for calibration standard $i$, $y_i$ is the mean $y$ -value (!) observed for standard $i$, $\hat{y_i}$ is the estimated value for -standard $i$, $n$ is the number calibration standards, $w_s$ is the weight -attributed to the sample $s$, $m$ is the number of replicate measurements of -sample $s$, $\bar{y_s}$ is the mean response for the sample, -$\bar{y_w} = \frac{\sum{w_i y_i}}{\sum{w_i}}$ is the weighted mean of responses -$y_i$, and $x_i$ is the given $x$ value for standard $i$. - -The weight $w_s$ for the sample should be estimated or calculated in accordance -to the weights used in the linear regression. - -I adjusted the above equation in order to be able to take a different -precisions in standards and samples into account. In analogy to Equation 8.26 -from \cite{massart97} we get - -\begin{equation} -s_{\hat{x_s}} = \frac{1}{b_1} \sqrt{\frac{{s_s}^2}{w_s m} + - {s_e}^2 \left( \frac{1}{\sum{w_i}} + - \frac{(\bar{y_s} - \bar{y_w})^2 \sum{w_i}} - {{b_1}^2 \left( \sum{w_i} \sum{w_i {x_i}^2} - {\left( \sum{ w_i x_i } \right)}^2 \right) } \right) } -\end{equation} - -where I interpret $\frac{{s_s}^2}{w_s}$ as an estimator of the variance at location -$\hat{x_s}$, which can be replaced by a user-specified value using the argument -\texttt{var.s} of the function \texttt{inverse.predict}. - -\section*{Fitting and using a variance function} - -In the R package \texttt{nlme} variance functions are used for weighted -regressions. But since the \texttt{predict.nlme} method does not calculate -prediction intervals, this is not useful for the \texttt{calplot} function. - -Two approaches could be used for fitting variance functions, one based on -residuals from an unweighted fit, and one based on just the variances -of the different samples along the x axis. If we used the residuals for -fitting, a bias of the model in a certain area would result in a higher -variance, so it seems preferable to choose the second approach. Of course, -a prerequisite is to have sufficient repetitions for each sample in any -case. - -Let's take the above example and estimate a variance function - -<<>>= -library(doBy) -massart97ex3 -massart97ex3$x <- factor(massart97ex3$x) -summary <- summaryBy(y~x, data = massart97ex3,FUN=c(mean,sd,var)) -summary$x <- as.numeric(as.vector((summary$x))) -plot(summary$x, summary$y.var, - xlim=c(0,50), - ylim=c(0,max(summary$y.var))) -varModel <- lm(y.var ~ I(x^2) + x, data=summary) -varCurve <- predict(varModel, newdata=data.frame(x=0:5000/100)) -lines(0:5000/100,varCurve) -@ - - -\begin{thebibliography}{1} -\bibitem{massart97} -Massart, L.M, Vandenginste, B.G.M., Buydens, L.M.C., De Jong, S., Lewi, P.J., -Smeyers-Verbeke, J. -\newblock Handbook of Chemometrics and Qualimetrics: Part A, -\newblock Elsevier, Amsterdam, 1997 -\end{thebibliography} - -\end{document} diff --git a/trunk/chemCal/vignettes/chemCal.aux b/trunk/chemCal/vignettes/chemCal.aux deleted file mode 100644 index 20bfc98..0000000 --- a/trunk/chemCal/vignettes/chemCal.aux +++ /dev/null @@ -1,18 +0,0 @@ -\relax -\ifx\hyper@anchor\@undefined -\global \let \oldcontentsline\contentsline -\gdef \contentsline#1#2#3#4{\oldcontentsline{#1}{#2}{#3}} -\global \let \oldnewlabel\newlabel -\gdef \newlabel#1#2{\newlabelxx{#1}#2} -\gdef \newlabelxx#1#2#3#4#5#6{\oldnewlabel{#1}{{#2}{#3}}} -\AtEndDocument{\let \contentsline\oldcontentsline -\let \newlabel\oldnewlabel} -\else -\global \let \hyper@last\relax -\fi - -\citation{massart97} -\citation{massart97} -\citation{massart97} -\citation{massart97} -\bibcite{massart97}{1} diff --git a/trunk/chemCal/vignettes/chemCal.bbl b/trunk/chemCal/vignettes/chemCal.bbl deleted file mode 100644 index e69de29..0000000 --- a/trunk/chemCal/vignettes/chemCal.bbl +++ /dev/null diff --git a/trunk/chemCal/vignettes/chemCal.blg b/trunk/chemCal/vignettes/chemCal.blg deleted file mode 100644 index 2904a63..0000000 --- a/trunk/chemCal/vignettes/chemCal.blg +++ /dev/null @@ -1,46 +0,0 @@ -This is BibTeX, Version 0.99c (TeX Live 2009/Debian) -The top-level auxiliary file: chemCal.aux -I found no \bibdata command---while reading file chemCal.aux -I found no \bibstyle command---while reading file chemCal.aux -You've used 1 entry, - 0 wiz_defined-function locations, - 84 strings with 497 characters, -and the built_in function-call counts, 0 in all, are: -= -- 0 -> -- 0 -< -- 0 -+ -- 0 -- -- 0 -* -- 0 -:= -- 0 -add.period$ -- 0 -call.type$ -- 0 -change.case$ -- 0 -chr.to.int$ -- 0 -cite$ -- 0 -duplicate$ -- 0 -empty$ -- 0 -format.name$ -- 0 -if$ -- 0 -int.to.chr$ -- 0 -int.to.str$ -- 0 -missing$ -- 0 -newline$ -- 0 -num.names$ -- 0 -pop$ -- 0 -preamble$ -- 0 -purify$ -- 0 -quote$ -- 0 -skip$ -- 0 -stack$ -- 0 -substring$ -- 0 -swap$ -- 0 -text.length$ -- 0 -text.prefix$ -- 0 -top$ -- 0 -type$ -- 0 -warning$ -- 0 -while$ -- 0 -width$ -- 0 -write$ -- 0 -(There were 2 error messages) diff --git a/trunk/chemCal/vignettes/chemCal.log b/trunk/chemCal/vignettes/chemCal.log deleted file mode 100644 index 2720411..0000000 --- a/trunk/chemCal/vignettes/chemCal.log +++ /dev/null @@ -1,399 +0,0 @@ -This is pdfTeX, Version 3.1415926-1.40.10 (TeX Live 2009/Debian) (format=pdflatex 2012.6.4) 11 SEP 2012 23:19 -entering extended mode - %&-line parsing enabled. -**chemCal.tex -(./chemCal.tex -LaTeX2e <2009/09/24> -Babel <v3.8l> and hyphenation patterns for english, usenglishmax, dumylang, noh -yphenation, ngerman, german, german-x-2009-06-19, ngerman-x-2009-06-19, loaded. - -(/usr/share/texmf-texlive/tex/latex/base/article.cls -Document Class: article 2007/10/19 v1.4h Standard LaTeX document class -(/usr/share/texmf-texlive/tex/latex/base/size10.clo -File: size10.clo 2007/10/19 v1.4h Standard LaTeX file (size option) -) -\c@part=\count79 -\c@section=\count80 -\c@subsection=\count81 -\c@subsubsection=\count82 -\c@paragraph=\count83 -\c@subparagraph=\count84 -\c@figure=\count85 -\c@table=\count86 -\abovecaptionskip=\skip41 -\belowcaptionskip=\skip42 -\bibindent=\dimen102 -) -(/usr/share/texmf-texlive/tex/latex/hyperref/hyperref.sty -Package: hyperref 2009/10/09 v6.79a Hypertext links for LaTeX - -(/usr/share/texmf-texlive/tex/latex/graphics/keyval.sty -Package: keyval 1999/03/16 v1.13 key=value parser (DPC) -\KV@toks@=\toks14 -) -(/usr/share/texmf-texlive/tex/generic/oberdiek/ifpdf.sty -Package: ifpdf 2009/04/10 v2.0 Provides the ifpdf switch (HO) -Package ifpdf Info: pdfTeX in pdf mode detected. -) -(/usr/share/texmf-texlive/tex/generic/oberdiek/ifvtex.sty -Package: ifvtex 2008/11/04 v1.4 Switches for detecting VTeX and its modes (HO) -Package ifvtex Info: VTeX not detected. -) -(/usr/share/texmf-texlive/tex/generic/ifxetex/ifxetex.sty -Package: ifxetex 2009/01/23 v0.5 Provides ifxetex conditional -) -(/usr/share/texmf-texlive/tex/latex/oberdiek/hycolor.sty -Package: hycolor 2009/10/02 v1.5 Code for color options of hyperref/bookmark (H -O) - -(/usr/share/texmf-texlive/tex/latex/oberdiek/xcolor-patch.sty -Package: xcolor-patch 2009/10/02 xcolor patch -)) -\@linkdim=\dimen103 -\Hy@linkcounter=\count87 -\Hy@pagecounter=\count88 - -(/usr/share/texmf-texlive/tex/latex/hyperref/pd1enc.def -File: pd1enc.def 2009/10/09 v6.79a Hyperref: PDFDocEncoding definition (HO) -) -(/usr/share/texmf-texlive/tex/generic/oberdiek/etexcmds.sty -Package: etexcmds 2007/12/12 v1.2 Prefix for e-TeX command names (HO) - -(/usr/share/texmf-texlive/tex/generic/oberdiek/infwarerr.sty -Package: infwarerr 2007/09/09 v1.2 Providing info/warning/message (HO) -) -Package etexcmds Info: Could not find \expanded. -(etexcmds) That can mean that you are not using pdfTeX 1.50 or -(etexcmds) that some package has redefined \expanded. -(etexcmds) In the latter case, load this package earlier. -) -(/usr/share/texmf-texlive/tex/latex/latexconfig/hyperref.cfg -File: hyperref.cfg 2002/06/06 v1.2 hyperref configuration of TeXLive -) -(/usr/share/texmf-texlive/tex/latex/oberdiek/kvoptions.sty -Package: kvoptions 2009/08/13 v3.4 Keyval support for LaTeX options (HO) - -(/usr/share/texmf-texlive/tex/generic/oberdiek/kvsetkeys.sty -Package: kvsetkeys 2009/07/30 v1.5 Key value parser with default handler suppor -t (HO) -)) -Package hyperref Info: Hyper figures OFF on input line 2975. -Package hyperref Info: Link nesting OFF on input line 2980. -Package hyperref Info: Hyper index ON on input line 2983. -Package hyperref Info: Plain pages OFF on input line 2990. -Package hyperref Info: Backreferencing OFF on input line 2995. - -Implicit mode ON; LaTeX internals redefined -Package hyperref Info: Bookmarks ON on input line 3191. -(/usr/share/texmf-texlive/tex/latex/ltxmisc/url.sty -\Urlmuskip=\muskip10 -Package: url 2006/04/12 ver 3.3 Verb mode for urls, etc. -) -LaTeX Info: Redefining \url on input line 3428. - -(/usr/share/texmf-texlive/tex/generic/oberdiek/bitset.sty -Package: bitset 2007/09/28 v1.0 Data type bit set (HO) - -(/usr/share/texmf-texlive/tex/generic/oberdiek/intcalc.sty -Package: intcalc 2007/09/27 v1.1 Expandable integer calculations (HO) -) -(/usr/share/texmf-texlive/tex/generic/oberdiek/bigintcalc.sty -Package: bigintcalc 2007/11/11 v1.1 Expandable big integer calculations (HO) - -(/usr/share/texmf-texlive/tex/generic/oberdiek/pdftexcmds.sty -Package: pdftexcmds 2009/09/23 v0.6 LuaTeX support for pdfTeX utility functions - (HO) - -(/usr/share/texmf-texlive/tex/generic/oberdiek/ifluatex.sty -Package: ifluatex 2009/04/17 v1.2 Provides the ifluatex switch (HO) -Package ifluatex Info: LuaTeX not detected. -) -(/usr/share/texmf-texlive/tex/generic/oberdiek/ltxcmds.sty -Package: ltxcmds 2009/08/05 v1.0 Some LaTeX kernel commands for general use (HO -) -) -Package pdftexcmds Info: LuaTeX not detected. -Package pdftexcmds Info: \pdf@primitive is available. -Package pdftexcmds Info: \pdf@ifprimitive is available. -))) -\Fld@menulength=\count89 -\Field@Width=\dimen104 -\Fld@charsize=\dimen105 -\Field@toks=\toks15 -Package hyperref Info: Hyper figures OFF on input line 4377. -Package hyperref Info: Link nesting OFF on input line 4382. -Package hyperref Info: Hyper index ON on input line 4385. -Package hyperref Info: backreferencing OFF on input line 4392. -Package hyperref Info: Link coloring OFF on input line 4397. -Package hyperref Info: Link coloring with OCG OFF on input line 4402. -Package hyperref Info: PDF/A mode OFF on input line 4407. - -(/usr/share/texmf-texlive/tex/generic/oberdiek/atbegshi.sty -Package: atbegshi 2008/07/31 v1.9 At begin shipout hook (HO) -) -\Hy@abspage=\count90 -\c@Item=\count91 -\c@Hfootnote=\count92 -) -*hyperref using default driver hpdftex* -(/usr/share/texmf-texlive/tex/latex/hyperref/hpdftex.def -File: hpdftex.def 2009/10/09 v6.79a Hyperref driver for pdfTeX -\Fld@listcount=\count93 -) -(/usr/share/texmf/tex/latex/R/tex/latex/Sweave.sty -Package: Sweave - -(/usr/share/texmf-texlive/tex/latex/base/ifthen.sty -Package: ifthen 2001/05/26 v1.1c Standard LaTeX ifthen package (DPC) -) -(/usr/share/texmf-texlive/tex/latex/graphics/graphicx.sty -Package: graphicx 1999/02/16 v1.0f Enhanced LaTeX Graphics (DPC,SPQR) - -(/usr/share/texmf-texlive/tex/latex/graphics/graphics.sty -Package: graphics 2009/02/05 v1.0o Standard LaTeX Graphics (DPC,SPQR) - -(/usr/share/texmf-texlive/tex/latex/graphics/trig.sty -Package: trig 1999/03/16 v1.09 sin cos tan (DPC) -) -(/etc/texmf/tex/latex/config/graphics.cfg -File: graphics.cfg 2009/08/28 v1.8 graphics configuration of TeX Live -) -Package graphics Info: Driver file: pdftex.def on input line 91. - -(/usr/share/texmf-texlive/tex/latex/pdftex-def/pdftex.def -File: pdftex.def 2010/03/12 v0.04p Graphics/color for pdfTeX -\Gread@gobject=\count94 -)) -\Gin@req@height=\dimen106 -\Gin@req@width=\dimen107 -) -(/usr/share/texmf-texlive/tex/latex/fancyvrb/fancyvrb.sty -Package: fancyvrb 2008/02/07 - -Style option: `fancyvrb' v2.7a, with DG/SPQR fixes, and firstline=lastline fix -<2008/02/07> (tvz) -\FV@CodeLineNo=\count95 -\FV@InFile=\read1 -\FV@TabBox=\box26 -\c@FancyVerbLine=\count96 -\FV@StepNumber=\count97 -\FV@OutFile=\write3 -) (/usr/share/texmf/tex/latex/R/tex/latex/upquote.sty -Package: upquote 2003/08/11 v1.1 Covington's upright-quote modification to verb -atim and verb - -(/usr/share/texmf-texlive/tex/latex/base/textcomp.sty -Package: textcomp 2005/09/27 v1.99g Standard LaTeX package -Package textcomp Info: Sub-encoding information: -(textcomp) 5 = only ISO-Adobe without \textcurrency -(textcomp) 4 = 5 + \texteuro -(textcomp) 3 = 4 + \textohm -(textcomp) 2 = 3 + \textestimated + \textcurrency -(textcomp) 1 = TS1 - \textcircled - \t -(textcomp) 0 = TS1 (full) -(textcomp) Font families with sub-encoding setting implement -(textcomp) only a restricted character set as indicated. -(textcomp) Family '?' is the default used for unknown fonts. -(textcomp) See the documentation for details. -Package textcomp Info: Setting ? sub-encoding to TS1/1 on input line 71. - -(/usr/share/texmf-texlive/tex/latex/base/ts1enc.def -File: ts1enc.def 2001/06/05 v3.0e (jk/car/fm) Standard LaTeX file -) -LaTeX Info: Redefining \oldstylenums on input line 266. -Package textcomp Info: Setting cmr sub-encoding to TS1/0 on input line 281. -Package textcomp Info: Setting cmss sub-encoding to TS1/0 on input line 282. -Package textcomp Info: Setting cmtt sub-encoding to TS1/0 on input line 283. -Package textcomp Info: Setting cmvtt sub-encoding to TS1/0 on input line 284. -Package textcomp Info: Setting cmbr sub-encoding to TS1/0 on input line 285. -Package textcomp Info: Setting cmtl sub-encoding to TS1/0 on input line 286. -Package textcomp Info: Setting ccr sub-encoding to TS1/0 on input line 287. -Package textcomp Info: Setting ptm sub-encoding to TS1/4 on input line 288. -Package textcomp Info: Setting pcr sub-encoding to TS1/4 on input line 289. -Package textcomp Info: Setting phv sub-encoding to TS1/4 on input line 290. -Package textcomp Info: Setting ppl sub-encoding to TS1/3 on input line 291. -Package textcomp Info: Setting pag sub-encoding to TS1/4 on input line 292. -Package textcomp Info: Setting pbk sub-encoding to TS1/4 on input line 293. -Package textcomp Info: Setting pnc sub-encoding to TS1/4 on input line 294. -Package textcomp Info: Setting pzc sub-encoding to TS1/4 on input line 295. -Package textcomp Info: Setting bch sub-encoding to TS1/4 on input line 296. -Package textcomp Info: Setting put sub-encoding to TS1/5 on input line 297. -Package textcomp Info: Setting uag sub-encoding to TS1/5 on input line 298. -Package textcomp Info: Setting ugq sub-encoding to TS1/5 on input line 299. -Package textcomp Info: Setting ul8 sub-encoding to TS1/4 on input line 300. -Package textcomp Info: Setting ul9 sub-encoding to TS1/4 on input line 301. -Package textcomp Info: Setting augie sub-encoding to TS1/5 on input line 302. -Package textcomp Info: Setting dayrom sub-encoding to TS1/3 on input line 303. -Package textcomp Info: Setting dayroms sub-encoding to TS1/3 on input line 304. - -Package textcomp Info: Setting pxr sub-encoding to TS1/0 on input line 305. -Package textcomp Info: Setting pxss sub-encoding to TS1/0 on input line 306. -Package textcomp Info: Setting pxtt sub-encoding to TS1/0 on input line 307. -Package textcomp Info: Setting txr sub-encoding to TS1/0 on input line 308. -Package textcomp Info: Setting txss sub-encoding to TS1/0 on input line 309. -Package textcomp Info: Setting txtt sub-encoding to TS1/0 on input line 310. -Package textcomp Info: Setting futs sub-encoding to TS1/4 on input line 311. -Package textcomp Info: Setting futx sub-encoding to TS1/4 on input line 312. -Package textcomp Info: Setting futj sub-encoding to TS1/4 on input line 313. -Package textcomp Info: Setting hlh sub-encoding to TS1/3 on input line 314. -Package textcomp Info: Setting hls sub-encoding to TS1/3 on input line 315. -Package textcomp Info: Setting hlst sub-encoding to TS1/3 on input line 316. -Package textcomp Info: Setting hlct sub-encoding to TS1/5 on input line 317. -Package textcomp Info: Setting hlx sub-encoding to TS1/5 on input line 318. -Package textcomp Info: Setting hlce sub-encoding to TS1/5 on input line 319. -Package textcomp Info: Setting hlcn sub-encoding to TS1/5 on input line 320. -Package textcomp Info: Setting hlcw sub-encoding to TS1/5 on input line 321. -Package textcomp Info: Setting hlcf sub-encoding to TS1/5 on input line 322. -Package textcomp Info: Setting pplx sub-encoding to TS1/3 on input line 323. -Package textcomp Info: Setting pplj sub-encoding to TS1/3 on input line 324. -Package textcomp Info: Setting ptmx sub-encoding to TS1/4 on input line 325. -Package textcomp Info: Setting ptmj sub-encoding to TS1/4 on input line 326. -)) -(/usr/share/texmf-texlive/tex/latex/base/fontenc.sty -Package: fontenc 2005/09/27 v1.99g Standard LaTeX package - -(/usr/share/texmf-texlive/tex/latex/base/t1enc.def -File: t1enc.def 2005/09/27 v1.99g Standard LaTeX file -LaTeX Font Info: Redeclaring font encoding T1 on input line 43. -)) -(/usr/share/texmf-texlive/tex/latex/ae/ae.sty -Package: ae 2001/02/12 1.3 Almost European Computer Modern - -(/usr/share/texmf-texlive/tex/latex/base/fontenc.sty -Package: fontenc 2005/09/27 v1.99g Standard LaTeX package - -(/usr/share/texmf-texlive/tex/latex/base/t1enc.def -File: t1enc.def 2005/09/27 v1.99g Standard LaTeX file -LaTeX Font Info: Redeclaring font encoding T1 on input line 43. -) -LaTeX Font Info: Try loading font information for T1+aer on input line 100. - -(/usr/share/texmf-texlive/tex/latex/ae/t1aer.fd -File: t1aer.fd 1997/11/16 Font definitions for T1/aer. -)))) (./chemCal.aux) -\openout1 = `chemCal.aux'. - -LaTeX Font Info: Checking defaults for OML/cmm/m/it on input line 9. -LaTeX Font Info: ... okay on input line 9. -LaTeX Font Info: Checking defaults for T1/cmr/m/n on input line 9. -LaTeX Font Info: ... okay on input line 9. -LaTeX Font Info: Checking defaults for OT1/cmr/m/n on input line 9. -LaTeX Font Info: ... okay on input line 9. -LaTeX Font Info: Checking defaults for OMS/cmsy/m/n on input line 9. -LaTeX Font Info: ... okay on input line 9. -LaTeX Font Info: Checking defaults for OMX/cmex/m/n on input line 9. -LaTeX Font Info: ... okay on input line 9. -LaTeX Font Info: Checking defaults for U/cmr/m/n on input line 9. -LaTeX Font Info: ... okay on input line 9. -LaTeX Font Info: Checking defaults for PD1/pdf/m/n on input line 9. -LaTeX Font Info: ... okay on input line 9. -LaTeX Font Info: Checking defaults for TS1/cmr/m/n on input line 9. -LaTeX Font Info: Try loading font information for TS1+cmr on input line 9. - -(/usr/share/texmf-texlive/tex/latex/base/ts1cmr.fd -File: ts1cmr.fd 1999/05/25 v2.5h Standard LaTeX font definitions -) -LaTeX Font Info: ... okay on input line 9. -Package hyperref Info: Link coloring OFF on input line 9. - -(/usr/share/texmf-texlive/tex/latex/hyperref/nameref.sty -Package: nameref 2007/05/29 v2.31 Cross-referencing by name of section - -(/usr/share/texmf-texlive/tex/latex/oberdiek/refcount.sty -Package: refcount 2008/08/11 v3.1 Data extraction from references (HO) -) -\c@section@level=\count98 -) -LaTeX Info: Redefining \ref on input line 9. -LaTeX Info: Redefining \pageref on input line 9. - (./chemCal.out) -(./chemCal.out) -\@outlinefile=\write4 -\openout4 = `chemCal.out'. - -\AtBeginShipoutBox=\box27 - (/usr/share/texmf-texlive/tex/context/base/supp-pdf.mkii -[Loading MPS to PDF converter (version 2006.09.02).] -\scratchcounter=\count99 -\scratchdimen=\dimen108 -\scratchbox=\box28 -\nofMPsegments=\count100 -\nofMParguments=\count101 -\everyMPshowfont=\toks16 -\MPscratchCnt=\count102 -\MPscratchDim=\dimen109 -\MPnumerator=\count103 -\everyMPtoPDFconversion=\toks17 -) -LaTeX Font Info: External font `cmex10' loaded for size -(Font) <12> on input line 11. -LaTeX Font Info: External font `cmex10' loaded for size -(Font) <8> on input line 11. -LaTeX Font Info: External font `cmex10' loaded for size -(Font) <6> on input line 11. -LaTeX Font Info: Try loading font information for T1+aett on input line 12. - (/usr/share/texmf-texlive/tex/latex/ae/t1aett.fd -File: t1aett.fd 1997/11/16 Font definitions for T1/aett. -) -<chemCal-001.pdf, id=5, 433.62pt x 433.62pt> -File: chemCal-001.pdf Graphic file (type pdf) - <use chemCal-001.pdf> [1 - -{/var/lib/texmf/fonts/map/pdftex/updmap/pdftex.map}] -<chemCal-002.pdf, id=17, 433.62pt x 433.62pt> -File: chemCal-002.pdf Graphic file (type pdf) - <use chemCal-002.pdf> [2 <./chemCal-001.pdf>] -LaTeX Font Info: Try loading font information for TS1+aett on input line 86. - - (/usr/share/texmf/tex/latex/R/tex/latex/ts1aett.fd -File: ts1aett.fd -) -LaTeX Font Info: Try loading font information for TS1+cmtt on input line 86. - - -(/usr/share/texmf-texlive/tex/latex/base/ts1cmtt.fd -File: ts1cmtt.fd 1999/05/25 v2.5h Standard LaTeX font definitions -) -LaTeX Font Info: Font shape `TS1/aett/m/n' in size <10> not available -(Font) Font shape `TS1/cmtt/m/n' tried instead on input line 86. - [3 <./chemCal-002.pdf>] - -LaTeX Font Warning: Font shape `T1/aett/bx/n' undefined -(Font) using `T1/aett/m/n' instead on input line 117. - -LaTeX Font Info: External font `cmex10' loaded for size -(Font) <7> on input line 119. -LaTeX Font Info: External font `cmex10' loaded for size -(Font) <5> on input line 119. -[4] [5] [6] (./chemCal.aux) - -LaTeX Font Warning: Some font shapes were not available, defaults substituted. - - ) -Here is how much of TeX's memory you used: - 5329 strings out of 495021 - 77642 string characters out of 1181036 - 141123 words of memory out of 3000000 - 8405 multiletter control sequences out of 15000+50000 - 21574 words of font info for 52 fonts, out of 3000000 for 9000 - 28 hyphenation exceptions out of 8191 - 38i,7n,21p,255b,444s stack positions out of 5000i,500n,10000p,200000b,50000s - </home/ranke/.texmf-var/fonts/pk/ljfour/jknappen/ec/tctt1000.600pk></usr/sha -re/texmf-texlive/fonts/type1/public/amsfonts/cm/cmbx12.pfb></usr/share/texmf-te -xlive/fonts/type1/public/amsfonts/cm/cmex10.pfb></usr/share/texmf-texlive/fonts -/type1/public/amsfonts/cm/cmmi10.pfb></usr/share/texmf-texlive/fonts/type1/publ -ic/amsfonts/cm/cmmi5.pfb></usr/share/texmf-texlive/fonts/type1/public/amsfonts/ -cm/cmmi7.pfb></usr/share/texmf-texlive/fonts/type1/public/amsfonts/cm/cmr10.pfb -></usr/share/texmf-texlive/fonts/type1/public/amsfonts/cm/cmr12.pfb></usr/share -/texmf-texlive/fonts/type1/public/amsfonts/cm/cmr17.pfb></usr/share/texmf-texli -ve/fonts/type1/public/amsfonts/cm/cmr5.pfb></usr/share/texmf-texlive/fonts/type -1/public/amsfonts/cm/cmr7.pfb></usr/share/texmf-texlive/fonts/type1/public/amsf -onts/cm/cmsltt10.pfb></usr/share/texmf-texlive/fonts/type1/public/amsfonts/cm/c -msy10.pfb></usr/share/texmf-texlive/fonts/type1/public/amsfonts/cm/cmtt10.pfb>< -/usr/share/texmf-texlive/fonts/type1/public/amsfonts/cm/cmtt12.pfb> -Output written on chemCal.pdf (6 pages, 185068 bytes). -PDF statistics: - 121 PDF objects out of 1000 (max. 8388607) - 14 named destinations out of 1000 (max. 500000) - 11 words of extra memory for PDF output out of 10000 (max. 10000000) - diff --git a/trunk/chemCal/vignettes/chemCal.out b/trunk/chemCal/vignettes/chemCal.out deleted file mode 100644 index e69de29..0000000 --- a/trunk/chemCal/vignettes/chemCal.out +++ /dev/null diff --git a/trunk/chemCal/vignettes/chemCal.pdf b/trunk/chemCal/vignettes/chemCal.pdf Binary files differdeleted file mode 100644 index 5a46025..0000000 --- a/trunk/chemCal/vignettes/chemCal.pdf +++ /dev/null |