A function for computing highest density regions (HDRs)
hdr(date, prob = 0.95)A calibrated Bchron date, via e.g. BchronCalibrate
The desired probability interval, in the range(0, 1)
A list where each element is one of the contiguous sets making up the HDR
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.
# Calibrate multiple ages and summarise them
ages <- BchronCalibrate(
ages = 11553, ageSds = 230,
calCurves = "intcal20"
)
# Get samples
hdr(ages$Date1)
#> $`0.4%`
#> [1] 13004 13025
#>
#> $`77.9%`
#> [1] 13059 13874
#>
#> $`16.4%`
#> [1] 13923 14012
#>