def mean_confidence_interval(data, confidence=0.95): a = 1.0*np.array(data) n = len(a) m, se = np.mean(a), st.sem(a) h = se * st.t._ppf((1+confidence)/2., n-1) return h