From f381f9a6a8a47b89ec25cd627833a7248da7932b Mon Sep 17 00:00:00 2001
From: ranke <ranke@5fad18fb-23f0-0310-ab10-e59a3bee62b4>
Date: Tue, 23 May 2006 07:33:22 +0000
Subject: 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
---
 tests/massart97.R | 12 ++++++++++++
 1 file changed, 12 insertions(+)
 create mode 100644 tests/massart97.R

(limited to 'tests/massart97.R')

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)
-- 
cgit v1.2.1