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

4 files changed, 25 insertions, 14 deletions

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 <jranke@uni-bremen.de>
 Maintainer: Johannes Ranke <jranke@uni-bremen.de>
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 
 
 >