用于重复测量方差分析的事后 HSD 使用 statsmodels 库（python）答案

【问题标题】：post hoc HSD for Repeated measures ANOVA be using statsmodels library (python)用于重复测量方差分析的事后 HSD 使用 statsmodels 库（python）
【发布时间】：2022-12-17 13:17:24
【问题描述】：

我尝试对我的重复测量数据进行事后 HSD。我使用 statsmodels.stats.multicomp.MultiComparison，但我不确定它是否解决了重复测量问题。

我的代码：

from statsmodels.stats.anova import AnovaRM
import statsmodels.stats.multicomp as mc

aovrm2way = AnovaRM(df, amp, 'subject', within=['cond', 'type'])
res2way = aovrm2way.fit()

print(res2way)


comp = mc.MultiComparison(df[amp], df['cond'])
post_hoc_res = comp.tukeyhsd()
post_hoc_res.summary()
print(post_hoc_res.summary())

重复测量这样做是否正确？如果没有，是否有另一个库提供重复测量事后测试的答案？

【问题讨论】：

我怀疑您在这里使用 MultiComparison 是否会给您正在寻找的答案。您在代码的任何地方都没有告诉 MultiComparison 使用重复测量模型。也就是说，虽然我不认为您在上面所做的是正确的，但遗憾的是我（还）不知道如何正确地做，或者当前版本的 sm 是否能够做到这一点。

标签： python statistics statsmodels anova

【解决方案1】：

我的建议是您从 pingouin 库中检查函数“pairwise_tests” https://pingouin-stats.org/generated/pingouin.pairwise_tests.html#pingouin.pairwise_tests。

与 statsmodels.stats.multicomp 不同，pg.pairwise_tests 允许您指定“包含主题内因素的列的名称，即重复测量”。

【讨论】：