np confidence interval

1import numpy as np
2import scipy.stats
3
4
5def mean_confidence_interval(data, confidence=0.95):
6    a = 1.0 * np.array(data)
7    n = len(a)
8    m, se = np.mean(a), scipy.stats.sem(a)
9    h = se * scipy.stats.t.ppf((1 + confidence) / 2., n-1)
10    return m, m-h, m+h