Alguém sabe como o valor F é calculado na tabela de Anova do modelo gls do pacote nlme?
Li a documentação, procurei na internet e não achei.
Desde já agradeço.
y=c(2370,1687,2592,2283,2910,3020,1282,1527,871,1025,825,920,562,321,636,317,485,842,173,127,132,150,129,227,193,71,82,62,96,44)
x=as.factor(c("A","A","A","A","A","A","B","B","B","B","B","B","C","C","C","C","C","C","D","D","D","D","D","D","E","E","E","E","E","E"))
library(nlme)
model=gls(y~x, weights = varExp())
anova(model)
Denom. DF: 25
numDF F-value p-value
(Intercept) 1 148.55052 <.0001
x 4 50.83288 <.0001