From 8d30b2cd951c992e4f9aa3055054091e18b8b4f0 Mon Sep 17 00:00:00 2001
From: ranke <ranke@5fad18fb-23f0-0310-ab10-e59a3bee62b4>
Date: Tue, 9 May 2006 20:28:46 +0000
Subject: Added the example data, and validated it by writing a demo producing
 the first four columns from Table 8.5 and the weighted regression.

git-svn-id: http://kriemhild.uft.uni-bremen.de/svn/chemCal@5 5fad18fb-23f0-0310-ab10-e59a3bee62b4
---
 demo/massart97ex8.R | 12 ++++++++++++
 1 file changed, 12 insertions(+)
 create mode 100644 demo/massart97ex8.R

(limited to 'demo/massart97ex8.R')

diff --git a/demo/massart97ex8.R b/demo/massart97ex8.R
new file mode 100644
index 0000000..332bd1d
--- /dev/null
+++ b/demo/massart97ex8.R
@@ -0,0 +1,12 @@
+library(chemCal)
+data(massart97ex3)
+attach(massart97ex3)
+xi <- levels(factor(x))
+yx <- split(y,factor(x))
+ybari <- sapply(yx,mean)
+si <- round(sapply(yx,sd),digits=2)
+wi <- round(1/(si^2),digits=3)
+weights <- wi[factor(x)]
+m <- lm(y ~ x,w=weights)
+inverse.predict(m,15)
+inverse.predict(m,90)
-- 
cgit v1.2.1