2016-08-30 40 views
0

我試圖在anova上運行TukeyHSD並且不斷收到此錯誤消息:使用方法(「TukeyHSD」)中的錯誤: 沒有將「TukeyHSD」應用於對象的適用方法class「c('anova','data.frame')」在R中的TukeyHSD post hoc測試中出現錯誤

任何人都可以告訴我這個消息的含義是什麼? 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) 
+1

您應該閱讀文檔:'aov < - aov(group.1)',但當然您需要RHS上的因子變量。 – Roland

回答

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