hpd {rice} | R Documentation |
Calculate highest posterior density
Description
Calculate highest posterior density ranges of calibrated distribution
Usage
hpd(
calib,
prob = 0.95,
return.raw = FALSE,
BCAD = FALSE,
age.round = 0,
prob.round = 1,
every = 0.1
)
Arguments
calib |
The calibrated distribution, as returned from caldist() |
prob |
Probability range which should be calculated. Default |
return.raw |
The raw data to calculate hpds can be returned, e.g. to draw polygons of the calibrated distributions. Defaults to |
BCAD |
Which calendar scale to use. Defaults to cal BP, |
age.round |
Rounding for ages. Defaults to 0 decimals. |
prob.round |
Rounding for reported probabilities. Defaults to 1 decimal. |
every |
Yearly precision (defaults to 0.1, as a compromise between speed and accuracy). |
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)
[Package rice version 1.0.0 Index]