<div dir="ltr"><div><div><div>Olhando por cima,<br><br></div>parece-me que você tem uma mistura de distribuições, não?<br><br></div><div>Talvez o pacote flexmix possa te ajudar nesse aspecto.<br><br><a href="http://rss.acs.unt.edu/Rdoc/library/flexmix/html/flexmix.html">http://rss.acs.unt.edu/Rdoc/library/flexmix/html/flexmix.html</a><br>
</div><div><br></div>abcs<br></div>M<br></div><div class="gmail_extra"><br><br><div class="gmail_quote">2013/9/22 Mauro Sznelwar <span dir="ltr"><<a href="mailto:sznelwar@uol.com.br" target="_blank">sznelwar@uol.com.br</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><u></u>
<div>
<div><font></font><font color="#0000ff" face="Arial">Não estão sendo
encontrados os objetos 'm1' e 'época'</font></div>
<div><font face="Arial"><br><font></font></font></div><font 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: <a href="tel:%28%2B55%29%2065%208132-8112" value="+556581328112" target="_blank">(+55) 65 8132-8112</a> (TIM) <a href="tel:%28%2B55%29%2065%209686-6970" value="+556596866970" target="_blank">(+55) 65 9686-6970</a>
(VIVO)<br><a href="mailto:e-mails%3Aalexandresantosbr@yahoo.com.br" target="_blank">e-mails:alexandresantosbr@yahoo.com.br</a><br>
</font><a href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank"><font face="Arial">alexandre.santos@cas.ifmt.edu.br</font></a><br><font face="Arial">Lattes: </font><a href="http://lattes.cnpq.br/1360403201088680" target="_blank"><font face="Arial">http://lattes.cnpq.br/1360403201088680</font></a><br>
</div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><br>-- <br>Manoel Galdino<br>
<a href="https://sites.google.com/site/galdinomcz/" target="_blank">https://sites.google.com/site/galdinomcz/</a><br>
</div>