From df80dceaccbac1387604059060ef75fec62d6557 Mon Sep 17 00:00:00 2001 From: ranke Date: Mon, 12 Sep 2011 08:42:38 +0000 Subject: New version, adding a comparison of the meaning of the alpha parameter between calplot and lod. Thanks to Jeremie Lasue from the Los Alamos National Lab (USA) for the hint and the suggestion how to address it. git-svn-id: http://kriemhild.uft.uni-bremen.de/svn/chemCal@27 5fad18fb-23f0-0310-ab10-e59a3bee62b4 --- branches/0.1/chemCal/DESCRIPTION | 4 ++-- branches/0.1/chemCal/man/calplot.lm.Rd | 9 ++++++++- branches/0.1/chemCal/tests/din32645.Rout.save | 9 +++++---- branches/0.1/chemCal/tests/massart97.Rout.save | 17 ++++++++++------- 4 files changed, 25 insertions(+), 14 deletions(-) (limited to 'branches') diff --git a/branches/0.1/chemCal/DESCRIPTION b/branches/0.1/chemCal/DESCRIPTION index 92bc658..323827b 100644 --- a/branches/0.1/chemCal/DESCRIPTION +++ b/branches/0.1/chemCal/DESCRIPTION @@ -1,6 +1,6 @@ Package: chemCal -Version: 0.1-26 -Date: 2007-10-01 +Version: 0.1-27 +Date: 2010-09-12 Title: Calibration functions for analytical chemistry Author: Johannes Ranke Maintainer: Johannes Ranke diff --git a/branches/0.1/chemCal/man/calplot.lm.Rd b/branches/0.1/chemCal/man/calplot.lm.Rd index 6f6d584..b60a032 100644 --- a/branches/0.1/chemCal/man/calplot.lm.Rd +++ b/branches/0.1/chemCal/man/calplot.lm.Rd @@ -30,7 +30,9 @@ The label of the y axis. } \item{alpha}{ - The error tolerance level for the confidence and prediction bands. + 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. @@ -48,6 +50,11 @@ 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) diff --git a/branches/0.1/chemCal/tests/din32645.Rout.save b/branches/0.1/chemCal/tests/din32645.Rout.save index c5ed5a7..00ae02e 100644 --- a/branches/0.1/chemCal/tests/din32645.Rout.save +++ b/branches/0.1/chemCal/tests/din32645.Rout.save @@ -1,7 +1,8 @@ -R : Copyright 2006, The R Foundation for Statistical Computing -Version 2.3.1 (2006-06-01) +R version 2.13.1 (2011-07-08) +Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 +Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. @@ -17,7 +18,6 @@ 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) @@ -39,7 +39,8 @@ $x [1] 0.132904 $y -[1] 3764.977 + 1 +3764.977 > loq <- loq(m, alpha = 0.01) > diff --git a/branches/0.1/chemCal/tests/massart97.Rout.save b/branches/0.1/chemCal/tests/massart97.Rout.save index cb113d0..740dd18 100644 --- a/branches/0.1/chemCal/tests/massart97.Rout.save +++ b/branches/0.1/chemCal/tests/massart97.Rout.save @@ -1,7 +1,8 @@ -R : Copyright 2006, The R Foundation for Statistical Computing -Version 2.3.1 (2006-06-01) +R version 2.13.1 (2011-07-08) +Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 +Platform: x86_64-pc-linux-gnu (64-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. @@ -17,7 +18,6 @@ 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 @@ -75,7 +75,7 @@ $Prediction [1] 5.865367 $`Standard Error` -[1] 0.892611 +[1] 0.8926109 $Confidence [1] 2.478285 @@ -103,7 +103,8 @@ $x [1] 5.406637 $y -[1] 13.63822 + 1 +13.63822 > > loq(m0) @@ -111,13 +112,15 @@ $x [1] 13.97767 $y -[1] 30.62355 + 1 +30.62355 > loq(m, w.loq = 1.67) $x [1] 7.346231 $y -[1] 17.90784 + 1 +17.90784 > -- cgit v1.2.1