These functions compute the highest posterior density intervals
(sometimes called minimum length confidence intervals) for a Bayesian
posterior distribution. The hpd function is used when you have
a function representing the inverse cdf (the common case with
conjugate families). The emp.hpd function is used when you
have realizations of the posterior (when you have results from an MCMC
run).