Question

我目前正在为我指导的某些学生进行测试。他们设计了一个野外实验，试图了解磷酸盐，硝酸盐和pH值对三个不同位置（食尸鬼）浮游生物丰富度的影响。我设计了这个方差分析（但是现在我意识到其他一些测试可能最合适），并且得到了我粘贴在下面的输出。方差分析似乎可以运行，但是在尝试了Tukey post-hoc之后，出现了一个错误。数据可能不是参数化的，所以我想知道多元回归模型是否最好？还是我只是错误地设计了ANOVA？

谢谢！

这是我的代码：

practiceanova <- aov(WQAbundancebyGuzzler$Abundance[factor(WQAbundancebyGuzzler$Guzzler..)] ~ WQAbundancebyGuzzler$Phospates, na.rm=TRUE)
anova(practiceanova)
TukeyHSD(practiceanova, na.rm=TRUE)

输出：

> practiceanova<-aov(WQAbundancebyGuzzler$Abundance[factor(WQAbundancebyGuzzler$Guzzler..)] ~ WQAbundancebyGuzzler$Phospates, na.rm=TRUE)
Warning message:
In lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) :
 extra argument ‘na.rm’ will be disregarded 
> anova(practiceanova)
Analysis of Variance Table

Response: WQAbundancebyGuzzler$Abundance[factor(WQAbundancebyGuzzler$Guzzler..)]
                               Df Sum Sq Mean Sq F value
WQAbundancebyGuzzler$Phospates  1 3700.8  3700.8  7.6667
Residuals                       8 3861.7   482.7        
                                Pr(>F)  
WQAbundancebyGuzzler$Phospates 0.02434 *
Residuals                               
---
Signif. codes:  
0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
> TukeyHSD(practiceanova, na.rm=TRUE)
Error in TukeyHSD.aov(practiceanova, na.rm = TRUE) : 
  no factors in the fitted model
In addition: Warning message:
In replications(paste("~", xx), data = mf) :
  non-factors ignored: WQAbundancebyGuzzler$Phospates

R Studio：ANOVA Factor和Tukey发行后的问题

0 个答案: