Indeterminate order rate equation kinetics

Function describing exponential decline from a defined starting value, with a concentration dependent rate constant.

IORE.solution(t, parent_0, k__iore, N)

Arguments

t

Time.

parent_0

Starting value for the response variable at time zero.

k__iore

Rate constant. Note that this depends on the concentration units used.

N

Exponent describing the nonlinearity of the rate equation

Value

The value of the response variable at time t.

Note

The solution of the IORE kinetic model reduces to the SFO.solution if N = 1. The parameters of the IORE model can be transformed to equivalent parameters of the FOMC mode - see the NAFTA guidance for details.

References

NAFTA Technical Working Group on Pesticides (not dated) Guidance for Evaluating and Calculating Degradation Kinetics in Environmental Media

See also

Other parent solutions: DFOP.solution(), FOMC.solution(), HS.solution(), SFO.solution(), SFORB.solution(), logistic.solution()

Examples

  plot(function(x) IORE.solution(x, 100, 0.2, 1.3), 0, 2, ylim = c(0, 100))

  # \dontrun{
    fit.fomc <- mkinfit("FOMC", FOCUS_2006_C, quiet = TRUE)
    fit.iore <- mkinfit("IORE", FOCUS_2006_C, quiet = TRUE)
    fit.iore.deS <- mkinfit("IORE", FOCUS_2006_C, solution_type = "deSolve", quiet = TRUE)

    print(data.frame(fit.fomc$par, fit.iore$par, fit.iore.deS$par, 
                     row.names = paste("model par", 1:4)))
#>             fit.fomc.par fit.iore.par fit.iore.deS.par
#> model par 1  85.87489063    85.874890        85.874891
#> model par 2   0.05192238    -4.826631        -4.826631
#> model par 3   0.65096665     1.949403         1.949403
#> model par 4   1.85744396     1.857444         1.857444
    print(rbind(fomc = endpoints(fit.fomc)$distimes, iore = endpoints(fit.iore)$distimes, 
                iore.deS = endpoints(fit.iore)$distimes))
#>              DT50    DT90 DT50back
#> fomc     1.785233 15.1479 4.559973
#> iore     1.785233 15.1479 4.559973
#> iore.deS 1.785233 15.1479 4.559973
  # }