SFO_SFO_SFO <- mkinmod(T245 = list(type = "SFO", to = "phenol"),
phenol = list(type = "SFO", to = "anisole"),
anisole = list(type = "SFO"))
if (FALSE) {
fit.1 <- mkinfit(SFO_SFO_SFO, subset(mccall81_245T, soil == "Commerce"), quiet = TRUE)
summary(fit.1)$bpar
endpoints(fit.1)
# 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)
}
Run the code above in your browser using DataLab