credible_prep function

Get credible intervals for a, R, b