1 files changed, 38 insertions, 0 deletions
diff --git a/tests/testthat/test_multistart.R b/tests/testthat/test_multistart.R
new file mode 100644
index 00000000..56eb140c
--- /dev/null
+++ b/tests/testthat/test_multistart.R
@@ -0,0 +1,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)
+  parplot_sfo <- function() parplot(saem_sfo_s_multi, ylim = c(0.5, 2))
+  vdiffr::expect_doppelganger("llhist for sfo fit", llhist_sfo)
+  vdiffr::expect_doppelganger("parplot for sfo fit", parplot_sfo)
+
+  llhist_biphasic <- function() llhist(saem_biphasic_m_multi)
+  parplot_biphasic <- function() parplot(saem_biphasic_m_multi,
+    ylim = c(0.5, 2))
+
+  vdiffr::expect_doppelganger("llhist for biphasic saemix fit", llhist_biphasic)
+  vdiffr::expect_doppelganger("parplot for biphasic saemix fit", parplot_biphasic)
+})