[R-br] F-value no modelo gls do pacote nlme

Fernando Souza nandodesouza em gmail.com
Ter Mar 17 15:15:23 -03 2020


No capítulo 2 do livro Pinheiro e Bates 2000 a partir da página 57 explica toda algebra utilizada nos modelos lme e gls. Você provavelmente encontrará a resposta lá.

https://books.google.com.br/books?id=y54QDUTmvDcC&pg=PA268&lpg=PA268&dq=variance+structura+nlme&source=bl&ots=L_VWkvrQRS&sig=ACfU3U2xReHee03_Lap44C-nNwq9xyyP4A&hl=pt-BR&sa=X&ved=2ahUKEwi__PbOoozoAhWKdN8KHa2RD7wQ6AEwAnoECAoQAQ#v=onepage&q=variance%20structura%20nlme&f=false

On Mar 17 2020, at 2:12 pm, Luiz Leal por (R-br) <r-br em listas.c3sl.ufpr.br> wrote:
> Prezados, boa tarde
>
> 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.
> Luiz
>
> 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
>
>
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200317/dbff5dd8/attachment.html>


Mais detalhes sobre a lista de discussão R-br