<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Walmes,<br>
<br>
Atendi a sua sugestão e resolvi ajustar os modelos, no
entanto, fiz as comparações múltiplas para saber a diferença entre
os tratamentos e tudo bem, todos foram significativamente
diferentes, porém na hora de fazer as curvas para cada tratamento no
tempo não consigo obter os coeficientes do modelo para cada um dos
tratamentos ou seja quando faço o summary() do modelo completo tenho
apenas um coeficiente para trat e outro para tempo. Gostrai o que
posso fazer para obter os coeficientes para cada tratamento? Segue
CRM:<br>
<br>
## Desdobramento de dados com distribuição de Poisson inflacionado
por zeros ---<br>
<br>
require(pscl)<br>
require(multicomp)<br>
<br>
y1<- c(mapply(rpois, lambda=c(4,45,93),
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>
summary(m1 <- zeroinfl(y ~ trat*tempo | 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>
summary(p1 <- glm(y ~ trat + tempo, family = poisson, data =
dados))<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>
### Interação entre os tratamamentos no
tempo-----------------------------------<br>
<br>
zi_m1<-zeroinfl(y ~ trat*tempo | trat*tempo, data = dados, link =
"logit")<br>
summary(zinb_edanx2) ## Modelos completo<br>
<br>
## Calculando os contrastes para o modelo inflacionado--------------------------<br>
nr <- length(levels(trat))<br>
contr <- matrix(0, nrow = nr, ncol = length(coef(zi_m1)))<br>
colnames(contr) <- names(coef(zi_m1))<br>
rownames(contr) <- paste(levels(trat)[c(2, 3, 3)],<br>
levels(trat)[c(1, 1, 2)], sep = " - ")<br>
contr[,3:4] <- contrMat(numeric(nrow(contr)), type =
"Tukey")[,-2]<br>
glht_zi <- glht(zi_m1, linfct = contr)<br>
<br>
## Comparações multiplas-------------------------------------------------------------------<br>
summary(glht_zi)<br>
<br>
## Realizando um ajuste para cada tratamento - Primeiro para Poisson<br>
<br>
sort(tapply(y,trat,mean))<br>
levels(trat)<br>
zi_final<-zeroinfl(y ~ trat+tempo | trat+tempo, data = dados,
link = "logit")<br>
summary(zi_final)<br>
<br>
## Aqui não consigo descobrir os coeficientes para cada tratamento
para fazer as curvas!!!!!<br>
<br>
Obrigado,<br>
<br>
Alexandre <br>
<div class="moz-cite-prefix">Em 26/09/2013 16:28, walmes . escreveu:<br>
</div>
<blockquote
cite="mid:CAFU=EkbuGMCHmCdwDtiYfCOza5hjd+uG1Zg1h1ZOZ9RUKP0y-w@mail.gmail.com"
type="cite">
<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 moz-do-not-send="true"
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
moz-do-not-send="true"
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
moz-do-not-send="true"
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
moz-do-not-send="true"
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
moz-do-not-send="true"
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 moz-do-not-send="true" href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a>
<a moz-do-not-send="true" 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 moz-do-not-send="true" 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 moz-do-not-send="true" href="tel:%28%2B55%29%2065%208132-8112" value="+556581328112" target="_blank">(+55) 65 8132-8112</a> (TIM) <a moz-do-not-send="true" href="tel:%28%2B55%29%2065%209686-6970" value="+556596866970" target="_blank">(+55) 65 9686-6970</a> (VIVO)
<a moz-do-not-send="true" href="mailto:e-mails:alexandresantosbr@yahoo.com.br" target="_blank">e-mails:alexandresantosbr@yahoo.com.br</a>
<a moz-do-not-send="true" href="mailto:alexandre.santos@cas.ifmt.edu.br" target="_blank">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a moz-do-not-send="true" 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 moz-do-not-send="true"
href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a moz-do-not-send="true"
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 moz-do-not-send="true"
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>
<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>