grow_gompertz {growthrates}

Growth Model According to Gompertz

Description

Gompertz growth model written as analytical solution of the differential equation system.

Usage

grow_gompertz(time, parms)

Arguments

time

vector of time steps (independent variable).

parms

named parameter vector of the Gompertz growth model with:

• y0 initial value of abundance,

• mumax maximum growth rate (1/time),

• K maximum abundance (carrying capacity).

Details

The equation used here is:

y = K * exp(log(y0 / K) * exp(-mumax * time))

Value

vector of dependent variable (y)

Note

The naming of parameter "mumax" was done in analogy to the other growth models, but it turned out that it was not consistent with the maximum growth rate of the population. This can be considered as bug. The function will be removed or replaced in future versions of the package. Please use grow_gompertz2 instead.

References

Tsoularis, A. (2001) Analysis of Logistic Growth Models. Res. Lett. Inf. Math. Sci, (2001) 2, 23-46.

See Also

Other growth models: grow_baranyi(), grow_exponential(), grow_gompertz2(), grow_huang(), grow_logistic(), grow_richards(), growthmodel, ode_genlogistic(), ode_twostep()

Examples

time <- seq(0, 30, length=200)
y    <- grow_gompertz(time, c(y0=1, mumax=.2, K=10))[,"y"]
plot(time, y, type="l", ylim=c(0, 20))