From e13189e6a34d87fabf650c777d0432057bcacbb9 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Wed, 11 Nov 2020 19:47:09 +0100 Subject: Fix example code that was broken by use_of_ff = "max" which is now the default --- man/mccall81_245T.Rd | 16 ++++++++++------ 1 file changed, 10 insertions(+), 6 deletions(-) (limited to 'man/mccall81_245T.Rd') diff --git a/man/mccall81_245T.Rd b/man/mccall81_245T.Rd index fe8dfcc6..44a0c2cf 100644 --- a/man/mccall81_245T.Rd +++ b/man/mccall81_245T.Rd @@ -29,16 +29,20 @@ } \examples{ SFO_SFO_SFO <- mkinmod(T245 = list(type = "SFO", to = "phenol"), - phenol = list(type = "SFO", to = "anisole"), - anisole = list(type = "SFO")) + phenol = list(type = "SFO", to = "anisole"), + anisole = list(type = "SFO")) \dontrun{ fit.1 <- mkinfit(SFO_SFO_SFO, subset(mccall81_245T, soil == "Commerce"), quiet = TRUE) summary(fit.1)$bpar endpoints(fit.1) - # k_phenol_sink is really small, therefore fix it to zero - fit.2 <- mkinfit(SFO_SFO_SFO, subset(mccall81_245T, soil == "Commerce"), - parms.ini = c(k_phenol_sink = 0), - fixed_parms = "k_phenol_sink", quiet = TRUE) + # formation fraction from phenol to anisol is practically 1. As we cannot + # fix formation fractions when using the ilr transformation, we can turn of + # the sink in the model generation + SFO_SFO_SFO_2 <- mkinmod(T245 = list(type = "SFO", to = "phenol"), + phenol = list(type = "SFO", to = "anisole", sink = FALSE), + anisole = list(type = "SFO")) + fit.2 <- mkinfit(SFO_SFO_SFO_2, subset(mccall81_245T, soil == "Commerce"), + quiet = TRUE) summary(fit.2)$bpar endpoints(fit.1) plot_sep(fit.2) -- cgit v1.2.1