[R-br] Modelo Misto + LSD de Fisher

Juliana Passos jullylp em yahoo.com.br
Quarta Novembro 21 13:24:20 BRST 2012


Olá,
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.
Os comandos que tentei são:


mod1 <- aov(resposta ~ fator + Error(sujeito/fator), data=dados)

require(nlme)
mod2 <- lme(resposta ~ fator, random = ~1|sujeito/fator,data=dados)


Usando o mod2 consigo rodar o teste de Tukey com a seguinte função:
glht(mod2,linfct=mcp(fator="Tukey"))
Porém eu preciso do LSD de Fisher. Encontrei a função LSD.test, mas ela não funciona para modelo misto

require(agricolae)

LSD.test(mod1,"fator", group=TRUE)
 

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? 
Desde já agradecida,
Juliana
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121121/96d034a0/attachment.html>


Mais detalhes sobre a lista de discussão R-br