<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Desdobrar as deviances tal como fazemos o desdobramento de somas de quadrados em modelos lineares gaussianos seria interessante mas eu acredito que os gráficos com os modelos ajutados (com bandas) já são suficientes para discutir os resultados. Isso é bem mais fácil de obter. Quando o leitor for ver os resultados ele vai gastar 5 segundos na tabela de desdobramento e 2 minutos vendo o seu gráfico.<br>
<br>À disposição.<br>Walmes.<br></div></div><div class="gmail_extra"><br clear="all"><div><div dir="ltr"><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Walmes Marques Zeviani</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Departamento de Estatística - Universidade Federal do Paraná</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">fone: (+55) 41 3361 3573</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">skype: walmeszeviani<br style="font-family:trebuchet ms,sans-serif"></span><span style="font-family:trebuchet ms,sans-serif">homepage: <a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">linux user number: 531218</span><br style="font-family:trebuchet ms,sans-serif"><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span></div>
</div>
<br><br><div class="gmail_quote">2013/9/23 ASANTOS <span dir="ltr"><<a href="mailto:alexandresantosbr@yahoo.com.br" target="_blank">alexandresantosbr@yahoo.com.br</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div text="#000000" bgcolor="#FFFFFF">
Walmes,<br>
<br>
Na verdade eu gostaria de saber quais tratamentos tiveram
mesmo desempenho no tempo e representá-las através de curvas, mas
não sei como lidar com a decomposição das somas de quadrados quando
tenho duas distribuições (Poisson e Binomial), pois espero dois
resultados: 1) Poisson: a quantidade de insetos representada para
cada um dos tratamentos no tempo e; 2)Binomial: a ocorrência ou não
dos insetos em cada um dos tratamentos no tempo. Teria alguma
abordagem para sugerir?<br>
<br>
Obrigado,<br>
<br>
Alexandre<br>
<br>
CRM:<div class="im"><br>
## Desdobramento de dados com distribuição de Poisson inflacionado
por zeros ---<br>
<br>
require(pscl)<br>
<br>
y1<- c(mapply(rpois, lambda=c(5,20,45),
MoreArgs=list(n=400)))##Criação da variável resposta Poisson<br>
y2<- c(mapply(rbinom, size=c(1,0,0), prob=c(0.5,1,1),
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></div><div class="im">
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)
## Teste de Chi o modelo completo foi significativo<br>
<br>
## Comparando com GLM Poisson sem inflação por zeros
---------------------------<br>
<br></div>
summary(p1 <- glm(y ~ trat + tempo, family = poisson, data =
dados))<div class="im"><br>
<br>
##Teste de
Vuong----------------------------------------------------------------<br>
<br>
vuong(p1, m1) ## O GLM Poisson estava tão mal ajustado que o Poisson
iflacionado é a única opção<br>
<br>
#<br>
<br>
## E agora, como desdobrar isto? Faço um desdobramento para parte de
poisson e outra para parte binomial?<br>
<br>
<br>
<br>
</div><div>Em 23/09/2013 09:26, walmes . escreveu:<br>
</div>
<blockquote type="cite"><div><div class="h5">
<div dir="ltr">
<div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Antes de desdobrar a interação é preciso saber
qual o objetivo da análise. É possível desdobrar de várias
formas, fazendo contrastes específicos, comparações múltiplas.
Eu gosto da abordagem de gráficos com os valores preditos
(médias) e intervalos de confiança (ou bandas). Todavia
existem pessoas da escola antiga que querer ver tabelas com
letras ao lado de médias um em CV baixo, infelizmente.<br>
<br>
À disposição.<br>
Walmes.<br>
</div>
<div class="gmail_extra"><br clear="all">
<div><span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Walmes
Marques Zeviani</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">LEG
(Laboratório de Estatística e Geoinformação, 25.450418 S,
49.231759 W)</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Departamento
de Estatística - Universidade Federal do Paraná</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">fone:
<a href="tel:%28%2B55%29%2041%203361%203573" value="+554133613573" target="_blank">(+55) 41 3361 3573</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">VoIP:
(3361 3600) 1053 1173</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">e-mail: <a href="mailto:walmes@ufpr.br" target="_blank">walmes@ufpr.br</a><br>
skype: walmeszeviani<br style="font-family:trebuchet ms,sans-serif">
</span><span style="font-family:trebuchet ms,sans-serif">twitter:
@walmeszeviani</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">homepage:
<a href="http://www.leg.ufpr.br/%7Ewalmes" target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">linux user
number: 531218</span><br style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">==========================================================================</span><br>
</div>
<br>
</div>
</div>
<br>
<fieldset></fieldset>
<br>
</div></div><div class="im"><pre>_______________________________________________
R-br mailing list
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<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>
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.</pre>
</div></blockquote><div class="im">
<br>
<pre 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: <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)
<a href="mailto:e-mails:alexandresantosbr@yahoo.com.br" target="_blank">e-mails:alexandresantosbr@yahoo.com.br</a>
<a href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a href="http://lattes.cnpq.br/1360403201088680" target="_blank">http://lattes.cnpq.br/1360403201088680</a>
======================================================================
</pre>
</div></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></div>