Fit Weibull curve for fecundity

Fits a function that returns daily fecundity based on cohort adult age

Usage

fit_weibull(K, iota, gamma, eta)

Arguments

K: fecundity rate scalar
iota: shape of fecundity curve
gamma: scale of fecundity curve
eta: locations of fecundity curve

Value

function that accepts x and returns y based on fitted model

Details

Fits a curve using the Weibull distribution to describe how daily fecundity changes with adult age: $$ K \cdot \frac{\iota}{\gamma} \cdot \left( \frac{\text{age} - \eta}{\gamma} \right)^{\iota - 1} \cdot e^{-\left( \frac{\text{age} - \eta}{\gamma} \right)^\iota} $$

References

Weibull, W., A statistical distribution function of wide applicability. Journal of Applied Mechanics, 18, 293-297 (1951)

Examples

my_func <- EndoSim:::fit_weibull(58, 1.885, 5.953, 0)
plot(seq(0, 50, by = 0.1), my_func(seq(0, 50, by = 0.1)), type = "l")