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