我尝试过查看Python的多个统计模块,但似乎找不到任何支持one-way ANOVA
事后测试的模块。
单因素方差分析
from scipy import stats
f_value, p_value = stats.f_oneway(data1, data2, data3, data4, ...)
这是单向方差分析,它返回F值和p值。
P值低于设置
Tukey-kramer HSD测试可以像
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
------------------------------------------------
请参考这个网站如何设置参数。
统计模型的tukeyhsd不返回P值。
所以,如果你想知道P的值,从这些输出值计算或使用r
我认为Pyvttbl库返回一个ANOVA表,包括事后测试(即TukeyHSD)。事实上,Pyvttbl的简洁之处在于,您也可以对重复测量进行方差分析。
请参阅此处的Anova1way文档