def mean_confidence_interval(data, confidence=0.95): mean = np.mean(data) conf_int = st.t.interval(confidence, len(data) - 1, loc=mean, scale=st.sem(data)) return mean, conf_int[0], conf_int[1]