<html><body><div style="color:#000; background-color:#fff; font-family:arial, helvetica, sans-serif;font-size:12pt"><div><span>Valeu Walmes, era exatamente isso que eu procurava. Agora entendi.<br></span></div><div><br><span></span></div><div><span>Muito obrigado!<br></span></div><div> </div><div style="background-color:transparent;" align="center"><font style="background-color:transparent;" face="comic sans ms" size="2"><font size="1"><b><br></b></font></font></div><div style="text-align:left;background-color:transparent;" align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif" size="2">\begin{signature}</font></font></div><div style="text-align:left;background-color:transparent;" align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif" size="2"><<>>=</font></font></div><div
style="text-align:left;background-color:transparent;" align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Prof. Dr. Ivan Bezerra Allaman</font></div><div style="text-align:left;background-color:transparent;" align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Universidade Estadual de Santa Cruz</font></div><div style="text-align:left;background-color:transparent;" align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Departamento de Ciências Exatas e Tecnológicas</font></div><div style="text-align:left;background-color:transparent;" align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">Ilhéus/BA - Brasil</font></div><div style="text-align:left;background-color:transparent;" align="center"><font class="Apple-style-span" face="courier, monaco, monospace,
sans-serif" size="2">Fone: +55 73 3680-5596</font></div><div style="text-align:left;background-color:transparent;" align="center"><font class="Apple-style-span" face="courier, monaco, monospace, sans-serif" size="2">E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com</font></div><div style="text-align:left;background-color:transparent;" align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif" size="2">@</font></font></div><div style="text-align:left;background-color:transparent;" align="center"><font style="background-color:transparent;"><font face="courier, monaco, monospace, sans-serif" size="2">\end{signature}</font></font></div><div><br></div> <div style="font-family: arial, helvetica, sans-serif; font-size: 12pt;"> <div style="font-family: times new roman, new york, times, serif; font-size: 12pt;"> <div dir="ltr"> <font face="Arial" size="2"> <hr size="1"> <b><span
style="font-weight:bold;">De:</span></b> Walmes Zeviani <walmeszeviani@gmail.com><br> <b><span style="font-weight: bold;">Para:</span></b> Ivan Bezerra Allaman <ivanalaman@yahoo.com.br> <br> <b><span style="font-weight: bold;">Enviadas:</span></b> Quinta-feira, 10 de Maio de 2012 0:14<br> <b><span style="font-weight: bold;">Assunto:</span></b> Re: [R-br] Teste de médias em GLMM com covariável!<br> </font> </div> <br><div id="yiv1478648885"><font face="trebuchet ms,sans-serif">Ivan,<br><br>Bem, de fato não estou seguro sobre ter entendido o que você quer, mas vamos lá. Partindo do fato que o nosso modelo aponta interação entre F e x nos teremos retas não paralelas (na escala do preditor linear), portando diferenças de médias não são constantes a medida que variamos o valor de x. Podemos comparar as médias em um valor fixo de x. Por exemplo, podemos comparar o peso final de suínos do sexo macho vs fêmea (F) supondo que os
animais tem mesma idade (x=100 dias). Outra possibilidade, que deve ser a que você quer, é comparar nas médias amostrais de x nos níveis de F. Assim, comparamos o peso final dos machos com idade média de 100 dias com as fêmeas com idade média final de 95 dias. Se isso é uma comparação justa ou não é outra história e depende do contexto. Eu fiz o CMR desse último caso. É isso que você procura? Eu prefiro sempre que possível evitar as comparações (principalmente seguidas de letras) e dar um gráfico com curvas ajustadas acompanhadas das bandas de confiança. Assim o leitor pode ver como são as diferenças em cada nível de x, entender a relação de x com F e apontar para quais intervalos em x existe ou se pronuncia a diferença entre os níveis de F. É sem dúvida mais informativo. Espero ter ajudado.<br>
<br><span style="font-family:courier new, monospace;">pred <- expand.grid(F=levels(da$F), x=seq(0,1,l=30))</span><br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;">pred$eta <- predict(ac0, newdata=pred)</span><br style="font-family:courier new, monospace;">
<span style="font-family:courier new, monospace;">pred$y <- predict(ac0, newdata=pred, type="response")</span><br style="font-family:courier new, monospace;"><br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;">require(lattice); require(latticeExtra)</span><br style="font-family:courier new, monospace;">
<br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;"># retas não paralelas. em que ponto vou compara-las?</span><br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;"># em cada lugar dá uma conclusão diferente!</span><br style="font-family:courier new, monospace;">
<span style="font-family:courier new, monospace;">xyplot(eta~x, groups=F, data=pred, type="l")</span><br style="font-family:courier new, monospace;"><br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;"># você pode comparar todas em um valor fixo, como x=0.5</span><br style="font-family:courier new, monospace;">
<span style="font-family:courier new, monospace;"># ou pode comparar ela no valor médio de x para cada nível de F</span><br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;">xm <- with(da, tapply(x, F, mean)); xm</span><br style="font-family:courier new, monospace;">
<span style="font-family:courier new, monospace;">mpm==0.5</span><br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;">mpm2 <- mpm</span><br style="font-family:courier new, monospace;">
<span style="font-family:courier new, monospace;">for(i in 1:length(xm)) mpm2[i,mpm2[i,]==0.5] <- xm[i]</span><br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;">mpm2</span><br style="font-family:courier new, monospace;">
<br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;">mpm2%*%betas</span><br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;">exp(mpm2%*%betas)</span><br style="font-family:courier new, monospace;">
<br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;"># diferença nas médias populacionais marginais duas à duas</span><br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;">nlevels(da$F)</span><br style="font-family:courier new, monospace;">
<span style="font-family:courier new, monospace;">comp <- t(combn(1:nlevels(da$F), 2))</span><br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;">rownames(comp) <- apply(comp, 1, paste, collapse="-")</span><br style="font-family:courier new, monospace;">
<br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;">dmpm <- t(apply(comp, 1, function(vs) mpm2[vs[1],]-mpm2[vs[2],]))</span><br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;">rownames(dmpm) <- rownames(comp)</span><br style="font-family:courier new, monospace;">
<br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;">apply(dmpm, 1, function(x) x%*%betas)</span><br style="font-family:courier new, monospace;"><br style="font-family:courier new, monospace;">
<span style="font-family:courier new, monospace;">require(multcomp)</span><br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;">summary(glht(ac0, linfct=mpm2)) # médias populacionais marginais</span><br style="font-family:courier new, monospace;">
<span style="font-family:courier new, monospace;">summary(glht(ac0, linfct=dmpm)) # contraste entre as médias populacionais marginais</span><br style="font-family:courier new, monospace;"><br style="font-family:courier new, monospace;">
<span style="font-family:courier new, monospace;">xyplot(y~x, groups=F, data=da, type="p")+</span><br style="font-family:courier new, monospace;"><span style="font-family:courier new, monospace;"> as.layer(xyplot(y~x, groups=F, data=pred, type="l"))</span><br style="font-family:courier new, monospace;">
<br>À disposição.<br>Walmes.<br><br clear="all"></font><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;">VoIP: (3361 3600) 1053 1173</span><br style="font-family:trebuchet ms, sans-serif;"><span style="font-family:trebuchet ms, sans-serif;">e-mail: <a rel="nofollow" ymailto="mailto:walmes@ufpr.br" target="_blank" href="mailto:walmes@ufpr.br">walmes@ufpr.br</a></span><br style="font-family:trebuchet ms, sans-serif;">
<span style="font-family:trebuchet ms, sans-serif;">twitter: @walmeszeviani</span><br style="font-family:trebuchet ms, sans-serif;"><span style="font-family:trebuchet ms, sans-serif;">homepage: <a rel="nofollow" target="_blank" href="http://www.leg.ufpr.br/%7Ewalmes">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><br>
</div><br><br> </div> </div> </div></body></html>