我正在研究混合设计方差分析,并希望运行TukeyHSD进行后Hoc测试 .
我一直收到错误,"Error in UseMethod(" TukeyHSD ") : no applicable method for 'TukeyHSD' applied to an object of class " c('aovlist','listof')“ .
我在一些F&Q网站上看到过关于这个错误的几个问题,但我仍然找不到解决方案 .
我的数据集如下所示:
subject response time group
1 S1 0.99676 time1 task1
2 S2 1.00220 time1 task1
3 S3 1.00420 time1 task1
4 S4 0.99467 time1 task1
5 S5 0.97906 time1 task1
6 S6 0.99162 time1 task1
7 S7 0.99771 time1 task1
8 S8 1.01780 time1 task2
9 S9 0.98682 time1 task2
10 S10 0.99124 time1 task2
11 S11 1.01670 time1 task2
12 S12 0.99769 time1 task2
13 S13 1.02090 time1 task2
14 S14 1.01740 time1 task2
15 S15 0.98851 time1 task3
16 S16 1.00690 time1 task3
17 S17 0.99717 time1 task3
18 S18 0.98945 time1 task3
19 S19 1.00270 time1 task3
20 S20 1.02690 time1 task3
21 S21 1.00050 time1 task3
22 S1 0.96908 time2 task1
23 S2 0.94024 time2 task1
......
然后,anova结果是:
anova = aov(响应〜(组*时间)错误(主题/时间),数据)摘要(anova)
Error: subject
Df Sum Sq Mean Sq F value Pr(>F)
group 2 0.00381 0.001907 0.701 0.509
Residuals 18 0.04896 0.002720
Error: subject:time
Df Sum Sq Mean Sq F value Pr(>F)
time 3 0.08205 0.027351 42.80 2.68e-14 ***
group:time 6 0.00272 0.000454 0.71 0.643
Residuals 54 0.03451 0.000639
如果我运行TukeyHSD:
TukeyHSD(anova)UseMethod(“TukeyHSD”)中的错误:“TukeyHSD”没有适用于“c('aovlist','listof')类对象的适用方法”
我的命令或数据集有什么问题吗?我知道这是一个非常重要的问题......但如果有人能帮助我,我们将不胜感激!谢谢 .
1 回答
您正在处理的是重复测量ANOVA,您需要对此进行正确的事后测试 . 有关其他信息,请参阅以下链接
Post hoc tests with ezANOVA output
Post hoc test after ANOVA with repeated measures using R
[R] Tukey HSD (or other post hoc tests) following repeated measures ANOVA
我认为你最好用这个指定的错误结构构建一个线性混合效果模型,如上面的链接所示 . 这是一个靠近你的人工示例数据集,以及使用nlme-package构建的模型的multcomp-package的临时测试:
这会输出相当长的可能组合列表,但是我们注意到task1,特别是task1.time3,与预期的完全不同: