diff options
| author | ranke <ranke@5fad18fb-23f0-0310-ab10-e59a3bee62b4> | 2006-05-23 07:33:22 +0000 |
|---|---|---|
| committer | ranke <ranke@5fad18fb-23f0-0310-ab10-e59a3bee62b4> | 2006-05-23 07:33:22 +0000 |
| commit | f381f9a6a8a47b89ec25cd627833a7248da7932b (patch) | |
| tree | 3155c1f5b2f5810a453aa8cb8a8f44f5920b01e8 /tests/massart97.R | |
| parent | e12be874ff477509b737ad09bf05144a7fbedac2 (diff) | |
Don't do calplot and lod for linear models from weighted
regression any more, since this is not supported (PR#8877).
git-svn-id: http://kriemhild.uft.uni-bremen.de/svn/chemCal@13 5fad18fb-23f0-0310-ab10-e59a3bee62b4
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) |