[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