rbacon (version 2.4.2)

accrate.age: Obtain estimated accumulation rates for any age of a core.

Description

Obtain accumulation rates (in years per cm, so actually sedimentation times) as estimated by the MCMC iterations for any age of a core.

Usage

accrate.age(age, set = get("info"), cmyr = FALSE, BCAD = set$BCAD)

Arguments

age

The age for which the accumulation rates need to be returned.

set

Detailed information of the current run, stored within this session's memory as variable info.

cmyr

Accumulation rates can be calculated in cm/year or year/cm. By default cmyr=FALSE and accumulation rates are calculated in year per cm.

BCAD

The calendar scale of graphs and age output-files is in cal BP by default, but can be changed to BC/AD using BCAD=TRUE.

Value

all MCMC estimates of accumulation rate of the chosen age.

Details

Considering accumulation rates is crucial for age-depth modelling, and even more so if they are subsequently used for calculating proxy influx values, or interpreted as proxy for environmental change such as carbon accumulation. See also accrate.age.ghost, accrate.depth and accrate.depth.ghost. Bacon deals explicitly with accumulation rate and its variability through defining prior distributions. This function obtains accumulation rates (in years per cm, so actually sedimentation times) as estimated by the MCMC iterations for any age of a core. Deals with only 1 age at a time. See also accrate.depth.

References

Blaauw, M. and Christen, J.A., Flexible paleoclimate age-depth models using an autoregressive gamma process. Bayesian Anal. 6 (2011), no. 3, 457--474. https://projecteuclid.org/euclid.ba/1339616472

See Also

http://www.qub.ac.uk/chrono/blaauw/manualBacon_2.3.pdf

Examples

# NOT RUN {
  Bacon(run=FALSE, coredir=tempfile())
  agedepth(yr.res=50, d.res=50, d.by=10)
  accrate.a5000 = accrate.age(5000)
  plot(accrate.a5000, pch='.')
  hist(accrate.a5000)
# }