运行anova时出现错误消息

时间:2014-01-11 14:01:56

标签: r anova

当我运行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

这背后有原因吗?

1 个答案:

答案 0 :(得分: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很小则没有太大的力量,并且进行ANOVA并没有多大意义。