# Comparación de varianzas, más de dos muestras
# Test de Levene
require(car)
## Loading required package: car
set.seed(10)
x1 <- rnorm(30,10,1)
x2 <- rnorm(30,10,1.1)
x3 <- rnorm(30,10,1.2)
f <- factor(c(rep("x1",30),rep("x2",30),rep("x3",30)))
x <- c(x1,x2,x3)
var <- leveneTest(x ~ f)
print(var)
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 2 2.666 0.07521 .
## 87
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1