From 1ba099d12ece70f77737505c41eaa7fd1c701140 Mon Sep 17 00:00:00 2001
From: Johannes Ranke <jranke@uni-bremen.de>
Date: Fri, 18 Nov 2022 06:50:00 +0100
Subject: Avoid ill-determined fit, it is platform dependent
The previous test lead to an error on the Fedora test systems (gcc and
clang) on CRAN:
...
Results have changed from known value recorded in 'print_dfop_saemix_1.txt'.
old[13:24] vs new[13:24]
""
"Fitted parameters:"
" estimate lower upper"
- "parent_0 100.09 98.94 101.25"
+ "parent_0 100.17 99.01 101.33"
- "log_k1 -2.68 -2.91 -2.45"
+ "log_k1 -2.66 -2.89 -2.43"
"log_k2 -4.12 -4.24 -4.00"
- "g_qlogis -0.41 -0.63 -0.20"
+ "g_qlogis -0.43 -0.65 -0.21"
- "a.1 0.91 0.67 1.15"
+ "a.1 0.91 0.67 1.16"
"b.1 0.05 0.04 0.06"
- "SD.log_k1 0.36 0.21 0.50"
+ "SD.log_k1 0.35 0.21 0.50"
and 2 more ...
---
tests/testthat/print_dfop_saemix_1.txt | 21 ++++++++++-----------
tests/testthat/setup_script.R | 2 +-
2 files changed, 11 insertions(+), 12 deletions(-)
(limited to 'tests')
diff --git a/tests/testthat/print_dfop_saemix_1.txt b/tests/testthat/print_dfop_saemix_1.txt
index f6fda37c..bdc40065 100644
--- a/tests/testthat/print_dfop_saemix_1.txt
+++ b/tests/testthat/print_dfop_saemix_1.txt
@@ -9,16 +9,15 @@ Data:
Likelihood computed by importance sampling
AIC BIC logLik
- 1409 1415 -695
+ 1409 1415 -696
Fitted parameters:
- estimate lower upper
-parent_0 100.09 98.94 101.25
-log_k1 -2.68 -2.91 -2.45
-log_k2 -4.12 -4.24 -4.00
-g_qlogis -0.41 -0.63 -0.20
-a.1 0.91 0.67 1.15
-b.1 0.05 0.04 0.06
-SD.log_k1 0.36 0.21 0.50
-SD.log_k2 0.22 0.13 0.30
-SD.g_qlogis 0.15 -0.09 0.40
+ estimate lower upper
+parent_0 99.92 98.77 101.06
+log_k1 -2.72 -2.95 -2.50
+log_k2 -4.14 -4.27 -4.01
+g_qlogis -0.35 -0.53 -0.16
+a.1 0.92 0.68 1.16
+b.1 0.05 0.04 0.06
+SD.log_k1 0.37 0.23 0.51
+SD.log_k2 0.23 0.14 0.31
diff --git a/tests/testthat/setup_script.R b/tests/testthat/setup_script.R
index 777c998a..362038c3 100644
--- a/tests/testthat/setup_script.R
+++ b/tests/testthat/setup_script.R
@@ -185,7 +185,7 @@ nlme_biphasic <- suppressWarnings(nlme(mmkin_biphasic))
sfo_saem_1 <- saem(mmkin_sfo_1, quiet = TRUE, transformations = "saemix")
sfo_saem_1_reduced <- update(sfo_saem_1, no_random_effect = "parent_0")
dfop_saemix_1 <- saem(mmkin_dfop_1, quiet = TRUE, transformations = "mkin",
- no_random_effect = "parent_0")
+ no_random_effect = c("parent_0", "g_qlogis"))
saem_biphasic_m <- saem(mmkin_biphasic, transformations = "mkin", quiet = TRUE)
saem_biphasic_s <- saem(mmkin_biphasic, transformations = "saemix", quiet = TRUE)
--
cgit v1.2.1