<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Seguem sugestões no CMR.<br><br><span style="font-family:monospace">#------------------------------------------------------------------<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>betas <- c(0.1,0.9,0.6,0.3,0) ## Simulando sem efeito de tempo.<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><br>#------------------------------------------------------------------<br># Ajuste do modelo.<br><br>m0 <- zeroinfl(y~trat+tempo|trat, data=da)<br>summary(m0)<br><br>length(coef(m0))==dim(vcov(m0))<br>names(coef(m0))==colnames(vcov(m0))<br><br>##-----------------------------------------------------------------------------<br>## Teste por meio de comparações entre modelos aninhados.<br><br>## Testar o efeito de tempo na parte Poisson.<br>## ~trat+tempo|trat vs ~trat|trat.<br><br>m1 <- zeroinfl(y~trat|trat, data=da)<br><br>class(m1)<br>anova.zeroinf <- function(m0, m1){<br>    ll0 <- c(logLik(m0))<br>    ll1 <- c(logLik(m1))<br>    np0 <- length(coef(m0))<br>    np1 <- length(coef(m1))<br>    dnp <- abs(np0-np1)<br>    dll <- 2*abs(ll0-ll1)<br>    pchi <- pchisq(dll, df=dnp, lower.tail=FALSE)<br>    list("twice difference em log-likelihood"=dll,<br>         "difference in number of parameters"=dnp,<br>         "Pr(>chi)"=pchi)<br>}<br><br>## Por razão de log-verossimilhanças.<br>anova.zeroinf(m0, m1)<br><br>## Usando Wald para testar a mesma hipótese.<br>names(coef(m0))<br>L <- rbind(c(0,0,0,0,1,0,0,0,0))<br><br>require(multcomp)<br><br>## Pela aproximação quadrática da log-verossimilhança.<br>summary(glht(m0, linfct=L), test=Chisqtest())<br><br>## Em modelos lineares de efeito fixo, as duas aboradagens dão mesmo<br>## resultado porque a aproximação quadrática de uma log-verossimilhança.<br></span><br></div>​<div class="gmail_default" style="font-family:trebuchet ms,sans-serif;display:inline">​À disposição.<br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif;display:inline">Walmes.​</div></div>