MultiHPD function

Bayesian highest posterior density regions for a series of MCMC chains