dDoubleExponentialGrowth: Double Exponential Growth Model

In nimbleCarbon: Bayesian Analyses of Radiocarbon Dates with NIMBLE

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)

nimbleCarbon documentation built on Aug. 14, 2023, 5:08 p.m.