diff options
Diffstat (limited to 'tests/massart97.R')
| -rw-r--r-- | tests/massart97.R | 12 |
1 files changed, 12 insertions, 0 deletions
diff --git a/tests/massart97.R b/tests/massart97.R new file mode 100644 index 0000000..7170ec4 --- /dev/null +++ b/tests/massart97.R @@ -0,0 +1,12 @@ +library(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) +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) |