我有一个数据集,我需要对多个变量执行 Tukey_test 而不是成对比较:
我的数据看起来像:
Subject State Type Success
A AZ 1 0.89
B TX 2 0.91
C NJ 3 0.78
A OH 4 0.91
我的代码:
from statsmodels.stats.multicomp import pairwise_tukeyhsd
from statsmodels.formula.api import ols
from statsmodels.base.model import Model
p1_tukey_test = pairwise_tukeyhsd(data["Success"], data["Subject"] )
f.write("\n\n" + str(p1_tukey_test._results_table))
我的输出:
===================================================================
group1 group2 meandiff lower upper reject
---------------------------------------------------------------------
A B 13.8883 12.493 15.2835 True
A C 10.2884 8.5795 11.9972 True
我希望输出如下:
===================================================================
group1 group2 group3 meandiff lower upper reject
---------------------------------------------------------------------
A B C 13.8883 12.493 15.2835 True
换句话说,我不想进行成对比较,而是想比较每个主题变量类别的总体结果