R Studio: ANOVA Factor & Tukey Post-Hoc Issue - PullRequest
Новая
       38

R Studio: ANOVA Factor & Tukey Post-Hoc Issue

0 голосов
/ 26 февраля 2019

В настоящее время я провожу тест для некоторых студентов, которых я наставник.Они разработали полевой эксперимент, в котором они пытаются увидеть влияние фосфатов, нитратов и рН на содержание планктона в трех разных местах (Guzzlers).Я разработал этот ANOVA (однако теперь я понимаю, что какой-то другой тест может подойти лучше всего) и получил вывод, который я вставил ниже.ANOVA, похоже, побежал, но после попытки получить 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
...