<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Bom dia Mauro,<br>
<br>
Na verdade troquei tempo por epoca e o m1 depende da instalação
do pacote pscl, ficando o CRM corrigido:<br>
<br>
<font face="Arial">require(pscl)<br>
<br>
y1<- c(mapply(rpois, lambda=c(5,20,45),
MoreArgs=list(n=400)))##Criação <br>
da variável resposta Poisson<br>
y2<- c(mapply(rbinom, size=c(1,0,0), prob=c(0.5,1,1), <br>
MoreArgs=list(n=200)))##Criação da variável resposta Binomial<br>
y<-c(y1,y2)<br>
<br>
trat <- as.factor(gl(3,600)) ##Criação dos tratamentos<br>
<br>
tempo<- as.factor(rep(gl(6,100),3)) ### Criação da variável
tempo<br>
<br>
dados<-as.data.frame(cbind(trat,tempo, y))<br>
<br>
#-------------------------------------------------------------------------------<br>
# Análise de variância do dados inflacionados<br>
<br>
summary(m1 <- zeroinfl(y ~ trat | tempo, data = dados)) ##
Modelo completo<br>
<br>
mnull <- update(m1, . ~ 1) ### Modelo nulo<br>
<br>
pchisq(2 * (logLik(m1) - logLik(mnull)), df = 2, lower.tail =
FALSE) ## <br>
Teste de Chi o modelo completo foi significativo<br>
<br>
## Comparando com GLM Poisson sem inflação por zeros <br>
---------------------------<br>
<br>
summary(p1 <- glm(y ~ trat + epoca, family = poisson, data =
dados))<br>
<br>
##Teste de <br>
Vuong----------------------------------------------------------------<br>
<br>
vuong(p1, m1) ## O GLM Poisson estava tão mal ajustado que o
Poisson <br>
iflacionado é a única opção<br>
<br>
#</font><br>
<div class="moz-cite-prefix">Em 22/09/2013 19:34, Mauro Sznelwar
escreveu:<br>
</div>
<blockquote cite="mid:AE3EA39A91A54BE7B53C6506FBEDCD9F@ACER4520"
type="cite">
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
<meta name="GENERATOR" content="MSHTML 8.00.6001.23520">
<style></style>
<div><font color="#0000ff" face="Arial">Não estão sendo
encontrados os objetos 'm1' e 'época'</font></div>
<div><font face="Arial"><br>
</font></div>
<font size="2" face="Arial">Bom dia Pessoal,<br>
<br>
Estou quebrando a cabeça para tentar desdobrar
interações em <br>
dados com distribuição de Poisson inflacionada por zeros, na
verdade é <br>
uma situação muito recorrente em entomologia, onde em algumas
situações <br>
temos zero insetos e em outras muitos, considerando-se dados de
<br>
contagens. Fiz uma pequena simulação abaixo para ilustrar o
problema <br>
(Variáveis observadas Trat e Epoca e variável resposta y), se
acho um <br>
modelo significativo que tem uma parte com distribuição de
Poisson e <br>
outra com Binomial, como faço o desdobramento? Uma para cada <br>
distribuição ou tem alguma abordagem diferente? Segue CRM
abaixo:<br>
<br>
## Desdobramento de dados com distribuição de Poisson
inflacionado por <br>
zeros ---<br>
<br>
require(pscl)<br>
<br>
y1<- c(mapply(rpois, lambda=c(5,20,45),
MoreArgs=list(n=400)))##Criação <br>
da variável resposta Poisson<br>
y2<- c(mapply(rbinom, size=c(1,0,0), prob=c(0.5,1,1), <br>
MoreArgs=list(n=200)))##Criação da variável resposta Binomial<br>
y<-c(y1,y2)<br>
<br>
trat <- as.factor(gl(3,600)) ##Criação dos tratamentos<br>
<br>
tempo<- as.factor(rep(gl(6,100),3)) ### Criação da variável
tempo<br>
<br>
dados<-as.data.frame(cbind(trat,tempo, y))<br>
<br>
#-------------------------------------------------------------------------------<br>
# Análise de variância do dados inflacionados<br>
<br>
summary(m1 <- zeroinfl(y ~ trat | epoca, data = dados)) ##
Modelo completo<br>
<br>
mnull <- update(m1, . ~ 1) ### Modelo nulo<br>
<br>
pchisq(2 * (logLik(m1) - logLik(mnull)), df = 2, lower.tail =
FALSE) ## <br>
Teste de Chi o modelo completo foi significativo<br>
<br>
## Comparando com GLM Poisson sem inflação por zeros <br>
---------------------------<br>
<br>
summary(p1 <- glm(y ~ trat + epoca, family = poisson, data =
dados))<br>
<br>
##Teste de <br>
Vuong----------------------------------------------------------------<br>
<br>
vuong(p1, m1) ## O GLM Poisson estava tão mal ajustado que o
Poisson <br>
iflacionado é a única opção<br>
<br>
#<br>
<br>
## E agora, como desdobrar isto? Faço um desdobramento para
parte de <br>
poisson e outra para parte binomial?<br>
<br>
Obrigado,<br>
<br>
-- <br>
======================================================================<br>
Alexandre dos Santos<br>
Proteção Florestal<br>
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de
Mato Grosso<br>
Campus Cáceres<br>
Caixa Postal 244<br>
Avenida dos Ramires, s/n<br>
Bairro: Distrito Industrial<br>
Cáceres - MT CEP: 78.200-000<br>
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO)<br>
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a><br>
</font><a moz-do-not-send="true"
href="mailto:alexandre.santos@cas.ifmt.edu.br"><font size="2"
face="Arial">alexandre.santos@cas.ifmt.edu.br</font></a><br>
<font size="2" face="Arial">Lattes: </font><a
moz-do-not-send="true"
href="http://lattes.cnpq.br/1360403201088680"><font size="2"
face="Arial">http://lattes.cnpq.br/1360403201088680</font></a><br>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">--
======================================================================
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)
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<a class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a class="moz-txt-link-freetext" href="http://lattes.cnpq.br/1360403201088680">http://lattes.cnpq.br/1360403201088680</a>
======================================================================
</pre>
</body>
</html>