diff options
author | ranke <ranke@5fad18fb-23f0-0310-ab10-e59a3bee62b4> | 2007-10-01 19:44:04 +0000 |
---|---|---|
committer | ranke <ranke@5fad18fb-23f0-0310-ab10-e59a3bee62b4> | 2007-10-01 19:44:04 +0000 |
commit | 14a5af60a36071f6a9b4471fdf183fd91e89e1cd (patch) | |
tree | 8c845109c3b3e7663b903f3a9d06f7094a4438d8 /trunk/tests/massart97.R | |
parent | 3dec3886b58f73427409d3ef9427c8440420cbc0 (diff) |
Moved everything into the trunk directory, in order to enable branching
git-svn-id: http://kriemhild.uft.uni-bremen.de/svn/chemCal@22 5fad18fb-23f0-0310-ab10-e59a3bee62b4
Diffstat (limited to 'trunk/tests/massart97.R')
-rw-r--r-- | trunk/tests/massart97.R | 25 |
1 files changed, 25 insertions, 0 deletions
diff --git a/trunk/tests/massart97.R b/trunk/tests/massart97.R new file mode 100644 index 0000000..00f837f --- /dev/null +++ b/trunk/tests/massart97.R @@ -0,0 +1,25 @@ +require(chemCal) +data(massart97ex1) +m <- lm(y ~ x, data = massart97ex1) +inverse.predict(m, 15) # 6.1 +- 4.9 +inverse.predict(m, 90) # 43.9 +- 4.9 +inverse.predict(m, rep(90,5)) # 43.9 +- 3.2 + +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) +#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) |