0

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

谢谢！

这是我的代码：

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
4

0 回答 0