Commit a1c3efc1 authored by 王志伟's avatar 王志伟

数据指标波动假设检验统计

parent 8bcfdc70
...@@ -192,7 +192,7 @@ def chi_cal(data): ...@@ -192,7 +192,7 @@ def chi_cal(data):
data['共计'] = data.apply(lambda x: x.sum(), axis=1) data['共计'] = data.apply(lambda x: x.sum(), axis=1)
# print(data) # print(data)
data.loc['共计'] = data.apply(lambda x: x.sum()) data.loc['共计'] = data.apply(lambda x: x.sum())
# print(data) print(data)
t1=data.iloc[0] t1=data.iloc[0]
t2=data.iloc[1] t2=data.iloc[1]
t11_count=t1[0] t11_count=t1[0]
...@@ -207,8 +207,8 @@ def chi_cal(data): ...@@ -207,8 +207,8 @@ def chi_cal(data):
temp2=data.iloc[2] temp2=data.iloc[2]
t11_theory=temp2[0]*rate1 t11_theory=temp2[0]*rate1
t12_theory=temp2[1]*rate1 t12_theory=temp2[1]*rate1
t21_theory = temp2[0] * rate2 t21_theory = temp2[0]*rate2
t22_theory = temp2[1] * rate2 t22_theory = temp2[1]*rate2
#计算卡方值 #计算卡方值
X=(((t11_count-t11_theory)**2)/t11_theory)+(((t12_count-t12_theory)**2)/t12_theory)+(((t21_count-t21_theory)**2)/t21_theory)+(((t22_count-t22_theory)**2)/t22_theory) X=(((t11_count-t11_theory)**2)/t11_theory)+(((t12_count-t12_theory)**2)/t12_theory)+(((t21_count-t21_theory)**2)/t21_theory)+(((t22_count-t22_theory)**2)/t22_theory)
print("卡方值为:{}".format(X)) print("卡方值为:{}".format(X))
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment