From 280d36230052de4f94e384648c1283031fbc9840 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Tue, 17 Jul 2018 17:29:14 +0200 Subject: Fix inverse predictions for replicate measurements For details, see NEWS.md --- man/din32645.Rd | 13 ++++++------- 1 file changed, 6 insertions(+), 7 deletions(-) (limited to 'man/din32645.Rd') diff --git a/man/din32645.Rd b/man/din32645.Rd index 12c641a..ffcbaed 100644 --- a/man/din32645.Rd +++ b/man/din32645.Rd @@ -10,19 +10,18 @@ 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)) +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) +round(prediction$Confidence, 5) ## Critical value: -(crit <- lod(m, alpha = 0.01, beta = 0.5)) +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") @@ -40,12 +39,12 @@ 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) +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 +3.04 * lod(m, alpha = 0.01, beta = 0.5)$x } \references{ DIN 32645 (equivalent to ISO 11843), Beuth Verlag, Berlin, 1994 -- cgit v1.2.1