accrate.depth: Obtain estimated accumulation rates as for any depth of a...

View source: R/accrate.R

accrate.depthR Documentation

Obtain estimated accumulation rates as for any depth of a core.

Description

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

Usage

accrate.depth(d, set = get("info"), cmyr = FALSE, na.rm = FALSE)

Arguments

d

The depth for which 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.

na.rm

Remove NA entries. These are NOT removed by default, ensuring that always the same amount of iterations is returned.

Details

should take into account hiatuses

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. 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 depth of a core. Deals with only 1 depth at a time. See also accrate.age.

Value

all MCMC estimates of accumulation rate of the chosen depth.

Author(s)

Maarten Blaauw, J. Andres Christen

Examples

## Not run: 
  Bacon(run=FALSE, coredir=tempfile())
  agedepth(yr.res=50, d.res=50, d.by=10)
  d20 <- accrate.depth(20)
  hist(d20)
  d20 <- accrate.depth(20, cmyr=TRUE) # to calculate accumulation rates in cm/yr
  mean(d20)

## End(Not run)

rbacon documentation built on Sept. 30, 2024, 9:40 a.m.