# NOT RUN {
SFO_SFO_SFO <- mkinmod(T245 = list(type = "SFO", to = "phenol"),
phenol = list(type = "SFO", to = "anisole"),
anisole = list(type = "SFO"))
# }
# NOT RUN {
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)
summary(fit.2)$bpar
endpoints(fit.1)
plot_sep(fit.2)
# }
Run the code above in your browser using DataLab