diff options
| author | Johannes Ranke <jranke@uni-bremen.de> | 2022-03-01 10:42:56 +0100 |
|---|---|---|
| committer | Johannes Ranke <jranke@uni-bremen.de> | 2022-03-01 10:42:56 +0100 |
| commit | c013e35de82825786a12fa98a202219dd241257e (patch) | |
| tree | d0bf144328f1a43b5cd5b5b33ed319829c8eab9d /tests | |
| parent | 1d581d24bc85402770a195b44e453a1449597faf (diff) | |
Some more tests for nlmixr.mmkin
Diffstat (limited to 'tests')
| -rw-r--r-- | tests/testthat/setup_script.R | 4 | ||||
| -rw-r--r-- | tests/testthat/test_mixed.R | 2 |
2 files changed, 6 insertions, 0 deletions
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") }) |