IORE.solution. mkin 0.9.44.9000

Usage

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

Description

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

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.

Value

The value of the response variable at time t.

References

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

Examples

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

  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(coef(fit.fomc), coef(fit.iore), coef(fit.iore.deS),
                   row.names = paste("model par", 1:3)))
            coef.fit.fomc. coef.fit.iore. coef.fit.iore.deS.
model par 1    85.87489063      85.874891          85.874890
model par 2     0.05192238      -4.826631          -4.826631
model par 3     0.65096665       1.949403           1.949403

  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