R中TukeyHSD事后测试中的错误

Question

我正在尝试在方差分析上运行TukeyHSD并不断收到此错误消息：UseMethod（“ TukeyHSD”）中的错误：没有适用于'TukeyHSD'的适用方法应用于类“ c（'anova'，'data.frame”的对象'）”

谁能告诉我这则消息是什么意思？ 我对R非常陌生，因此不胜感激！ 谢谢！

产生此错误的脚本：

a <- data.frame(list(rep(1,100),rnorm(100,50)))
colnames(a) <- c("ID","FeretMinimum")

b <- data.frame(list(rep(2,100),rnorm(100,50)))
colnames(b) <- c("ID","FeretMinimum")

c <- data.frame(list(rep(3,100),rnorm(100,50)))
colnames(c) <- c("ID","FeretMinimum")

d <- data.frame(list(rep(4,100),rnorm(100,50)))
colnames(d) <- c("ID","FeretMinimum")

ColumnFilter = "FeretMinimum"
IDFilter = "ID"

groups <- rbind(a, b ,c ,d)

group.1 <- lm(groups[[ColumnFilter]]~groups[[IDFilter]], data=groups) 

aov<- anova(group.1) 
print(aov)

posthoc <- TukeyHSD(x=aov, as.factor(groups[[ColumnFilter]]), conf.level=0.95)

Answer 1

错误消息说明了一切。 阅读帮助页面。 该函数期望的是aov-fit对象而不是anova对象。

aov.fit1 <- aov(FeretMinimum ~ factor(ID), data=groups)
posthoc <- TukeyHSD(x=aov.fit1, conf.level=0.95)
posthoc

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

Fit: aov(formula = FeretMinimum ~ factor(ID), data = groups)

$`factor(ID)`
            diff        lwr       upr     p adj
2-1 -0.139461021 -0.4988176 0.2198956 0.7487518
3-1 -0.125291199 -0.4846478 0.2340654 0.8050820
4-1 -0.138264399 -0.4976210 0.2210922 0.7537054
3-2  0.014169822 -0.3451868 0.3735264 0.9996226
4-2  0.001196622 -0.3581600 0.3605532 0.9999998
4-3 -0.012973200 -0.3723298 0.3463834 0.9997102