[R-br] Problema com abordagem inflacionada de zeros

[Leonard de Assis]

http://www.statisticalhorizons.com/zero-inflated-models

-----Original Message-----
From: R-br [mailto:r-br-bounces em listas.c3sl.ufpr.br] On Behalf Of ASANTOS
Sent: Tuesday, April 29, 2014 8:41 AM
To: r-br em 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 em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================

_______________________________________________
R-br mailing list
R-br em 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.