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

[ASANTOS]

Muito Obrigado Walmes!!!!
      Problema resolvido,
abraço,
Alexandre


On 09/09/2011 02:24 PM, Walmes Zeviani wrote:
> 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 <mailto:walmes em ufpr.br>
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes 
> <http://www.leg.ufpr.br/%7Ewalmes>
> linux user number: 531218
> ==========================================================================
>
>
> 2011/9/9 ASANTOS <alexandresantosbr em yahoo.com.br 
> <mailto: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 <mailto: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.
>
>
>
> _______________________________________________
> 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/20110909/742eb488/attachment.html>