Bchron (version 4.7.3)

hdr: Calculate highest density regions for Bchron calibrated ages

Description

A function for computing highest density regions (HDRs)

Usage

hdr(date, prob = 0.95)

Arguments

date

A calibrated Bchron date, via e.g. BchronCalibrate

prob

The desired probability interval, in the range(0, 1)

Value

A list where each element is one of the contiguous sets making up the HDR

Details

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.

See Also

BchronCalibrate

Examples

Run this code
# NOT RUN {
# Calibrate multiple ages and summarise them
ages = BchronCalibrate(ages=11553,ageSds=230,
                       calCurves='intcal20')
# Get samples
hdr(ages$Date1)
# }

Run the code above in your browser using DataLab