Fit Wang model for development

Fits a function that returns aphid nymph development rate based on daily temperature

Usage

fit_wang(H, r, Topt, Tupper, Tlower, bound)

Arguments

H: development rate scalar
r: development rate scalar
Topt: optimum development temperature
Tupper: upper development threshold
Tlower: lower development threshold
bound: boundary width of thresholds

Value

function that accepts x and returns y based on fitted model

Details

A development rate function to pass on to the Endosymbiont model. Fits a curve using the Wang model (Wang et al. 1982): $$ \left(\frac{H}{{1 + e^{-r \cdot (T - T_{\text{opt}})}}}\right) \cdot \left(1 - e^{- \frac{T - T_{\text{lower}}}{\text{bound}}}\right) \cdot \left(1 - e^{- \frac{T_{\text{upper}} - T}{\text{bound}}}\right) $$

References

Wang, R.S., Lan, Z.X., Ting, Y.C., Studies on mathematical models of the relationships between insect development and temperature. Acta Ecologica Sinica, 2, 47-57 (1982)

Examples

my_func <- EndoSim:::fit_wang(24.1, 0.152, 18.09, 32.79, 2.84, 1.8977)
plot(seq(0, 50, by = 0.1), my_func(seq(0, 50, by = 0.1)), type = "l")