[R-br] Modelo Misto + LSD de Fisher
Fernando Colugnati
fcolugnati em gmail.com
Quarta Novembro 21 17:54:55 BRST 2012
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 em 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 em 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 em 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 em 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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121121/4c2eb41b/attachment.html>
Mais detalhes sobre a lista de discussão R-br