<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 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 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>