1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
|
context("Multistart method for saem.mmkin models")
test_that("multistart works for saem.mmkin models", {
set.seed(123456)
saem_sfo_s_multi <- multistart(sfo_saem_1_reduced, n = 8, cores = n_cores,
no_random_effect = "parent_0")
anova_sfo <- anova(sfo_saem_1,
sfo_saem_1_reduced,
best(saem_sfo_s_multi),
test = TRUE
)
expect_true(anova_sfo[3, "Pr(>Chisq)"] > 0.5)
skip_on_cran() # Save CRAN time
set.seed(123456)
saem_biphasic_m_multi <- multistart(saem_biphasic_m, n = 8,
cores = n_cores)
expect_known_output(print(saem_biphasic_m_multi),
file = "print_multistart_biphasic.txt")
anova_biphasic <- anova(saem_biphasic_m,
best(saem_biphasic_m_multi))
expect_true(anova_biphasic[2, "AIC"] < anova_biphasic[1, "AIC"])
skip_on_travis() # Plots are platform dependent
llhist_sfo <- function() llhist(saem_sfo_s_multi)
parhist_sfo <- function() parhist(saem_sfo_s_multi, ylim = c(0.5, 2))
vdiffr::expect_doppelganger("llhist for sfo fit", llhist_sfo)
vdiffr::expect_doppelganger("parhist for sfo fit", parhist_sfo)
llhist_biphasic <- function() llhist(saem_biphasic_m_multi)
parhist_biphasic <- function() parhist(saem_biphasic_m_multi,
ylim = c(0.5, 2))
vdiffr::expect_doppelganger("llhist for biphasic saemix fit", llhist_biphasic)
vdiffr::expect_doppelganger("parhist for biphasic saemix fit", parhist_biphasic)
})
|
