Get sample ages from a set of Bchron calibrated dates
A function for extracting sample ages from Bchron calibrated dates
sampleAges(CalDates, n_samp = 10000)
CalDates |
A list created from either |
n_samp |
The desired number of samples |
Sometimes it is useful to have a set of sample calendar ages for your calibrated dates. For example the samples might be required to create a credible/confidence interval, or to create another non-trivial function of calibrated dates, such as differences. By default the BchronCalibrate function provides a grid of ages and an associated density, similar to OxCal. This function extracts that information and uses the sample function to output the desired number of samples
A vector of length n_samp containing sample ages for the specified date
# Calibrate multiple ages and summarise them
ages <- BchronCalibrate(
ages = c(3445, 11553, 7456), ageSds = c(50, 230, 110),
calCurves = c("intcal20", "intcal20", "shcal20")
)
# Get samples
age_samples <- sampleAges(ages)
# Create a credible interval and the median for each date
apply(age_samples, 2, quantile, probs = c(0.05, 0.5, 0.95))Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.