Classical logistic growth model written as analytical solution of the differential equation.
Usage
grow_logistic(time, parms)
Arguments
- time
vector of time steps (independent variable)
- parms
named parameter vector of the logistic growth model with:
y0 initial value of population measure
mumax intrinsic growth rate (1/time)
K carrying capacity (max. total concentration of cells)
Value
vector of dependent variable (y).
Details
The equation used is:
$$y = (K * y0) / (y0 + (K - y0) * exp(-mumax * time))$$
Examples
time <- seq(0, 30, length=200)
y <- grow_logistic(time, c(y0=1, mumax=0.5, K=10))[,"y"]
plot(time, y, type="l")