diff options
author | Johannes Ranke <jranke@uni-bremen.de> | 2020-11-11 19:47:09 +0100 |
---|---|---|
committer | Johannes Ranke <jranke@uni-bremen.de> | 2020-11-11 19:47:09 +0100 |
commit | e13189e6a34d87fabf650c777d0432057bcacbb9 (patch) | |
tree | 4c8e054246a05cc0b5af80ecf247863b47b9931f /man/mccall81_245T.Rd | |
parent | f9a5da41b6bcb8a34e27574ede50ddef2255430f (diff) |
Fix example code that was broken by use_of_ff = "max"
which is now the default
Diffstat (limited to 'man/mccall81_245T.Rd')
-rw-r--r-- | man/mccall81_245T.Rd | 16 |
1 files changed, 10 insertions, 6 deletions
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)
|