kruskal_Wallis.R

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

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

test1 <- kruskal.test(x~f)

print(test1)

## 
##  Kruskal-Wallis rank sum test
## 
## data:  x by f
## Kruskal-Wallis chi-squared = 15.3198, df = 2, p-value = 0.0004713

# Segunda opción: data.frame
x<- data.frame(x1=x1,x2=x2,x3=x3)

test2 <- kruskal.test(x)

print(test2)

## 
##  Kruskal-Wallis rank sum test
## 
## data:  x
## Kruskal-Wallis chi-squared = 15.3198, df = 2, p-value = 0.0004713