http://www.statisticalhorizons.com/zero-inflated-models -----Original Message----- From: R-br [mailto:r-br-bounces@listas.c3sl.ufpr.br] On Behalf Of ASANTOS Sent: Tuesday, April 29, 2014 8:41 AM To: r-br@listas.c3sl.ufpr.br Subject: [R-br] Problema com abordagem inflacionada de zeros Boa dia Pessoal, Tenho alguns dados que estou com dificuldade de encontrar uma abordagem de análise, tenho uma variável n.insetos (são dados de contagem de insetos, mas tenho várias parcelas com 0 insetos) e a variável resposta n.ataques (são dados de contagem do número de plantas atacadas, também com zeros), o meu problema é que: quando não tenho insetos, não tenho ataque, mas também tenho as vezes insetos, mas não tenho ataque. Tentei pela abordagem de Poisson inflacionado por zeros, sem sucesso (minha variável explicativa tem inúmeros zeros e é uma informação com implicação biológica que não pode ser descartada). Alguém saberia me dizer que abordagem eu poderia tentar para modelar essa relação? Segue CRM com dados parciais e a tentativa de ajuste inflacionada por zeros: #--------------------------------------------------------------------------- ---- # Definições da sessão. rm(list=ls()) require(pscl) require(VGAM) #--------------------------------------------------------------------------- ---- #Dados n.insetos<-c(4,3,1,2,2,2,1,4,1,1,1,0,2,2,3,3,2,6,2,3,3,4,3,0,2,5,6,3,3,3,5,2 ,1, 2,1,1,0,2,2,1,3,6,3,1,2,0,3,1,0,4,2,2,3,2,0,2,2,2,1,2,4,4,3,5,2,1,5,2,3,2,5, 3,3, 3,4,3,2,5,3,4,1,2,2,2,1,2,3,1,4,2,4,4,1,2,2,9,3,1,2,2,2,1,3,3,3,1,1,2,5,2,1, 1,0, 5,0,4,0,3,2,2,3,4,1,1,2,1,1,2,3,3,2,5,2,2,2,3,1,2,6,1,2,3,1,2) n.ataques<-c(0,1,3,0,1,2,1,1,2,2,0,0,1,0,1,1,0,0,2,0,1,0,2,0,0,1,3,0,1,0,0,1 ,1, 1,0,1,0,0,1,1,3,1,1,1,2,0,1,1,0,0,0,1,0,0,0,1,1,0,0,1,2,0,1,2,0,1,1,2,0,1,0, 0,1, 2,1,1,3,1,3,2,0,1,0,0,1,0,2,2,2,2,2,0,1,0,2,3,2,2,3,0,1,2,1,1,0,1,2,1,2,0,1, 1,0, 1,0,0,0,1,0,0,1,0,1,0,0,2,2,3,1,0,2,1,4,1,1,0,0,3,0,2,0,2,1,0) plot(n.insetos, n.ataques) #--------------------------------------------------------------------------- ---- #--------------------------------------------------------------------------- ---- # Modelo completo. compl.mod <- zeroinfl(n.ataques~n.insetos|n.insetos) coef(compl.mod) summary(compl.mod) #--------------------------------------------------------------------------- ---- Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== _______________________________________________ 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.