ANOVA_un_nivel.R

# Comparación de medias, ANOVA
# Prueba 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

test <- aov(x~f)

print(summary(test))

##              Df Sum Sq Mean Sq F value Pr(>F)    
## f             2  15.41   7.704   8.337  3e-04 ***
## Residuals   297 274.46   0.924                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1