diff options
author | Johannes Ranke <jranke@uni-bremen.de> | 2019-03-04 14:09:45 +0100 |
---|---|---|
committer | Johannes Ranke <jranke@uni-bremen.de> | 2019-03-04 14:09:45 +0100 |
commit | 9d6a975a0e3850a36889a18b5d9929d27b901ab3 (patch) | |
tree | 2130084703b45ba08a7f94e0d446a081a4ba2894 /man | |
parent | d0b124992232dde30fbff83dd4eeae8fc977d3ba (diff) |
Make some example output more reproducible
In order to avoid some unnecessary documentation rebuilds
Diffstat (limited to 'man')
-rw-r--r-- | man/Extract.mmkin.Rd | 5 | ||||
-rw-r--r-- | man/mccall81_245T.Rd | 6 | ||||
-rw-r--r-- | man/transform_odeparms.Rd | 18 |
3 files changed, 18 insertions, 11 deletions
diff --git a/man/Extract.mmkin.Rd b/man/Extract.mmkin.Rd index e27e83a0..973dc28f 100644 --- a/man/Extract.mmkin.Rd +++ b/man/Extract.mmkin.Rd @@ -37,9 +37,4 @@ # This extracts an mkinfit object with lots of components fits[["FOMC", "B"]] ) - - head( - # The same can be achieved by - fits["SFO", "B", drop = TRUE] - ) } diff --git a/man/mccall81_245T.Rd b/man/mccall81_245T.Rd index 0ea8af63..22368d6d 100644 --- a/man/mccall81_245T.Rd +++ b/man/mccall81_245T.Rd @@ -33,13 +33,15 @@ anisole = list(type = "SFO"))
\dontrun{
fit.1 <- mkinfit(SFO_SFO_SFO, subset(mccall81_245T, soil == "Commerce"), quiet = TRUE)
- summary(fit.1, data = FALSE)
+ summary(fit.1)$bpar
+ endpoints(fit.1)
# No convergence, no covariance matrix ...
# 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)
- summary(fit.2, data = FALSE)
+ summary(fit.2)$bpar
+ endpoints(fit.1)
}
}
\keyword{datasets}
diff --git a/man/transform_odeparms.Rd b/man/transform_odeparms.Rd index a1ccdf26..006897d6 100644 --- a/man/transform_odeparms.Rd +++ b/man/transform_odeparms.Rd @@ -66,11 +66,17 @@ SFO_SFO <- mkinmod( m1 = list(type = "SFO")) # Fit the model to the FOCUS example dataset D using defaults fit <- mkinfit(SFO_SFO, FOCUS_2006_D, quiet = TRUE) -summary(fit, data=FALSE) # See transformed and backtransformed parameters +fit.s <- summary(fit) +# Transformed and backtransformed parameters +print(fit.s$par, 3) +print(fit.s$bpar, 3) \dontrun{ +# Compare to the version without transforming rate parameters fit.2 <- mkinfit(SFO_SFO, FOCUS_2006_D, transform_rates = FALSE, quiet = TRUE) -summary(fit.2, data=FALSE) +fit.2.s <- summary(fit.2) +print(fit.2.s$par, 3) +print(fit.2.s$bpar, 3) } initials <- fit$start$value @@ -88,7 +94,9 @@ SFO_SFO.ff <- mkinmod( use_of_ff = "max") fit.ff <- mkinfit(SFO_SFO.ff, FOCUS_2006_D, quiet = TRUE) -summary(fit.ff, data = FALSE) +fit.ff.s <- summary(fit.ff) +print(fit.ff.s$par, 3) +print(fit.ff.s$bpar, 3) initials <- c("f_parent_to_m1" = 0.5) transformed <- transform_odeparms(initials, SFO_SFO.ff) backtransform_odeparms(transformed, SFO_SFO.ff) @@ -101,7 +109,9 @@ SFO_SFO.ff.2 <- mkinmod( fit.ff.2 <- mkinfit(SFO_SFO.ff.2, FOCUS_2006_D, quiet = TRUE) -summary(fit.ff.2, data = FALSE) +fit.ff.2.s <- summary(fit.ff.2) +print(fit.ff.2.s$par, 3) +print(fit.ff.2.s$bpar, 3) } } \keyword{ manip } |