compute_ci function

Compute confidence interval