Fits a function that returns aphid mortality rate based on divergence of maximum and minimum temperatures from development thresholds
Usage
fit_supergaus(Tupper, Tlower, Piupper, Pilower, Psi, Chi)
Arguments
- Tupper
upper development threshold
- Tlower
lower development threshold
- Piupper
mortality rate beyond upper threshold
- Pilower
mortality rate beyond lower threshold
- Psi
survival breadth parameter
- Chi
survival breadth parameter
Value
function that accepts Tmax and Tmin and returns y based on fitted model
Details
A temperature-dependent mortality rate function to pass on to the Endosymbiont model.
Fits a super-gaussian function:
$$
1 - e^ ( -\Pi_{\text{upper}} \cdot (\chi \cdot (\frac{T_{\text{max}} - T_{\text{upper}}}{\Psi})^2) -\Pi_{\text{lower}} \cdot (\chi \cdot (\frac{T_{\text{lower}} - T_{\text{min}}}{\Psi})^2) )
$$
Examples
my_func <- EndoSim:::fit_supergaus(32.79, 2.84, 0.00083, 0.00017, 7.66, 17.66)
plot(seq(-10, 50, by = 0.1), sapply(seq(-10, 50, by = 0.1), function(x) my_func(x, x)), type = "l")
