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.