Bom dia pessoal, Estava tentando comparar modelos baseados em distribuição quasi Poisson com Poisson, dentre outros critérios, também realizar a comparação através de Quasi AIC com AIC (Não sei a aceitabilidade dos senhores quanto ao uso do quasi AIC), porém dentre os adeptos, gostaria de saber porque a função não consegue calcular o quasi AIC quando não há sobredispersão, como em:
QAIC(m.2,chat=deviance(m.2)/df.residual(m.2))## Calculo do quasi AIC [1] NA Mensagens de aviso perdidas: In QAIC(m.2, chat = deviance(m.2)/df.residual(m.2)) : 'chat' given is < 1, increased to 1
Segue CRM completo: require(MuMIn) ##Pacote que calcula quasi AIC set.seed(122) y<- c(mapply(rpois, lambda=c(5,20), MoreArgs=list(n=100)))##Criação da variável resposta trat<- gl(2,100) ##Criação dos tratamentos m.1<- glm(y~trat,family="poisson")### Distribuição Poisson summary(m.1) m.2<- glm(y~trat,family="quasipoisson")### Distribuição quasi-Poisson summary(m.2) QAIC(m.2,chat=deviance(m.2)/df.residual(m.2))## Calculo do quasi AIC # Obrigado, -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304