Caros membros, Tentei criar dois tratamentos diferentes simulados através de rzipois() do pacote VGAM, mas estranhamente o modelo nulo não é diferente do modelo completo, o que pode estar acontecendo? Segue CRM: #------------------------------------------------------------------ # Definições da sessão. rm(list=ls()) require(pscl) require(VGAM) require(multcomp) require(lattice) require(latticeExtra) #------------------------------------------------------------------ # Dados artificiais. trat <- gl(3,100) ## 3 tratamentos tempo<- rep(sort(rep(1:10,10)),3) ## 10 tempos #Simulação de 3 distribuições inflacionadas de Poisson usando pacote VGAM lambda = 10 phi = 0.1 y1<-rzipois(100, lambda, phi) lambda = 4 phi = 0.3 y2<-rzipois(100, lambda, phi) lambda = 8 phi = 0.5 y3<-rzipois(100, lambda, phi) y <- c(y1,y2,y3) da<-as.data.frame(cbind(trat,tempo,y)) 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) ## Testando nulo e o completo compl.mod <- zeroinfl(y~trat+tempo|trat, data=da) null.mod <- update(compl.mod, . ~ 1) pchisq(-2*(logLik(compl.mod)-logLik(null.mod)),df=2,lower.tail=FALSE) ## A não significância indica que os tratamentos são iguais, por que? 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 ====================================================================== Em 03/11/2013 23:46, walmes . escreveu:
Eu não corri o CMR mas a princípio está correto. Só o grau de liberdade do teste de razão de verossimilhanças que deveria ser 1 e não 3, certo? Porque a diferença entre os espaços dos modelos é de 1 parâmetro.
À disposição. Walmes.
========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 skype: walmeszeviani homepage: http://www.leg.ufpr.br/~walmes <http://www.leg.ufpr.br/%7Ewalmes> linux user number: 531218 ==========================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br 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.