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?

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)?

Abs



Em 21 de novembro de 2012 15:33, Edson Lira <edinhoestat@yahoo.com.br> escreveu:
No pacote agrícolae, digite:
?Kruskal
Usa LSD.

Veja se ajuda.
[. ]'s.
Edson Lira
Estatístico
Ma-Am

Em 21/11/2012, às 11:24, Juliana Passos <jullylp@yahoo.com.br> escreveu:

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

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.



--
Fernando A.B. Colugnati