rbacon (version 3.3.1)

Bacon.d.Age: Output all depths for a single age.

Description

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

Usage

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

Value

Outputs all MCMC-derived depths for a given age.

Arguments

age

The age estimate for which depths are to be returned. Has to be a single age.

set

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

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.

its

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

na.rm

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

Author

Maarten Blaauw, J. Andres Christen

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.d.Age can be used to list all MCMC-derived depths belonging to a given (single) age, for example to calculate mean depths belonging to a modelled depth. This function was kindly written and provided by Timon Netzel (Bonn University). See also Bacon.Age.d, which calculates the ages for a single depth.

Examples

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