運行anova時出現錯誤消息

Question

當我運行anova函數時，我正在從以下數據集中收到此錯誤消息：summary(aov(Site~Chlo,data=alldata))。數據集保存在txt文件下。

Error in levels(x)[x] : only 0's may be mixed with negative subscripts 

In addition: Warning messages: 

1: In model.response(mf, "numeric") : 
    using type = "numeric" with a factor response will be ignored 

2: In Ops.factor(y, z$residuals) : - not meaningful for factors 

Year Site Chlo DAC PARD SST 
2003 Seych 2.95 0.24 -39.36 0.40 
2003 Brazil -2.35 -0.14 22.97 4.03 
2003 Indo 0.42 0.04 6.82 0.60 
2004 Seych 0.20 0.02 -2.30 -0.63 
2004 Brazil -0.22 -0.01 -10.28 -1.22 
2004 Indo 0.32 0.03 15.82 -1.72 

這是背後的原因嗎？

Answer 1

我假設你想測試Chlo是否與Site不同。所以Chlo屬於LHS和Site在公式中的RHS。與您的數據：

DF <- read.table(text="Year Site Chlo DAC PARD SST 
2003 Seych 2.95 0.24 -39.36 0.40 
2003 Brazil -2.35 -0.14 22.97 4.03 
2003 Indo 0.42 0.04 6.82 0.60 
2004 Seych 0.20 0.02 -2.30 -0.63 
2004 Brazil -0.22 -0.01 -10.28 -1.22 
2004 Indo 0.32 0.03 15.82 -1.72", header=TRUE) 

summary(aov(Chlo~Site, data=DF)) 
#   Df Sum Sq Mean Sq F value Pr(>F) 
#Site   2 8.247 4.124 2.043 0.275 
#Residuals 3 6.055 2.018 

所以，它不是顯著不同，但如果n是小沒有太多的權力，並進行方差分析並沒有真正意義。