卡方检验

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)); #卡方检验结果

obs= [21 20 20 18 32 56 38]
exp= [18.63636364 18.63636364 18.63636364 18.63636364 37.27272727 55.90909091 37.27272727]
Power_divergenceResult(statistic=1.2813008130081303, pvalue=0.9726787782729751)