Caros membros,
Tentei criar três 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:[email protected]
[email protected]
Lattes:http://lattes.cnpq.br/1360403201088680
======================================================================
Em 04/11/2013 16:15, ASANTOS escreveu:
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
[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.
_______________________________________________
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.
--
======================================================================
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.
