def get_ci(theta_star, blockratio=<NUM_LIT:1.0>):
<EOL>b_star = np.sort(theta_star[~np.isnan(theta_star)])<EOL>se = np.std(b_star) * np.sqrt(blockratio)<EOL>ci = [b_star[int(len(b_star) * <NUM_LIT>)], b_star[int(len(b_star) * <NUM_LIT>)]]<EOL>return ci<EOL>
Get the confidence interval.