Caros membros, Gostaria de fazer a junção de termos qualitativos não significativos, através de análises de contraste de modelos, para verificar a semelhança entre tratamentos, no entanto, usando modelos de regressão de poisson inflacionados por zeros, aproveitando um modelo do Walmes fiz: #------------------------------------------------------------------ # Definições da sessão. rm(list=ls()) require(pscl) require(multcomp) require(lattice) require(latticeExtra) #------------------------------------------------------------------ # Dados artificiaiscom 3 tratamentos e 10 tempos da <- expand.grid(trat=gl(3,1), tempo=1:10) X <- model.matrix(~trat+tempo, da); ncol(X) betas <- c(0.1,0.1,0.3,0.5) eta <- X%*%betas y1 <- rpois(da$trat, lambda=exp(eta)) y2 <- rbinom(y1, size=1, prob=0.7) da$y <- y1*y2 str(da) xyplot(y~tempo|trat, data=da, jitter.x=TRUE) #------------------------------------------------------------------ # Ajuste do modelo. compl.mod <- zeroinfl(y~trat+tempo|trat, data=da) summary(compl.mod) #------------------------------------------------------------------ # Contrate de modelos de poisson inflacionado com zeros sort(tapply(da$y,da$trat,mean))#Ordena as médias trat2<-da$trat levels(trat2) levels(trat2)[1]<-"T1eT3" levels(trat2)[3]<-"T1eT3" levels(trat2) reduc.mod<-zeroinfl(y~trat2+tempo|trat, data=da) # Comparando o modelo completo e o ajunte com junção dos tratamentos 1 e 3 pchisq(-2*(logLik(compl.mod)-logLik(reduc.mod)),df=3,lower.tail=FALSE) # Na minha simulação: #[1] 1 #attr(,"df") #[1] 7 #attr(,"class") #[1] "logLik" Então pergunto: T1 e T3 são iguais e podem ser representados por um único modelo, isto esta correto? Ou existe alguma outra abordagem? Obrigado, -- ====================================================================== 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) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================