Comparação entre duas populações com distribuição de Poisson
Bom dia pessoal, Tenho duas populações com distribuição de Poisson e gostaria de saber qual a abordagem mais correta para comprar essa duas populações, se empregando análise de variância sobre um GLM ajustado ou teste de Chi quadrado sendo que os dois testes fornecem o mesmo resultado, sendo o CRM: # Dados y<- c(mapply(rpois, lambda=c(5,5), MoreArgs=list(n=15)))##Criação da variável resposta trat <- gl(2,15) ##Criação da variável resposta #Usando GLM m.1 <- glm(y~trat,family="poisson") anova.m.1 <- anova(m.1, test="Chi") anova.m.1 # #Usando teste Chi^2 chisq.test(y,correct=TRUE) # Obrigado, Alexandre
No GLM você especifica que y ~ lambda_i e que lambda_i = X*beta. No caso em que simulou X terá duas colunas para representar os parâmetros envolvidos nos dois grupos. O GLM é uma abordagem paramétrica, você estima lambda e pode comparar submodelos. No chisq.test você compara se as distribuições de frequência observadas são iguais e não são poisson e iguais. Essa é a diferença. ========================================================================== 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 skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Obrigado Walmes, Dúvida esclarecida, Alexandre --- Em ter, 12/2/13, Walmes Zeviani <walmeszeviani@gmail.com> escreveu: De: Walmes Zeviani <walmeszeviani@gmail.com> Assunto: Re: [R-br] Comparação entre duas populações com distribuição de Poisson Para: r-br@listas.c3sl.ufpr.br Data: Terça-feira, 12 de Fevereiro de 2013, 16:35 No GLM você especifica que y ~ lambda_i e que lambda_i = X*beta. No caso em que simulou X terá duas colunas para representar os parâmetros envolvidos nos dois grupos. O GLM é uma abordagem paramétrica, você estima lambda e pode comparar submodelos. No chisq.test você compara se as distribuições de frequência observadas são iguais e não são poisson e iguais. Essa é a diferença. ========================================================================== 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 skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== -----Anexo incorporado----- _______________________________________________ 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)
-
Alexandre Santos -
Walmes Zeviani