<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body 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">jullylp@yahoo.com.br</a>> escreveu:<br><br></div><blockquote type="cite"><div><div style="color:#000; background-color:#fff; font-family:verdana, helvetica, sans-serif;font-size:10pt"><div style="font-family: verdana, helvetica, sans-serif; font-size: 10pt;"><span>Olá,</span></div><div style="font-family: verdana, helvetica, sans-serif; font-size: 13px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;">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-family: verdana, helvetica, sans-serif; font-size: 13px; color: rgb(0, 0, 0); background-color: transparent; font-style: normal;">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" size="2">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="background-color: transparent; color: rgb(0, 0, 0); font-size: 13px; font-family: verdana, helvetica, sans-serif; font-style: normal;"><span style="font-family: verdana, helvetica, sans-serif; font-size: 13px; background-color: transparent;">require(nlme)</span></div><div style="background-color: transparent; color: rgb(0, 0, 0); font-size: 13px; font-family: verdana, helvetica, sans-serif; font-style: normal;"><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="background-color: transparent; color: rgb(0, 0,
0); font-size: 16px; font-family: 'Times New Roman'; font-style: normal;">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" size="2">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="background-color: transparent; color: rgb(0, 0, 0); font-size: 16px; font-family: 'Times New Roman'; font-style: normal;"><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" size="2">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><blockquote type="cite"><div><span>_______________________________________________</span><br><span>R-br mailing list</span><br><span><a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a></span><br><span><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">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">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</span></div></blockquote></body></html>