From 280d36230052de4f94e384648c1283031fbc9840 Mon Sep 17 00:00:00 2001
From: Johannes Ranke <jranke@uni-bremen.de>
Date: Tue, 17 Jul 2018 17:29:14 +0200
Subject: Fix inverse predictions for replicate measurements

For details, see NEWS.md
---
 tests/testthat/test_compare_investr.R | 30 ++++++++++++++++++++++++++++++
 1 file changed, 30 insertions(+)
 create mode 100644 tests/testthat/test_compare_investr.R

(limited to 'tests/testthat/test_compare_investr.R')

diff --git a/tests/testthat/test_compare_investr.R b/tests/testthat/test_compare_investr.R
new file mode 100644
index 0000000..6191c89
--- /dev/null
+++ b/tests/testthat/test_compare_investr.R
@@ -0,0 +1,30 @@
+context("Compare with investr::calibrate")
+
+library(chemCal)
+library(investr)
+
+test_that("Unweighted regressions give same results as investr::calibrate using the Wald method", {
+  compare_investr <- function(object, y_sample) {
+    pred_chemCal <- inverse.predict(object, y_sample)
+    pred_investr <- calibrate(object, y_sample, interval = "Wald")
+    expect_equivalent(pred_chemCal[["Prediction"]], 
+                      pred_investr$estimate)
+    expect_equivalent(pred_chemCal[["Standard Error"]], 
+                      pred_investr$se)
+    expect_equivalent(pred_chemCal[["Confidence Limits"]][1], 
+                      pred_investr$lower)  
+    expect_equivalent(pred_chemCal[["Confidence Limits"]][2],
+                      pred_investr$upper)  
+  }
+  m_tol <- lm(peak_area ~ amount, data = rl95_toluene)
+  compare_investr(m_tol, 1000)
+
+  m_din <- lm(y ~ x, din32645)
+  compare_investr(m_din, 5000)
+
+  m_m1 <- lm(y ~ x, massart97ex1)
+  compare_investr(m_m1, 15)
+
+  m_m3 <- lm(y ~ x, massart97ex3)
+  compare_investr(m_m3, 15)
+})
-- 
cgit v1.2.1