
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

Alexandre, A mensagem de aviso informa "'chat' given is < 1", ou seja, é um caso de subdispersão e não superdispersão. À disposição. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

Obrigado Walmes, Mas o meu problema é que tento simular várias situações possíveis, tanto com ou sem sobredispersão e a função só me retorna [1] NA, para ser sincero a função comigo só funcionou com o exemplo do help(QAIC) do pacote, não consigo obter o valor de quasi AIC mesmo quando mudo o número de penalizações por parâmetro ou altero o parâmetro de dispersão, existe alguma outra forma de calcular o quasi AIC, Obrigado, Alexandre Em 03/07/2012 13:47, Walmes Zeviani escreveu:
Alexandre,
A mensagem de aviso informa "'chat' given is < 1", ou seja, é um caso de subdispersão e não superdispersão.
À disposição.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br <mailto:walmes@ufpr.br> twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes <http://www.leg.ufpr.br/%7Ewalmes> linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.
-- 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

Boa noite: Primeiro, quero dizer que eu sou Argentino e estou no nível inicial das aulas de Português, por isso peço perdão por meu portunhol. Eu uso para fazer os modelos que tem resposta com distribuição binomial negativa a função glmmadmb [install.packages("glmmADMB", repos="http://r-forge.r-project.org")] porque ele me permite trabalhar com modelos mistos e o resultado tem o valor de parâmetro de dispersão binomial negativa (alpha). Depois, pode usar a função qAIC do pacote bblme. Agora eu faço uma pergunta. É o mesmo alpha (parâmetro de disperção binomial negativa) que chat? Até mais. Darío Ezequiel Manzoli Laboratorio de Ecología de Enfermedades Facultad de Ciencias Veterinarias Universidad Nacional del Litoral Kreder 2805 - Esperanza (S3080HOF) Santa Fe 03496-426304(int 111-353) web: www.fcv.unl.edu.ar/archivos/servicios/laboratorios/lecen/lecen1.htm From: ASANTOS Sent: Tuesday, July 03, 2012 3:28 PM To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Quasi AIC Obrigado Walmes, Mas o meu problema é que tento simular várias situações possíveis, tanto com ou sem sobredispersão e a função só me retorna [1] NA, para ser sincero a função comigo só funcionou com o exemplo do help(QAIC) do pacote, não consigo obter o valor de quasi AIC mesmo quando mudo o número de penalizações por parâmetro ou altero o parâmetro de dispersão, existe alguma outra forma de calcular o quasi AIC, Obrigado, Alexandre Em 03/07/2012 13:47, Walmes Zeviani escreveu: Alexandre, A mensagem de aviso informa "'chat' given is < 1", ou seja, é um caso de subdispersão e não superdispersão. À disposição. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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. -- 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 R-br@listas.c3sl.ufpr.br 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.
participantes (3)
-
ASANTOS
-
Darío Manzoli
-
Walmes Zeviani