aboutsummaryrefslogtreecommitdiff
path: root/tests/massart97.R
diff options
context:
space:
mode:
authorJohannes Ranke <jranke@uni-bremen.de>2018-07-05 18:44:22 +0200
committerJohannes Ranke <jranke@uni-bremen.de>2018-07-05 18:44:22 +0200
commit9411139beee167c5339e96db448e5dbed19e06bc (patch)
tree1bba625b3a2288c1ede4eefe1c23d3bcfa0f9fcc /tests/massart97.R
parent801e95f1cc7bfcc8480f6f49b9da2331be71d1a9 (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.R25
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)

Contact - Imprint