diff options
Diffstat (limited to 'demo')
| -rw-r--r-- | demo/00Index | 2 | ||||
| -rw-r--r-- | demo/massart97ex3.R | 12 | ||||
| -rw-r--r-- | demo/massart97ex8.R | 12 |
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) |