<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Obrigado Walmes e uma última pergunta: A representação matemática
dos meus modelos segundo a análise abaixo:<br>
<br>
<br>
Count model coefficients (poisson with log link):<br>
Estimate Std. Error z value Pr(>|z|) <br>
(Intercept) 0.870118 0.092920 9.364 < 2e-16 ***<br>
trat2 0.129159 0.029110 4.437 9.13e-06 ***<br>
tempo 0.701863 0.009766 71.871 < 2e-16 ***<br>
<br>
<br>
<br>
fica para tratamento 1: y = Intercept + e^(0.701863*tempo) <br>
e tratamento 2: y= 0.870118 + e^(0.701863*tempo) + 0.129159 ?<br>
<br>
Segue abaixo CRM para futuras consultas de membros da lista:<br>
<br>
#------------------------------------------------------------------<br>
# Definições da sessão.<br>
<br>
rm(list=ls())<br>
require(pscl)<br>
require(multcomp)<br>
require(lattice)<br>
require(latticeExtra)<br>
<br>
#------------------------------------------------------------------<br>
<br>
## Para ter curvas você precisa que tempo seja númerico com mais<br>
## de 3 níveis, caso contrário não faz sentido.<br>
<br>
#------------------------------------------------------------------<br>
# Dados artificiais.<br>
<br>
da <- expand.grid(trat=gl(4,1), tempo=1:10)<br>
X <- model.matrix(~trat+tempo, da); ncol(X)<br>
betas <- c(0.1,0.9,0.6,0.3,0.7)<br>
eta <- X%*%betas<br>
y1 <- rpois(da$trat, lambda=exp(eta))<br>
y2 <- rbinom(y1, size=1, prob=0.7)<br>
da$y <- y1*y2<br>
str(da)<br>
xyplot(y~tempo|trat, data=da, jitter.x=TRUE)<br>
<br>
#------------------------------------------------------------------<br>
# Ajuste do modelo.<br>
<br>
m0 <- zeroinfl(y~trat+tempo|trat, data=da)<br>
summary(m0)<br>
<br>
#------------------------------------------------------------------<br>
# Predição do modelo considerando as duas porções.<br>
<br>
X <- model.matrix(~trat+tempo, da)<br>
i <- grep("^count\\_", names(coef(m0)))<br>
eta <- X%*%coef(m0)[i]<br>
da$y.pois <- exp(eta)<br>
<br>
X <- model.matrix(~trat, da)<br>
i <- grep("^zero\\_", names(coef(m0)))<br>
eta <- X%*%coef(m0)[i]<br>
da$y.zero <- exp(eta)/(1+exp(eta))<br>
<br>
xyplot(y~tempo|trat, data=da, jitter.x=TRUE)+<br>
as.layer(xyplot(y.pois~tempo|trat, data=da, type="l"))+<br>
as.layer(xyplot(y.zero~tempo|trat, data=da,<br>
type="l", lty=2, lwd=2))+<br>
layer(panel.abline(h=1, lty=2))<br>
<br>
# contínua: média da contagem ~ Poisson.<br>
# tracejada: probabilidade de um zero não Poisson.<br>
# abline: linha no 1, referência.<br>
<br>
#------------------------------------------------------------------<br>
<br>
<br>
<br>
Em 30/09/2013 10:53, walmes . escreveu:<br>
<blockquote
cite="mid:CAFU=Ekb63sts=2HgG9M+vbHSBRaSPB_LU+x3OJFdTagjGnPBXw@mail.gmail.com"
type="cite">
<div dir="ltr">
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif">Alexandre, segue CMR.<br>
<br>
<span style="font-family:courier new,monospace">#------------------------------------------------------------------<br>
# Definições da sessão.<br>
<br>
rm(list=ls())<br>
require(pscl)<br>
require(multcomp)<br>
require(lattice)<br>
require(latticeExtra)<br>
<br>
#------------------------------------------------------------------<br>
<br>
## Para ter curvas você precisa que tempo seja númerico com
mais<br>
## de 3 níveis, caso contrário não faz sentido.<br>
<br>
#------------------------------------------------------------------<br>
# Dados artificiais.<br>
<br>
da <- expand.grid(trat=gl(2,1), tempo=1:10)<br>
X <- model.matrix(~trat+tempo, da); ncol(X)<br>
betas <- c(0.1,0.1,0.3)<br>
eta <- X%*%betas<br>
y1 <- rpois(da$trat, lambda=exp(eta))<br>
y2 <- rbinom(y1, size=1, prob=0.7)<br>
da$y <- y1*y2<br>
str(da)<br>
xyplot(y~tempo|trat, data=da, jitter.x=TRUE)<br>
<br>
#------------------------------------------------------------------<br>
# Ajuste do modelo.<br>
<br>
m0 <- zeroinfl(y~trat+tempo|trat, data=da)<br>
summary(m0)<br>
<br>
#------------------------------------------------------------------<br>
# Predição do modelo considerando as duas porções.<br>
<br>
X <- model.matrix(~trat+tempo, da)<br>
i <- grep("^count\\_", names(coef(m0)))<br>
eta <- X%*%coef(m0)[i]<br>
da$y.pois <- exp(eta)<br>
<br>
X <- model.matrix(~trat, da)<br>
i <- grep("^zero\\_", names(coef(m0)))<br>
eta <- X%*%coef(m0)[i]<br>
da$y.zero <- exp(eta)/(1+exp(eta))<br>
<br>
xyplot(y~tempo|trat, data=da, jitter.x=TRUE)+<br>
as.layer(xyplot(y.pois~tempo|trat, data=da, type="l"))+<br>
as.layer(xyplot(y.zero~tempo|trat, data=da,<br>
type="l", lty=2, lwd=2))+<br>
layer(panel.abline(h=1, lty=2))<br>
<br>
# contínua: média da contagem ~ Poisson.<br>
# tracejada: probabilidade de um zero não Poisson.<br>
# abline: linha no 1, referência.<br>
<br>
#------------------------------------------------------------------</span><br>
<br>
</div>
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif">À disposição.<br>
Walmes.<br>
</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: <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">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><br>
</div>
</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>