dExponentialLogisticGrowth function

Exponential-Logistic Growth Model

Density and random generation of a exponential-logistic growth model distribution.

dExponentialLogisticGrowth(x, a, b, k, r1, r2, mu, log)

rExponentialLogisticGrowth(n, a, b, k, r1, r2, mu)

Arguments

• x: vector of calendar years (in BP).
• a: lower (earliest) limit of the distribution (in BP).
• b: upper (latest) limit of the distribution (in BP).
• k: initial proportion of the carrying capacity (must be between 0 and 1).
• r1: growth rate of the exponential phase.
• r2: growth rate of logistic phase.
• mu: change point (in BP).
• log: TRUE or 1 to return log probability. FALSE or 0 to return probability.
• n: number of random draws. Currently only n = 1 is supported, but the argument exists for standardization of "r" functions.

Returns

For dExponentialLogisticGrowth: the probability (or likelihood) or log probability of an observed date x (in Cal BP). For rExponentialLogisticGrowth a simulated date in Cal BP.

Examples

p = list(r1=-0.001,r2=0.01,mu=5200,k=0.2)
modelPlot(model = dExponentialLogisticGrowth,a=6000,b=4000,params=p,alpha = 1)

Author(s)

Enrico Crema