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
_______________________________________________
R-br mailing list
[email protected]
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.
