HPDint {rmsb} | R Documentation |
Highest Posterior Density Interval
Description
Adapts code from coda::HPDinterval()
to compute a highest posterior density interval from posterior samples for a single parameter. Quoting from the coda
help file, for each parameter the interval is constructed from the empirical cdf of the sample as the shortest interval for which the difference in the ecdf values of the endpoints is the nominal probability. Assuming that the distribution is not severely multimodal, this is the HPD interval.
Usage
HPDint(x, prob = 0.95)
Arguments
x |
a vector of posterior draws |
prob |
desired probability coverage |
Value
a 2-vector with elements Lower
and Upper
Author(s)
Douglas Bates and Frank Harrell
[Package rmsb version 1.1-0 Index]