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:[email protected]
[email protected]
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
