R中HSD Tukey测试出错

我正在使用单向方差分析，并希望运行一个事后测试。我一直收到一个错误：

UseMethod（“TukeyHSD”）中出错： 没有适用于“function”类对象的“TukeyHSD”方法

我仍然找不到解决办法

我的数据如下所示：

治疗IND T 17 T 17 T110 T15 T110 T110 T112 T18 T1 1 T18 T114 T19 T110 T110 T1 6 T19 T19 T111 T 1 2 T1 6 T15 T19 T111 T19 T 17 T112 T111 T18 T110 T19 T111 T19 T1 4 T19 T111 T111 T19 T112 T113 T111 T19 T110 T 17 T 17 T18 T111 T1 1 T 2 7 T 2 8 T 2 5 T 2 8 T 2 4 T 2 5 T 2 3 T 2 3 T 2 4 T 2 4 T 2 5 T 2 4 T 2 5 t26 T 2 4 T 2 8 T 2 7 T 2 5 t26 t26 T 2 3 T 2 7 T 2 4 T 2 4 T 2 4 t26 T 2 5 t26 t26 T 2 3 T 2 5 T 2 5 T 2 7 T 2 7 T 2 5 T 2 3 t26 t26 T 2 7 T 2 7 T 2 5 T 2 3 T 2 7 t26 T 2 8 T 2 5 T 2 7 T 2 5 t26 t37 T 3 11 t38 T310 t37 T310 T310 t36 t39 t38 t37 T 3 14 t39 t38 T 3 15 T 3 13 t35 t39 t39 T310 T310 T 3 12 T 3 13 T310 t39 T310 t37 t39 t39 T 3 11 t37 T 3 11 t37 T 3 11 t39 T310 t37 t35 t39 T310 T 3 11 T 3 12 T 3 11 t39 t39 t34 t37 t36

---

T3 4TukeyHSD使用函数

aov

产生的aov类对象。函数

oneway.test

返回htest类的对象。这就是你犯错误的原因。如果要运行

TukeyHSD

，则需要使用

aov

使用您的数据：

TukeyHSD(aov(lm(IND ~ Treatment, data = df1)))

  Tukey multiple comparisons of means
    95% family-wise confidence level

Fit: aov(formula = IND ~ Treatment, data = df1)

$`Treatment`
            diff        lwr       upr    p adj
T2-T1 -3.2726878 -4.3935451 -2.151831 0.000000
T3-T1  0.3803734 -0.7404838  1.501231 0.701296
T3-T2  3.6530612  2.5439410  4.762181 0.000000

数据：

---

df1感谢您的快速帮助！
df1 <- structure(list(Treatment = c("T1", "T1", "T1", "T1", "T1", "T1", 
"T1", "T1", "T1", "T1", "T1", "T1", "T1", "T1", "T1", "T1", "T1", 
"T1", "T1", "T1", "T1", "T1", "T1", "T1", "T1", "T1", "T1", "T1", 
"T1", "T1", "T1", "T1", "T1", "T1", "T1", "T1", "T1", "T1", "T1", 
"T1", "T1", "T1", "T1", "T1", "T1", "T1", "T1", "T2", "T2", "T2", 
"T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", 
"T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", 
"T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", 
"T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", "T2", 
"T2", "T2", "T3", "T3", "T3", "T3", "T3", "T3", "T3", "T3", "T3", 
"T3", "T3", "T3", "T3", "T3", "T3", "T3", "T3", "T3", "T3", "T3", 
"T3", "T3", "T3", "T3", "T3", "T3", "T3", "T3", "T3", "T3", "T3", 
"T3", "T3", "T3", "T3", "T3", "T3", "T3", "T3", "T3", "T3", "T3", 
"T3", "T3", "T3", "T3", "T3", "T3", "T3"), 
                      IND = c(7L, 7L, 10L, 
5L, 10L, 10L, 12L, 8L, 1L, 8L, 14L, 9L, 10L, 10L, 6L, 9L, 9L, 
11L, 2L, 6L, 5L, 9L, 11L, 9L, 7L, 12L, 11L, 8L, 10L, 9L, 11L, 
9L, 4L, 9L, 11L, 11L, 9L, 12L, 13L, 11L, 9L, 10L, 7L, 7L, 8L, 
11L, 1L, 7L, 8L, 5L, 8L, 4L, 5L, 3L, 3L, 4L, 4L, 5L, 4L, 5L, 
6L, 4L, 8L, 7L, 5L, 6L, 6L, 3L, 7L, 4L, 4L, 4L, 6L, 5L, 6L, 6L, 
3L, 5L, 5L, 7L, 7L, 5L, 3L, 6L, 6L, 7L, 7L, 5L, 3L, 7L, 6L, 8L, 
5L, 7L, 5L, 6L, 7L, 11L, 8L, 10L, 7L, 10L, 10L, 6L, 9L, 8L, 7L, 
14L, 9L, 8L, 15L, 13L, 5L, 9L, 9L, 10L, 10L, 12L, 13L, 10L, 9L, 
10L, 7L, 9L, 9L, 11L, 7L, 11L, 7L, 11L, 9L, 10L, 7L, 5L, 9L, 
10L, 11L, 12L, 11L, 9L, 9L, 4L, 7L, 6L, 4L)), 
class = "data.frame", 
row.names = c(NA, -145L))