A function for extracting sample ages from Bchron calibrated dates
sampleAges(CalDates, n_samp = 10000)A list created from either BchronCalibrate.
The desired number of samples
A vector of length n_samp containing sample ages for the specified date
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
# 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))
#> Date1 Date2 Date3
#> 5% 3584 13104 8035
#> 50% 3706 13434 8229
#> 95% 3825 13839 8382