From 764531edb7c5598c7b1e401d6e2028ec832db1c4 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Tue, 17 Jul 2018 19:24:17 +0200 Subject: Canonicalize link into bug tracking system Static documentation rebuilt by pkgdown::build_site() --- docs/articles/chemCal.html | 2 +- docs/news/index.html | 5 +++-- vignettes/chemCal.Rmd | 2 +- 3 files changed, 5 insertions(+), 4 deletions(-) diff --git a/docs/articles/chemCal.html b/docs/articles/chemCal.html index 29db7c8..3da10fd 100644 --- a/docs/articles/chemCal.html +++ b/docs/articles/chemCal.html @@ -80,7 +80,7 @@ Basic calibration functions

The 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 Massart et al. (1997) since it also covers the case of weighted regression. Studies of the IUPAC orange book and of DIN 32645 (equivalent to ISO 11843), publications by Currie (1997) and the Analytical Method Committee of the Royal Society of Chemistry (Analytical Methods Committee 1989) and a nice paper by Castells and Castillo (Castells and Castillo 2000) provided some further understanding of the matter.

At the moment, the package consists of four functions (calplot, lod, loq and inverse.predict), working on univariate linear models of class lm or rlm, plus several datasets for validation.

-

A bug report 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.

+

A bug report 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.

In fact, in June 2018 I was made aware of the fact that the inverse prediction method implemented in chemCal version 0.1.37 and before did not take the variance of replicate calibration standards about their means into account, nor the number of replicates when calculating the degrees of freedom. Thanks to PhD student Anna Burniol Figols for reporting this issue!

As a consequence, I rewrote inverse.predict not to automatically work with the mean responses for each calibration standard any more. The example calculations from Massart et al. (1997) can still be reproduced when the regression model is calculated using the means of the calibration data as shown below.

diff --git a/docs/news/index.html b/docs/news/index.html index 62ab315..1266127 100644 --- a/docs/news/index.html +++ b/docs/news/index.html @@ -104,11 +104,12 @@ chemCal 0.2.1 (2018-07-17) Unreleased
diff --git a/vignettes/chemCal.Rmd b/vignettes/chemCal.Rmd index ccbbdc7..de490f2 100644 --- a/vignettes/chemCal.Rmd +++ b/vignettes/chemCal.Rmd @@ -35,7 +35,7 @@ At the moment, the package consists of four functions working on univariate linear models of class `lm` or `rlm`, plus several datasets for validation. -A [bug report](http://bugs.r-project.org/bugzilla3/show_bug.cgi?id=8877) +A [bug report](https://bugs.r-project.org/bugzilla/show_bug.cgi?id=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 -- cgit v1.2.1