Comparações múltiplas no R - Monitoria MAE5762

Para fazer comparações múltiplas no R, há várias funções disponíveis. No sistema básico, há a função TukeyHSD(), que faz as comparações múltiplas de Tukey. Para fazer comparações por Bonferroni e Scheffé, conforme descrito em Neter et al (1994), veja Multiple comparisons JIMRC.

Um pacote que permite fazer outros tipos de contrastes, utilizando uma metodologia mais robusta (através da distribuição t-multivariada), é o pacote multcomp: Multiple Tests and Simultaneous Confidence Intervals. Veja o artigo On Multiple Comparisons in R, que acompanha o pacote multcomp.

As funções básicas nesse pacote são a simint() e simtest(). A primeira faz intervalos de confiança para contrastes, e a segunda faz testes dehipótese. Vários tipos de testes são possíveis, entre eles: "Dunnett", "Tukey", "Sequen", "AVE", "Changepoint", "Williams", "Marcus", "McDermott" e "Tetrade". Também é possível especificar exatamente quais contrastes você quer obter.

Como exemplo de uso, considere o exemplo da página de ajuda da simint():

>   data(warpbreaks)
>   mcHSD <- simint(breaks ~ wool + tension, data = warpbreaks, 
                     whichf="tension", type="Tukey")
>   print(mcHSD)

        Simultaneous confidence intervals: Tukey contrasts

Call: 
simint.formula(formula = breaks ~ wool + tension, data = warpbreaks, 
    whichf = "tension", type = "Tukey")

        95 % confidence intervals

                  Estimate   2.5 % 97.5 %
tensionM-tensionL  -10.000 -19.354 -0.646
tensionH-tensionL  -14.722 -24.076 -5.368
tensionH-tensionM   -4.722 -14.076  4.632

Envie seus comentários a respeito