tukey_test.R

# Comparación de pares de medias, test de Tukey
# Prueba paramétrica
set.seed(10)
x1 <- rnorm(100,10)
x2 <- rnorm(100,10.5)
x3 <- rnorm(100,10)

x <- c(x1,x2,x3) # Unión de las 3 variables en una
f <- c(rep("x1",100),rep("x2",100),rep("x3",100)) # Factores


test <- TukeyHSD(aov(x~f))

print(test)

##   Tukey multiple comparisons of means
##     95% family-wise confidence level
## 
## Fit: aov(formula = x ~ f)
## 
## $f
##             diff        lwr        upr     p adj
## x2-x1  0.5415864  0.2213564  0.8618163 0.0002513
## x3-x1  0.1653008 -0.1549292  0.4855308 0.4446870
## x3-x2 -0.3762856 -0.6965156 -0.0560556 0.0164737