Сообщение об ошибке при запуске anova

В настоящее время я получаю это сообщение об ошибке из следующего набора данных при запуске функции anova: summary(aov(Site~Chlo,data=alldata)). Набор данных сохраняется в текстовом файле.

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

Есть ли причина этого?


r anova
person user3170629    schedule 11.01.2014    source источник
comment
Я совершенно уверен, что вы хотите Chlo~Site вместо Site~Chlo. Вам следует прочитать документацию по синтаксису формул. Зависимая переменная находится слева, а предикторы - справа.   -  person Roland    schedule 11.01.2014

Ответы (1)


arrow_upward
1
arrow_downward

Я предполагаю, что вы хотите проверить, отличается ли Chlo между Sites. Таким образом, Chlo относится к левой стороне, а Site — к правой в формуле. С вашими данными:

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 такое маленькое, то не так много мощности, и проведение дисперсионного анализа на самом деле не имеет смысла.

person Roland    schedule 11.01.2014
comment
Вы правильно предположили! В этой таблице гораздо больше информации, год идет до 2012. Какой тест вы бы порекомендовали? - person user3170629; 11.01.2014
comment
Поскольку у вас есть данные временных рядов, вам, вероятно, потребуется провести некоторый анализ временных рядов/регрессии. Но это здесь не по теме. Проведите небольшое исследование или проконсультируйтесь со статистиком. Вы также можете задать вопрос на stats.stackexchange.com. - person Roland; 11.01.2014