Given an ord_res
object, this function returns the linearized model
that best fits the data in the phase space.
ord_res
object can be obtained using the function det_order()
.
Examples
t <- c(0, 4, 8, 12, 16, 20)
conc <- c(1, 0.51, 0.24, 0.12, 0.07, 0.02)
dframe <- data.frame(t, conc)
res <- det_order(dframe)
#> Reaction order estimated: 1
phase_space(res)
#>
#> Call:
#> stats::lm(formula = log_dx_dt ~ log_x, data = ldframe)
#>
#> Coefficients:
#> (Intercept) log_x
#> -1.9148 0.8886
#>