dDoubleExponentialGrowth {nimbleCarbon}

Double Exponential Growth Model

Description

Density and random generation of an exponential growth model distribution.

Usage

dDoubleExponentialGrowth(x, a, b, r1, r2, mu, log)

rDoubleExponentialGrowth(n, a, b, 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).

r1

growth rate before change point mu.

r2

growth rate after change point mu.

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.

Value

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

Author(s)

Enrico Crema

Examples

p = list(r1=0.003,r2=-0.001,mu=5200)
modelPlot(model = dDoubleExponentialGrowth,a=6000,b=4000,params=p,alpha = 1)