blob: 5243f971cf3e4c7939de4b2376e385b5f010c155 (
plain) (
tree)
|
|
context("Batch fitting and diagnosing hierarchical kinetic models")
test_that("Multiple hierarchical kinetic models can be fitted and diagnosed", {
skip_on_cran()
fits_synth_const <- suppressWarnings(
mmkin(c("SFO", "FOMC"), ds_sfo[1:6], cores = n_cores, quiet = TRUE))
fits_synth_tc <- suppressWarnings(
update(fits_synth_const, error_model = "tc"))
expect_known_output(
convergence(fits_synth_const),
"convergence_fits_synth_const.txt")
hfits <- mhmkin(list(fits_synth_const, fits_synth_tc))
expect_known_output(
convergence(hfits),
"convergence_hfits_synth.txt")
expect_known_output(
illparms(hfits),
"illparms_hfits_synth.txt")
expect_equal(which.min(AIC(hfits)), 3)
hfit_sfo_tc <- update(hfits[["SFO", "tc"]],
covariance.model = diag(c(0, 1)))
expect_equal(illparms(hfit_sfo_tc), character(0))
})
|