From c013e35de82825786a12fa98a202219dd241257e Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Tue, 1 Mar 2022 10:42:56 +0100 Subject: Some more tests for nlmixr.mmkin --- tests/testthat/setup_script.R | 4 ++++ tests/testthat/test_mixed.R | 2 ++ 2 files changed, 6 insertions(+) (limited to 'tests') diff --git a/tests/testthat/setup_script.R b/tests/testthat/setup_script.R index 08d899a1..bfa70005 100644 --- a/tests/testthat/setup_script.R +++ b/tests/testthat/setup_script.R @@ -186,6 +186,8 @@ ds_biphasic <- lapply(ds_biphasic_mean, function(ds) { mmkin_sfo_1 <- mmkin("SFO", ds_sfo, quiet = TRUE, error_model = "tc", cores = n_cores) mmkin_dfop_1 <- mmkin("DFOP", ds_dfop, quiet = TRUE, cores = n_cores) mmkin_biphasic <- mmkin(list("DFOP-SFO" = DFOP_SFO), ds_biphasic, quiet = TRUE, cores = n_cores) +mmkin_biphasic_tc <- mmkin(list("DFOP-SFO" = DFOP_SFO), ds_biphasic, quiet = TRUE, cores = n_cores, + error_model = "tc") mmkin_biphasic_mixed <- mixed(mmkin_biphasic) # nlme @@ -204,6 +206,8 @@ saem_biphasic_s <- saem(mmkin_biphasic, transformations = "saemix", quiet = TRUE # nlmixr saem tmp <- capture_output(nlmixr_saem_biphasic <- nlmixr(mmkin_biphasic, est = "saem", control = nlmixr::saemControl(nBurn = 300, nEm = 100, nmc = 9, print = 0))) +tmp <- capture_output(nlmixr_saem_biphasic_tc <- nlmixr(mmkin_biphasic_tc, est = "saem", + control = nlmixr::saemControl(nBurn = 300, nEm = 100, nmc = 9, print = 0))) # The FOCEI fit takes too long... #tmp <- capture_output(nlmixr_focei_biphasic <- nlmixr(mmkin_biphasic, est = "focei", # control = nlmixr::foceiControl(print = 0))) diff --git a/tests/testthat/test_mixed.R b/tests/testthat/test_mixed.R index 2e241931..d30dcaa8 100644 --- a/tests/testthat/test_mixed.R +++ b/tests/testthat/test_mixed.R @@ -83,6 +83,8 @@ test_that("nlmixr results are reproducible for biphasic fits", { test_summary$date.summary <- "Dummy date for testing" test_summary$time <- c(elapsed = "test time 0") + expect_known_output(print(nlmixr_saem_biphasic, digits = 1), "print_nlmixr_saem_biphasic.txt") + expect_known_output(print(nlmixr_saem_biphasic_tc, digits = 1), "print_nlmixr_saem_biphasic_tc.txt") expect_known_output(print(test_summary, digits = 1), "summary_nlmixr_saem_biphasic.txt") }) -- cgit v1.2.1