[R-br] Quasi AIC

[ASANTOS]

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