diff options
| author | Johannes Ranke <jranke@uni-bremen.de> | 2020-05-11 13:43:40 +0200 |
|---|---|---|
| committer | Johannes Ranke <jranke@uni-bremen.de> | 2020-05-11 13:43:40 +0200 |
| commit | b36ae3d710858ee3ff2907eb2d780e0dff48a4f3 (patch) | |
| tree | b9e075d38233106465481c25b0a777ef043fb1c7 /tests/testthat/setup_script.R | |
| parent | 576fbc9d86f4db3d1be2fbd4e97b3fcd58f43c2b (diff) | |
Analytical solutions for all SFO variants
Diffstat (limited to 'tests/testthat/setup_script.R')
| -rw-r--r-- | tests/testthat/setup_script.R | 22 |
1 files changed, 12 insertions, 10 deletions
diff --git a/tests/testthat/setup_script.R b/tests/testthat/setup_script.R index 86d5089f..58e328cd 100644 --- a/tests/testthat/setup_script.R +++ b/tests/testthat/setup_script.R @@ -39,25 +39,27 @@ fits <- mmkin(models, quiet = TRUE, cores = n_cores) # One metabolite -SFO_SFO <- mkinmod(parent = list(type = "SFO", to = "m1"), - m1 = list(type = "SFO"), +SFO_SFO <- mkinmod(parent = mkinsub("SFO", to = "m1"), + m1 = mkinsub("SFO"), use_of_ff = "min", quiet = TRUE) -SFO_SFO.ff <- mkinmod(parent = list(type = "SFO", to = "m1"), - m1 = list(type = "SFO"), +SFO_SFO.ff <- mkinmod(parent = mkinsub("SFO", to = "m1"), + m1 = mkinsub("SFO"), use_of_ff = "max", quiet = TRUE) SFO_SFO.ff.nosink <- mkinmod( parent = mkinsub("SFO", "m1", sink = FALSE), m1 = mkinsub("SFO"), quiet = TRUE, use_of_ff = "max") +FOMC_SFO <- mkinmod(parent = mkinsub("FOMC", to = "m1"), + m1 = mkinsub("SFO"), quiet = TRUE) -f_sfo_sfo_desolve <- mkinfit(SFO_SFO, - subset(FOCUS_2006_D, value != 0), +# Avoid warning when fitting a dataset where zero value is removed +FOCUS_D <- subset(FOCUS_2006_D, value != 0) + +f_sfo_sfo_desolve <- mkinfit(SFO_SFO, FOCUS_D, solution_type = "deSolve", quiet = TRUE) -f_sfo_sfo_eigen <- mkinfit(SFO_SFO, - subset(FOCUS_2006_D, value != 0), +f_sfo_sfo_eigen <- mkinfit(SFO_SFO, FOCUS_D, solution_type = "eigen", quiet = TRUE) -f_sfo_sfo.ff <- mkinfit(SFO_SFO.ff, - subset(FOCUS_2006_D, value != 0), +f_sfo_sfo.ff <- mkinfit(SFO_SFO.ff, FOCUS_D, quiet = TRUE) SFO_lin_a <- synthetic_data_for_UBA_2014[[1]]$data |