web-dev-qa-db-ja.com

pythonは、事後検定（Tukey、Scheffeなど）で一元配置分散分析をサポートする統計モジュールはどれですか？

複数の統計モジュールでPythonを調べてみましたが、one-way ANOVA事後テストをサポートするものが見つからないようです。

pythonstatisticsscipyipythonstatsmodels
18
david_adler

一元配置分散分析は次のように使用できます

from scipy import stats
f_value, p_value = stats.f_oneway(data1, data2, data3, data4, ...)

これは一元配置分散分析であり、F値とP値を返します。
P値が設定を下回っている場合は大きな違いがあります。


Tukey-kramerHSD検定は次のように使用できます

from statsmodels.stats.multicomp import pairwise_tukeyhsd
print pairwise_tukeyhsd(Data, Group)

これは複数比較です。出力は次のようになります

Multiple Comparison of Means - Tukey HSD,FWER=0.05
================================================
group1 group2 meandiff   lower    upper   reject
------------------------------------------------
  0      1    -35.2153 -114.8741 44.4434  False 
  0      2     46.697   -40.4993 133.8932 False 
  0      3    -7.5709    -87.49  72.3482  False 
  1      2    81.9123    5.0289  158.7956  True 
  1      3    27.6444   -40.8751  96.164  False 
  2      3    -54.2679 -131.4209 22.8852  False 
------------------------------------------------

引数の設定方法については、thisサイトを参照してください。
statsmodelsのtukeyhsdはP値を返しません。
したがって、P値を知りたい場合は、これらの出力値から計算するか、Rを使用します。

34
ami_GS

ライブラリ Pyvttbl は、事後テスト（つまり、TukeyHSD）を含むANOVAテーブルを返すと思います。実際、Pyvttblの優れている点は、反復測定でもANOVAを実行できることです。

Anova1wayのドキュメントを参照してください ここ

2
Freddy