<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body 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:<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>
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>
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>
<br>
## E agora, como desdobrar isto? Faço um desdobramento para parte de
poisson e outra para parte binomial?<br>
<br>
<br>
<br>
<div class="moz-cite-prefix">Em 23/09/2013 09:26, walmes . escreveu:<br>
</div>
<blockquote
cite="mid:CAFU=EkbGfEP+3FFRpP-Hzc8yehS=9OXEphRb3wO0XDyhvWfDbg@mail.gmail.com"
type="cite">
<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:
(+55) 41 3361 3573</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 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>