Python |
import numpy as np;
from scipy import stats;
obs=np.array([21,20,20,18,32,56,38]); #观察数
bili=np.array([1,1,1,1,2,3,2]); #期望比例
exp=obs.sum()*bili/bili.sum(); #理论数,注意也可以外界直接给出理论数。
print('obs=',obs);
print('exp=',exp);
print(stats.chisquare(obs,exp)); #卡方检验结果 |