credint function

Calculate Bayesian credible intervals