Ana Paula,
O seu exemplo não é reproduzível, você apenas enviou a saída do seu modelo mas eu não tenho como reconstruir os objetos para manipulá-los e estudar uma solução.
Segue uma possível solução com os exemplos disponíveis na documentação da glm().
da <- data.frame(counts = c(18, 17, 15, 20, 10, 20, 25, 13, 12),
outcome = gl(3, 1, 9),
treatment = gl(3, 3))
m0 <- glm(counts ~ outcome + treatment,
data = da,
family = poisson())
deviance(m0)
df.residual(m0)
summary(m0)$dispersion
À disposição.
Walmes.