diff options
Diffstat (limited to 'vignettes')
| -rw-r--r-- | vignettes/chemCal.R | 36 | ||||
| -rw-r--r-- | vignettes/chemCal.pdf | bin | 0 -> 165600 bytes |
2 files changed, 36 insertions, 0 deletions
diff --git a/vignettes/chemCal.R b/vignettes/chemCal.R new file mode 100644 index 0000000..701db7b --- /dev/null +++ b/vignettes/chemCal.R @@ -0,0 +1,36 @@ +### R code from vignette source '/home/jranke/git/chemCal/vignettes/chemCal.Rnw' + +################################################### +### code chunk number 1: chemCal.Rnw:38-42 +################################################### +library(chemCal) +data(massart97ex3) +m0 <- lm(y ~ x, data = massart97ex3) +calplot(m0) + + +################################################### +### code chunk number 2: chemCal.Rnw:49-50 +################################################### +plot(m0,which=3) + + +################################################### +### code chunk number 3: chemCal.Rnw:56-63 +################################################### +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) + + +################################################### +### code chunk number 4: chemCal.Rnw:69-71 +################################################### +inverse.predict(m, 15, ws=1.67) +inverse.predict(m, 90, ws = 0.145) + + diff --git a/vignettes/chemCal.pdf b/vignettes/chemCal.pdf Binary files differnew file mode 100644 index 0000000..cf1a413 --- /dev/null +++ b/vignettes/chemCal.pdf |