[R-br] Obtenção do Chi calculado em anova()

[Walmes Zeviani]

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 em ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================


2011/9/9 ASANTOS <alexandresantosbr em 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 em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/**cgi-bin/mailman/listinfo/r-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<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/20110909/25240c3f/attachment.html>