Obtenção do Chi calculado em anova()

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

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<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.

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@ufpr.br <mailto:walmes@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@yahoo.com.br <mailto: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 <mailto: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.
participantes (2)
-
ASANTOS
-
Walmes Zeviani