diff options
author | Johannes Ranke <jranke@uni-bremen.de> | 2018-07-05 18:44:22 +0200 |
---|---|---|
committer | Johannes Ranke <jranke@uni-bremen.de> | 2018-07-05 18:44:22 +0200 |
commit | 9411139beee167c5339e96db448e5dbed19e06bc (patch) | |
tree | 1bba625b3a2288c1ede4eefe1c23d3bcfa0f9fcc /tests/massart97.R | |
parent | 801e95f1cc7bfcc8480f6f49b9da2331be71d1a9 (diff) |
Maintenance in preparation of improvements
- Switch vignette to html
- Switch tests to testthat
- NEWS.md instead of ChangeLog
- Remove names of y in lists returned by lod and loq
Diffstat (limited to 'tests/massart97.R')
-rw-r--r-- | tests/massart97.R | 25 |
1 files changed, 0 insertions, 25 deletions
diff --git a/tests/massart97.R b/tests/massart97.R deleted file mode 100644 index 00f837f..0000000 --- a/tests/massart97.R +++ /dev/null @@ -1,25 +0,0 @@ -require(chemCal) -data(massart97ex1) -m <- lm(y ~ x, data = massart97ex1) -inverse.predict(m, 15) # 6.1 +- 4.9 -inverse.predict(m, 90) # 43.9 +- 4.9 -inverse.predict(m, rep(90,5)) # 43.9 +- 3.2 - -data(massart97ex3) -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) -#calplot(m) - -inverse.predict(m, 15, ws = 1.67) # 5.9 +- 2.5 -inverse.predict(m, 90, ws = 0.145) # 44.1 +- 7.9 - -m0 <- lm(y ~ x) -lod(m0) - -loq(m0) -loq(m, w.loq = 1.67) |