BchronMCMC {Bchron}R Documentation

Run Bchron model with MCMC

Description

Runs the compound Poisson-Gamma model for radiocarbon and non-radiocarbon dates

Usage

BchronMCMC(Bchrondata, iterations = 0, burnin = 0, thinby = 1, howmany = 1)

Arguments

Bchrondata A Bchron data object with fields for the the input and output files, the name of the core being used, the details of the calibration curve, and details of the model run required.
iterations Number of MCMC iterations required.
burnin Number of discarded burn-in iterations.
thinby Number of iterations to thin out.
howmany Number of iterations before R reports back and updates the screen.

Details

The Bchron MCMC function fits a compound Poisson-Gamma distribution to the increments between the dated levels. This involves a stochastic linear interpolation where the age gaps are Gamma(2,beta) distributed, and the depth gaps are Exponential(lambda). Radiocarbon dates are updated within the function also by MCMC.

References

Haslett and Parnell (2008), A simple monotone process with application to radiocarbon-dated depth chronologies. Applied Statistics 57 (5), 1–20.

See Also

Bchronloaddata


[Package Bchron version 2.1 Index]