Skip to contents

Growth Model According to Gompertz

Source: R/grow_gompertz.R
grow_gompertz.Rd

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).

Value

vector of dependent variable (y)

Details

The equation used here is: $$y = K * exp(log(y0 / K) * exp(-mumax * time))$$

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))