[R-br] Luz em glm() poisson vs quasipoisson

Marcelo Cardoso mello marcelo_cm32 em yahoo.com.br
Sexta Março 16 15:04:04 BRT 2012


Gostaria de uma luz na análise abaixo. Usando poisson e
quasipoisson tenho significância no modelo (ANOVA). Porem, NO SUMMARY NÃO TENHO
SIGNIFICÂNCIA PARA OS PARÂMETROS DE QUASIPOISSON.  Então é o seguinte, existe ou não diferença
entre os tratamentos na quasipoisson? 
 
dados<-structure(list(Momento
= structure(c(2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L,
2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 
2L, 2L, 2L, 2L, 2L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L,
2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label =
c("Alternativa", 
"Atual"),
class = "factor"), CT = c(13000L, 12000L, 31000L, 39000L, 
1250000L, 13600000L,
6880000L, 47450000L, 890000L, 33000L, 15400000L, 
130000L, 6075L,
850000L, 690000L, 1400000L, 41300000L, 4900000L, 
14900000L, 19000L,
31600000L, 1900000L, 600000L, 2000L, 400000L, 
840000L, 900L, 900L,
4700000L, 910000L, 1000L, 210L, 90L, 4725L, 
200000L, 480000L, 460000L,
107300000L, 8700000L, 3700000L, 4700000L, 
15000000L, 20000L,
180000L, 10000L, 30500000L, 59000L, 260000L, 
1000L, 100L, 550000L,
300L, 2750L, 46000L, 21500L, 5000L, 1050000L, 
9L, 600L, 110L, 310L,
190L, 50L, 18L, 230L, 135L, 100L, 1000L, 
108L, 230L, 100L,
3100L, 350L, 30L, 2400L, 19500L, 10000000L, 
200L, 10L, 10L, 9L,
9L, 9L, 310000L, 9L, 24000L, 10L, 150L, 9L, 
9L, 20L, 80L, 10L, 20L, 9L, 10L, 240L, 40L, 9L, 9L)), .Names =
c("Momento", 
"CT"),
row.names = c(NA, -100L), class = "data.frame")
 
 
m1<-glm(CT~Momento, family=poisson , data=dados)
m2<-glm(CT~Momento, family=quasipoisson,
data=dados)
anova(m1,
test= "Chisq")
anova(m2,
test= "Chisq")
summary(m1)
summary(m2)
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120316/137b059a/attachment.html>


Mais detalhes sobre a lista de discussão R-br