dMetselaar_model {bioinactivation} | R Documentation |
Calculates the first derivative of Metselaar model at a given time for the model parameters provided and the environmental conditions given.
dMetselaar_model(t, x, parms, temp_profile)
t |
numeric vector indicating the time of the experiment. |
x |
list with the value of N at t. |
parms |
parameters for the secondary model. No explicit check of their validity is performed (see section Model Parameters). |
temp_profile |
a function that provides the temperature at a given time. |
The model is developed from the isothermal Metselaar model without
taking into
account in the derivation the time dependence of \delta_T
for
non-isothermal temperature profiles.
This function is compatible with the function
predict_inactivation
.
The value of the first derivative of N at time t
as a list.
\frac{dN}{dt} = -N \cdot p \cdot (1/D)^p \cdot (t/Delta)^{p-1}
D(T) = D_{ref} \cdot 10^{- (T-T_ref)/z}
temp_ref: Reference temperature for the calculation.
D_R: D-value at the reference temperature.
z: z-value.
p: shape factor of the Weibull distribution.
Delta: Scaling parameter
For t=0, dN = 0 unless n=1. Hence, a small shift needs to be introduced to t.