hpd: Calculate highest posterior density

In Maarten14C/IntCal: Radiocarbon Calibration Curves

Calculate highest posterior density

Description

Calculate highest posterior density ranges of calibrated distribution

Usage

hpd(calib, prob = 0.95, return.raw = FALSE, rounded = 1)

Arguments

calib	The calibrated distribution, as returned from caldist()
prob	Probability range which should be calculated. Default prob=0.95.
return.raw	The raw data to calculate hpds can be returned, e.g. to draw polygons of the calibrated distributions. Defaults to return.raw=FALSE.
rounded	Rounding for reported probabilities. Defaults to 1 decimal.

Value

The highest posterior density ranges, as three columns: from age, to age, and the corresponding percentage(s) of the range(s)

Examples

hpd(caldist(130,20))
plot(tmp <- caldist(2450,50), type='l')
abline(v=hpd(tmp)[,1:2], col=4)

Maarten14C/IntCal documentation built on Oct. 11, 2023, 1:25 a.m.