Imported revision 18, because this was last released on CRAN

with version number 0.1-18



git-svn-id: http://kriemhild.uft.uni-bremen.de/svn/chemCal@24 5fad18fb-23f0-0310-ab10-e59a3bee62b4

1 files changed, 25 insertions, 0 deletions

diff --git a/branches/0.1/chemCal/tests/massart97.R b/branches/0.1/chemCal/tests/massart97.R
new file mode 100644
index 0000000..00f837f
--- /dev/null
+++ b/branches/0.1/chemCal/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)