[R-br] Problema com abordagem inflacionada de zeros

ASANTOS alexandresantosbr em yahoo.com.br
Terça Abril 29 08:40:46 BRT 2014


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 em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================



Mais detalhes sobre a lista de discussão R-br