calc_CrIs function

Calculate Credible Intervals