| This post was kindly contributed by fernandohrosa.com.br » R - go there to comment and to read the full post. |
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():
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 | > 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 |