aboutsummaryrefslogtreecommitdiff
path: root/branches
diff options
context:
space:
mode:
authorranke <ranke@5fad18fb-23f0-0310-ab10-e59a3bee62b4>2011-09-12 08:42:38 +0000
committerranke <ranke@5fad18fb-23f0-0310-ab10-e59a3bee62b4>2011-09-12 08:42:38 +0000
commitdf80dceaccbac1387604059060ef75fec62d6557 (patch)
treef81e3469079cc2b3be00dc48a2c6162e006ca4ee /branches
parenta34c0a65edde328590dad5b91738d387f62690d8 (diff)
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
Diffstat (limited to 'branches')
-rw-r--r--branches/0.1/chemCal/DESCRIPTION4
-rw-r--r--branches/0.1/chemCal/man/calplot.lm.Rd9
-rw-r--r--branches/0.1/chemCal/tests/din32645.Rout.save9
-rw-r--r--branches/0.1/chemCal/tests/massart97.Rout.save17
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
>

Contact - Imprint