Alexandre,

Deviance é o chi calculado, veja

y1 <- rpois(50, lambda = 4)
trat <- rep(1:10,5)
mod.y1 <- glm(y1~1, family="poisson")
mod.y2 <- glm(y1~trat, family="poisson")
an <- anova(mod.y1, mod.y2, test="Chisq")
an
pchisq(an$Deviance[2], df=an$Df[2], lower=FALSE)

À disposição.
Walmes.


==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes@ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================

2011/9/9 ASANTOS <alexandresantosbr@yahoo.com.br>

Boa tarde pessoal,

    Gostaria de saber como obter o valor de Chi calculado em uma anova() entre dois modelos usando GLM, exemplicando:

 y1<-rpois(50, lambda = 4)
trat<-rep(1:10,5)
mod.y1<-glm(y1~1, family="poisson")
mod.y2<-glm(y1~trat, family="poisson")
 names(anova(mod.y1,mod.y2,test="Chi"))
[1] "Resid. Df"  "Resid. Dev" "Df"         "Deviance"   "P(>|Chi|)"

O problema é que a anova() me dá o valor do  P(>|Chi|), mas não o valor de Chi calculado, onde eu gostaria de saber se é possível obter esse valor?
Obrigado,

Alexandre

_______________________________________________
R-br mailing list
R-br@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.

_______________________________________________
R-br mailing list
R-br@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.