Я делаю анову и хотел бы также получить значение F. Пока что он все еще делает тест на квадрат Вальда вместо F-теста. Я не понимаю, что я делаю неправильно.
Версия пакета автомобилей 3.0.2.
Мой скрипт:
library(nlme)
library(car)
data<-read.csv("season_flux.csv")
flux = data$CO2_ER #CO2_ER, CO2_NEE, GEP, or CH4
model = lme(flux ~ W_add * P_add * DOY, random = ~1|block/plot,
data = data, na.action = na.omit, method = "REML")
contrasts(data$W_add) <- contr.sum
contrasts(data$P_add) <- contr.sum
Anova <- Anova(model, type = "III", test.statistic="F")
summary(Anova)
print(Anova)
И заголовок моих данных:
DOE DOY plot type CO2_ER CO2_NEE GEP CH4 W_add P_add
1 0.4831486 179 25 D 4.615219 NA -0.1560684 -2.4904458 NW NP
2 0.4885537 179 25 L NA 4.459151 NA NA NW NP
3 0.4948533 179 26 D 4.178014 NA -1.3169272 -0.9756073 W NP
4 0.5003533 179 26 L NA 2.861087 NA NA W NP
5 0.5069753 179 27 D 3.524708 NA -0.9260853 -1.2356268 NW P
6 0.5123025 179 27 L NA 2.598623 NA NA NW P
treatment block PAR soil_moist soil_T_2 soil_T_5
1 C 1 NA 21.80000 4.6 3.0
2 C 1 206.2 21.80000 4.6 3.0
3 W 1 NA 32.13333 4.2 3.2
4 W 1 278.8 32.13333 4.2 3.2
5 P 1 NA 26.93333 4.2 2.6
6 P 1 303.4 26.93333 4.2 2.6
Эксперимент проводится на 9 кампаниях (DOY имеет 9 различных значений).
Большое спасибо за любую помощь с этим!
Приветствия