accrate.age {rbacon} | R Documentation |
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,
ages = c(),
BCAD = set$BCAD,
silent = TRUE,
na.rm = FALSE
)
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 |
cmyr |
Accumulation rates can be calculated in cm/year or year/cm. By default |
ages |
The ages of the age-depth model. Not provided by default, but can be provided to speed things up if the function is called repeatedly |
BCAD |
The calendar scale of graphs and age output-files is in |
silent |
Warn when ages are outside the core's range. Default |
na.rm |
Remove NA entries. These are NOT removed by default, so that always the same amount of iterations is returned. |
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
.
Value
all MCMC estimates of accumulation rate of the chosen age.
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)
accrate.a5000 = accrate.age(5000)
plot(accrate.a5000, pch='.')
hist(accrate.a5000)
## End(Not run)