python で、独立性の検定、期待度数の計算を行う方法を調べてみました。
結果を記載します。


使用するライブラリ

scipy.stats.chi2_contingency使用します。

計算する

統計検定2級 2017年6月の、問13 のデータを使って計算してみます。

from scipy import stats
import numpy as np
array = np.array([[19, 8], [30, 43]])
result = stats.chi2_contingency(array)

print("カイ二乗の統計量",result[0])
print("p値",result[1])
print("自由度",result[2])
print("期待度数",result[3])

カイ二乗の統計量 5.63856186818
p値 0.0175696130794
自由度 1
期待度数 [[ 13.23  13.77]
 [ 35.77  37.23]]

p値が、0.01757 なので、Aが好きか嫌いかは、性別と関係があるになります。


参考

以下、参考にした記事になります。
* Python: SciPyを使った仮説検定 - け日記
* データサイエンスをPythonで試す(1_検定) - あれもPython,これもPython

以上です。

コメント