aboutsummaryrefslogtreecommitdiff
path: root/demo
diff options
context:
space:
mode:
authorranke <ranke@5fad18fb-23f0-0310-ab10-e59a3bee62b4>2006-05-09 20:28:46 +0000
committerranke <ranke@5fad18fb-23f0-0310-ab10-e59a3bee62b4>2006-05-09 20:28:46 +0000
commit8d30b2cd951c992e4f9aa3055054091e18b8b4f0 (patch)
tree82c054517c44d2756690de73381791afba6c4e07 /demo
parentd00de8b829e1a426742df43c25e8982e2f9426af (diff)
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
Diffstat (limited to 'demo')
-rw-r--r--demo/00Index2
-rw-r--r--demo/massart97ex3.R12
-rw-r--r--demo/massart97ex8.R12
3 files changed, 26 insertions, 0 deletions
diff --git a/demo/00Index b/demo/00Index
new file mode 100644
index 0000000..8749adf
--- /dev/null
+++ b/demo/00Index
@@ -0,0 +1,2 @@
+massart97ex3 Analysis of example 3 in Massart (1997)
+massart97ex8 Analysis of example 8 in Massart (1997)
diff --git a/demo/massart97ex3.R b/demo/massart97ex3.R
new file mode 100644
index 0000000..7bf9633
--- /dev/null
+++ b/demo/massart97ex3.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)
+data.frame(xi,ybari,si,wi)
+
+weights <- wi[factor(x)]
+m <- lm(y ~ x,w=weights)
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)

Contact - Imprint