<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
Caros membros,<br>
<br>
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:<br>
<br>
#------------------------------------------------------------------<br>
# Definições da sessão.<br>
<br>
rm(list=ls())<br>
require(pscl)<br>
require(VGAM)<br>
require(multcomp)<br>
require(lattice)<br>
require(latticeExtra)<br>
<br>
#------------------------------------------------------------------<br>
# Dados artificiais.<br>
<br>
trat <- gl(3,100) ## 3 tratamentos <br>
tempo<- rep(sort(rep(1:10,10)),3) ## 10 tempos<br>
<br>
#Simulação de 3 distribuições inflacionadas de Poisson usando pacote
VGAM<br>
lambda = 10<br>
phi = 0.1<br>
y1<-rzipois(100, lambda, phi)<br>
lambda = 4<br>
phi = 0.3<br>
y2<-rzipois(100, lambda, phi)<br>
lambda = 8<br>
phi = 0.5<br>
y3<-rzipois(100, lambda, phi)<br>
y <- c(y1,y2,y3)<br>
<br>
da<-as.data.frame(cbind(trat,tempo,y))<br>
str(da)<br>
xyplot(y~tempo|trat, data=da, jitter.x=TRUE)<br>
<br>
#------------------------------------------------------------------<br>
# Ajuste do modelo.<br>
<br>
compl.mod <- zeroinfl(y~trat+tempo|trat, data=da)<br>
summary(compl.mod)<br>
<br>
## Testando nulo e o completo <br>
compl.mod <- zeroinfl(y~trat+tempo|trat, data=da)<br>
null.mod <- update(compl.mod, . ~ 1)<br>
pchisq(-2*(logLik(compl.mod)-logLik(null.mod)),df=2,lower.tail=FALSE)
<br>
<br>
## A não significância indica que os tratamentos são iguais, por
que?<br>
<br>
<br>
Obrigado,<br>
<br>
<br>
<pre class="moz-signature" cols="72">--
======================================================================
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)
<a class="moz-txt-link-abbreviated" href="mailto:e-mails:alexandresantosbr@yahoo.com.br">e-mails:alexandresantosbr@yahoo.com.br</a>
<a class="moz-txt-link-abbreviated" href="mailto:alexandre.santos@cas.ifmt.edu.br">alexandre.santos@cas.ifmt.edu.br</a>
Lattes: <a class="moz-txt-link-freetext" href="http://lattes.cnpq.br/1360403201088680">http://lattes.cnpq.br/1360403201088680</a>
======================================================================</pre>
<br>
<br>
<br>
<br>
<br>
<br>
<div class="moz-cite-prefix">Em 03/11/2013 23:46, walmes . escreveu:<br>
</div>
<blockquote
cite="mid:CAFU=Eka=5W=-3p9jF3=BTg34VrqbK_BhoMhiSryGYMfsQYcrMA@mail.gmail.com"
type="cite">
<div dir="ltr">
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif">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.</div>
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif"><br>
</div>
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif">À disposição.</div>
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif">
Walmes.</div>
<div class="gmail_default" style="font-family:trebuchet
ms,sans-serif"><br>
</div>
<div class="gmail_extra">
<div>
<div dir="ltr"><span style="font-family:trebuchet
ms,sans-serif">==========================================================================</span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Walmes
Marques Zeviani</span><br style="font-family:trebuchet
ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">LEG
(Laboratório de Estatística e Geoinformação, 25.450418
S, 49.231759 W)</span><br style="font-family:trebuchet
ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">Departamento
de Estatística - Universidade Federal do Paraná</span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">fone:
(+55) 41 3361 3573</span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">skype:
walmeszeviani<br style="font-family:trebuchet
ms,sans-serif">
</span><span style="font-family:trebuchet ms,sans-serif">homepage:
<a moz-do-not-send="true"
href="http://www.leg.ufpr.br/%7Ewalmes"
target="_blank">http://www.leg.ufpr.br/~walmes</a></span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">linux
user number: 531218</span><br
style="font-family:trebuchet ms,sans-serif">
<span style="font-family:trebuchet ms,sans-serif">==========================================================================</span></div>
</div>
</div>
</div>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
R-br mailing list
<a class="moz-txt-link-abbreviated" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a>
<a class="moz-txt-link-freetext" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a>
Leia o guia de postagem (<a class="moz-txt-link-freetext" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</pre>
</blockquote>
<br>
<pre class="moz-signature" cols="72">
</pre>
</body>
</html>