getCredibleIntervals function

Calculate confidence region from an MCMC or similar sample