From 1ba099d12ece70f77737505c41eaa7fd1c701140 Mon Sep 17 00:00:00 2001 From: Johannes Ranke 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