rbacon (version 2.3.9.1)

Bacon.Age.d: Output all ages for a single depth.

Description

Output all MCMC-derived age estimates for a given depth.

Usage

Bacon.Age.d(d, set = get("info"), its = set$output, BCAD = set$BCAD,
  remove = FALSE)

Arguments

d

The depth of which Bacon age estimates are to be returned. Has to be a single depth.

set

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

its

The set of MCMC iterations to be used. Defaults to the entire MCMC output, its=set$output.

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.

remove

Whether or not to remove NA values (ages within slumps)

Value

Outputs all MCMC-derived ages for a given depth.

Details

Obtaining an age-depth model is often only a step towards a goal, e.g., plotting a core's fossil series ('proxies') against calendar time. Bacon.Age.d can be used to list all MCMC-derived age estimates for a given (single) depth, for example to calculate mean ages for a 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.chrono.qub.ac.uk/blaauw/manualBacon_2.3.pdf

Examples

# NOT RUN {
  Bacon(run=FALSE, coredir=tempfile())
  agedepth(age.res=50, d.res=50, d.by=10)
  ages.d20 = Bacon.Age.d(20)
  mean(ages.d20)
# }