Uma dúvida....faz sentido este tipo de teste após um modelo misto? Se na modelagem vc escrever o modelo por meio de contrastes adequados, ainda seria necessária esta correção?<div><br></div><div>Estas correções são baseadas no MSW (quadrados médios intra), no respectivo Grau de Liberdade que dará o GL da t para o percentil desejado. Daí minha dúvida....este MSW "vale" em modelos mistos (principalmente se estimados via REML)?</div>
<div><br></div><div>Abs<br><div><br></div></div><div class="gmail_extra"><br><br><div class="gmail_quote">Em 21 de novembro de 2012 15:33, Edson Lira <span dir="ltr"><<a href="mailto:edinhoestat@yahoo.com.br" target="_blank">edinhoestat@yahoo.com.br</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="auto"><div>No pacote agrícolae, digite:</div><div>?Kruskal<br>Usa LSD.</div><div><br></div><div>Veja se ajuda.<br>
<div>[. ]'s.</div>Edson Lira<div>Estatístico</div><div>Ma-Am</div></div><div><br>Em 21/11/2012, às 11:24, Juliana Passos <<a href="mailto:jullylp@yahoo.com.br" target="_blank">jullylp@yahoo.com.br</a>> escreveu:<br>
<br></div><div><div class="h5"><blockquote type="cite"><div><div style="font-size:10pt;font-family:verdana,helvetica,sans-serif"><div style="font-family:verdana,helvetica,sans-serif;font-size:10pt"><span>Olá,</span></div>
<div style="font-style:normal;font-size:13px;background-color:transparent;font-family:verdana,helvetica,sans-serif">eu gostaria de usar o teste de comparações múltiplas LSD de Fisher após rodar um modelo misto, porém não encontrei nenhuma função para tal.</div>
<div style="font-style:normal;font-size:13px;background-color:transparent;font-family:verdana,helvetica,sans-serif">Os comandos que tentei são:</div><div style="background-color:transparent"><div style="background-color:transparent">
<br></div><div style="background-color:transparent"><font face="verdana, helvetica, sans-serif">mod1 <- aov(resposta ~ fator + Error(sujeito/fator),
 data=dados)</font></div><div style="background-color:transparent"><span style="font-family:verdana,helvetica,sans-serif;font-size:13px;background-color:transparent"><br></span></div><div style="font-style:normal;font-size:13px;background-color:transparent;font-family:verdana,helvetica,sans-serif">
<span style="font-family:verdana,helvetica,sans-serif;font-size:13px;background-color:transparent">require(nlme)</span></div><div style="font-style:normal;font-size:13px;background-color:transparent;font-family:verdana,helvetica,sans-serif">
<span style="font-family:verdana,helvetica,sans-serif;font-size:13px;background-color:transparent">mod2 <- lme(<span style="font-size:13px">resposta ~ fator</span>, random = ~1|sujeito/fator,data=dados)</span><br></div>
<div style="background-color:transparent"><br></div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'Times New Roman'">Usando o mod2 consigo rodar o teste de Tukey com a seguinte função:</div>
<div style="background-color:transparent"><font face="verdana, helvetica, sans-serif">glht(mod2,linfct=mcp(fator="Tukey"))</font></div><div style="background-color:transparent"><span style="background-color:transparent;font-size:16px">Porém eu preciso do LSD de Fisher. Encontrei a função LSD.test, mas ela não funciona para modelo misto</span><br>
</div><div style="font-style:normal;font-size:16px;background-color:transparent;font-family:'Times New Roman'"><span style="font-family:verdana,helvetica,sans-serif;font-size:13px">require(agricolae)</span><br></div>
<div style="background-color:transparent"><font face="verdana, helvetica, sans-serif">LSD.test(mod1,"fator", group=TRUE)</font></div></div><div style="font-family:verdana,helvetica,sans-serif;font-size:10pt"></div>
<div style="font-family:verdana,helvetica,sans-serif;font-size:10pt"> </div><div style="font-family:verdana,helvetica,sans-serif;font-size:10pt"><br></div><div style="font-family:verdana,helvetica,sans-serif;font-size:10pt">
Alguém sabe me indicar alguma função que faça o teste LSD de Fisher considerando o mod1, mod2, ou outro modelo com mesmos resultados? </div><div style="font-family:verdana,helvetica,sans-serif;font-size:10pt">Desde já agradecida,</div>
<div style="font-family:verdana,helvetica,sans-serif;font-size:10pt">Juliana</div><div style="font-family:verdana,helvetica,sans-serif;font-size:10pt"><br></div></div></div></blockquote></div></div><blockquote type="cite">
<div><span>_______________________________________________</span><br><span>R-br mailing list</span><br><span><a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a></span><br><span><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a></span><br>
<span>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</span></div></blockquote></div><br>_______________________________________________<br>

R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br><br clear="all"><div><br></div>-- <br>
Fernando A.B. Colugnati<br><br><br>
</div>