Calculate highest density regions for Bchron calibrated ages
A function for computing highest density regions (HDRs)
hdr(date, prob = 0.95)
date |
A calibrated Bchron date, via e.g. |
prob |
The desired probability interval, in the range(0, 1) |
The output of this function is a list of contiguous ranges which cover the probability interval requested. A highest density region might have multiple such ranges if the calibrated date is multi-modal. These differ from credible intervals, which are always contiguous but will not be a good representation of a multi-modal probability distribution.
A list where each element is one of the contiguous sets making up the HDR
# Calibrate multiple ages and summarise them ages <- BchronCalibrate( ages = 11553, ageSds = 230, calCurves = "intcal20" ) # Get samples hdr(ages$Date1)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.