From 9a5d2378845e584a2a2c6c09f7ef99a82d09273b Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Fri, 15 May 2015 13:04:36 +0200 Subject: Add synthetic data generated for UBA expertise, tests The new tests are based on the results documented in the expertise I wrote for the UBA (currently unpublished). --- tests/testthat/test_synthetic_data_for_UBA_2014.R | 67 +++++++++++++++++++++++ 1 file changed, 67 insertions(+) create mode 100644 tests/testthat/test_synthetic_data_for_UBA_2014.R (limited to 'tests/testthat/test_synthetic_data_for_UBA_2014.R') diff --git a/tests/testthat/test_synthetic_data_for_UBA_2014.R b/tests/testthat/test_synthetic_data_for_UBA_2014.R new file mode 100644 index 00000000..2bf01075 --- /dev/null +++ b/tests/testthat/test_synthetic_data_for_UBA_2014.R @@ -0,0 +1,67 @@ +# Copyright (C) 2015 Johannes Ranke +# Contact: jranke@uni-bremen.de + +# This file is part of the R package mkin + +# mkin is free software: you can redistribute it and/or modify it under the +# terms of the GNU General Public License as published by the Free Software +# Foundation, either version 3 of the License, or (at your option) any later +# version. + +# This program is distributed in the hope that it will be useful, but WITHOUT +# ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS +# FOR A PARTICULAR PURPOSE. See the GNU General Public License for more +# details. + +# You should have received a copy of the GNU General Public License along with +# this program. If not, see + +context("Results for synthetic data established in expertise for UBA (Ranke 2014)") + + +m_synth_SFO_lin <- mkinmod(parent = list(type = "SFO", to = "M1"), + M1 = list(type = "SFO", to = "M2"), + M2 = list(type = "SFO"), use_of_ff = "max") + + +m_synth_DFOP_par <- mkinmod(parent = list(type = "DFOP", to = c("M1", "M2")), + M1 = list(type = "SFO"), + M2 = list(type = "SFO"), use_of_ff = "max") + +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", { + 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)) +}) + +# References: +# Ranke (2014) Prüfung und Validierung von Modellierungssoftware als Alternative +# zu ModelMaker 4.0, Umweltbundesamt Projektnummer 27452 -- cgit v1.2.1