diff options
Diffstat (limited to 'tests/testthat/test_synthetic_data_for_UBA_2014.R')
| -rw-r--r-- | tests/testthat/test_synthetic_data_for_UBA_2014.R | 49 |
1 files changed, 22 insertions, 27 deletions
diff --git a/tests/testthat/test_synthetic_data_for_UBA_2014.R b/tests/testthat/test_synthetic_data_for_UBA_2014.R index ea762353..fd3e8bba 100644 --- a/tests/testthat/test_synthetic_data_for_UBA_2014.R +++ b/tests/testthat/test_synthetic_data_for_UBA_2014.R @@ -1,4 +1,4 @@ -# Copyright (C) 2015 Johannes Ranke +# Copyright (C) 2015, 2018 Johannes Ranke # Contact: jranke@uni-bremen.de # This file is part of the R package mkin @@ -19,46 +19,41 @@ context("Results for synthetic data established in expertise for UBA (Ranke 2014)") -m_synth_SFO_lin <- mkinmod(parent = mkinsub("SFO", "M1"), - M1 = mkinsub("SFO", "M2"), - M2 = mkinsub("SFO"), - use_of_ff = "max", quiet = TRUE) +test_that("Results are correct for SFO_lin_a", { + m_synth_SFO_lin <- mkinmod(parent = mkinsub("SFO", "M1"), + M1 = mkinsub("SFO", "M2"), + M2 = mkinsub("SFO"), + use_of_ff = "max", quiet = TRUE) + fit_SFO_lin_a <- mkinfit(m_synth_SFO_lin, + synthetic_data_for_UBA_2014[[1]]$data, + quiet = TRUE) + # Results for SFO_lin_a from p. 48 -m_synth_DFOP_par <- mkinmod(parent = mkinsub("DFOP", c("M1", "M2")), - M1 = mkinsub("SFO"), - M2 = mkinsub("SFO"), - use_of_ff = "max", quiet = TRUE) - -fit_SFO_lin_a <- mkinfit(m_synth_SFO_lin, - synthetic_data_for_UBA_2014[[1]]$data, - quiet = TRUE) -fit_DFOP_par_c <- mkinfit(m_synth_DFOP_par, - synthetic_data_for_UBA_2014[[12]]$data, - quiet = TRUE) - -# Results for SFO_lin_a from p. 48 - -test_that("Fitted parameters are correct for SFO_lin_a", { parms <- round(fit_SFO_lin_a$bparms.optim, c(1, 4, 4, 4, 4, 4)) expect_equivalent(parms, c(102.1, 0.7393, 0.2992, 0.0202, 0.7687, 0.7229)) -}) - -test_that("FOCUS chi2 error levels are correct for SFO_lin_a", { errmin <- round(100 * mkinerrmin(fit_SFO_lin_a)$err.min, 2) expect_equivalent(errmin, c(8.45, 8.66, 10.58, 3.59)) }) # Results for DFOP_par_c from p. 54 -test_that("Fitted parameters are correct for DFOP_par_c", { +test_that("Results are correct for DFOP_par_c", { + skip_on_cran() + m_synth_DFOP_par <- mkinmod(parent = mkinsub("DFOP", c("M1", "M2")), + M1 = mkinsub("SFO"), + M2 = mkinsub("SFO"), + use_of_ff = "max", quiet = TRUE) + + + fit_DFOP_par_c <- mkinfit(m_synth_DFOP_par, + synthetic_data_for_UBA_2014[[12]]$data, + quiet = TRUE) + parms <- round(fit_DFOP_par_c$bparms.optim, c(1, 4, 4, 4, 4, 4, 4, 4)) expect_equal(parms, c(parent_0 = 103.0, k_M1 = 0.0389, k_M2 = 0.0095, f_parent_to_M1 = 0.5565, f_parent_to_M2 = 0.3784, k1 = 0.3263, k2 = 0.0202, g = 0.7130)) -}) - -test_that("FOCUS chi2 error levels are correct for DFOP_par_c", { errmin <- round(100 * mkinerrmin(fit_DFOP_par_c)$err.min, 2) expect_equivalent(errmin, c(4.03, 3.05, 5.07, 3.17)) }) |