我正在进行方差分析，也希望获得F值。到目前为止，它仍在进行Wald卡方检验而不是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个不同的值）。
非常感谢您对此提供的任何帮助！
欢呼