calc_CrI function

Calculate Credible Interval