Python下数理统计T检验和卡方检验

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))