From 8d30b2cd951c992e4f9aa3055054091e18b8b4f0 Mon Sep 17 00:00:00 2001 From: ranke 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/00Index | 2 ++ demo/massart97ex3.R | 12 ++++++++++++ demo/massart97ex8.R | 12 ++++++++++++ 3 files changed, 26 insertions(+) create mode 100644 demo/00Index create mode 100644 demo/massart97ex3.R create mode 100644 demo/massart97ex8.R (limited to 'demo') 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) -- cgit v1.2.1