diff options
Diffstat (limited to 'tests/massart97.R')
-rw-r--r-- | tests/massart97.R | 25 |
1 files changed, 16 insertions, 9 deletions
diff --git a/tests/massart97.R b/tests/massart97.R index 7170ec4..58119d9 100644 --- a/tests/massart97.R +++ b/tests/massart97.R @@ -1,12 +1,19 @@ -library(chemCal) +require(chemCal) data(massart97ex3) attach(massart97ex3) -yx <- split(y,x) -ybar <- sapply(yx,mean) -s <- round(sapply(yx,sd),digits=2) -w <- round(1/(s^2),digits=3) +yx <- split(y, x) +ybar <- sapply(yx, mean) +s <- round(sapply(yx, sd), digits = 2) +w <- round(1 / (s^2), digits = 3) weights <- w[factor(x)] -m <- lm(y ~ x,w=weights) -# The following concords with the book -inverse.predict(m, 15, ws = 1.67) -inverse.predict(m, 90, ws = 0.145) +m <- lm(y ~ x, w = weights) +#calplot(m) + +inverse.predict(m, 15, ws = 1.67) # 5.9 +- 2.5 +inverse.predict(m, 90, ws = 0.145) # 44.1 +- 7.9 + +m0 <- lm(y ~ x) +lod(m0) + +loq(m0) +loq(m, w.loq = 1.67) |