I have learned how to find the 95% confidence interval with scipy.stats.t
like so
In [1]: from scipy.stats import t
In [2]: t.interval(0.95, 10, loc
You don't need .interval
method, to get the size of confidence interval, you just need the .ppf
method.
import numpy as np
import scipy.stats as ss
data_m=np.array([1,2,3,4]) #(Means of your data)
data_df=np.array([5,6,7,8]) #(Degree-of-freedoms of your data)
data_sd=np.array([11,12,12,14]) #(Standard Deviations of your data)
import matplotlib.pyplot as plt
plt.errorbar([0,1,2,3], data_m, yerr=ss.t.ppf(0.95, data_df)*data_sd)
plt.xlim((-1,4))
ss.t.ppf(0.95, data_df)*data_sd
is a fully vectorize way to get the (half) size of interval, given the degrees of freedom and standard deviation.