Given a vector of draws from the posterior (predictive) distribution,
compute a prob
their cumulative probability exceeds prob.
HPDregion(post_pred, prob = 0.95, merge_gaps = FALSE)An ordered vector of values comprising the HPD region
vector of draws from the posterior (predictive) distribution
numeric scalar in (0,1) giving the target probability content of the region
logical; if TRUE, add singleton points that are skipped over