from scipy import stats
1 sample test
x = stats.norm.rvs(loc=5,scale=10, size=50)
print(x.mean())
print(stats.ttest_1samp(x, 5.0))
print(stats.ttest_1samp(x, 1.0))
2 sample ttest
x =stats.norm.rvs(loc=5,scale=1,size=50)
y = stats.norm.rvs(loc=2, scale=10,size=50)
print(stats.ttest_ind(x, y)) #方差相等
print(stats.ttest_ind(x, y, equal_var=False)) # 方差不等
#chisquare 卡方
A=[16,18,16,14,12,12]
B=[16,16,16,16,16,8]
print(stats.chisquare(A)) # sum: (Oj-Ej)^2/Ej 是否是平均分布
print(stats.chisquare(A, f_exp=B))