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

January 12, 2006
By

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

Tags: , , ,

Comments are closed.