Cómo realizar el ANOVA de Welch en R (paso a paso)

El ANOVA de Welch es una alternativa al ANOVA unidireccional típico cuando se viola el supuesto de varianzas iguales .

El siguiente ejemplo paso a paso muestra cómo realizar el ANOVA de Welch en R.

Paso 1: crear los datos

Para determinar si tres técnicas de estudio diferentes conducen a puntajes de exámenes diferentes, un profesor asigna aleatoriamente a 10 estudiantes para que usen cada técnica (Técnica A, B o C) durante una semana y luego hace que cada estudiante tome un examen de igual dificultad.

Los puntajes de los exámenes de los 30 estudiantes se muestran a continuación:

#create data frame 
df <-data. marco (grupo = rep (c (' A ', ' B ', ' C '), cada uno = 10),
                puntuación = c (64, 66, 68, 75, 78, 94, 98, 79, 71, 80,
                          91, 92, 93, 85, 87, 84, 82, 88, 95, 96,
                          79, 78, 88, 94, 92, 85, 83, 85, 82, 81))

#ver las primeras seis filas del encabezado del marco de datos
 (df)

   puntuación de grupo
1 A 64
2 A 66
3 A 68
4 A 75
5 A 78
6 A 94

Paso 2: Prueba de varianzas iguales

A continuación, podemos realizar la prueba de Bartlett para determinar si las varianzas entre cada grupo son iguales.

Si el valor p del estadístico de prueba es menor que algún nivel de significancia (como α = .05), entonces podemos rechazar la hipótesis nula y concluir que no todos los grupos tienen la misma varianza.

Para realizar la prueba de Bartlett podemos usar la función bartlett.test en base R, que usa la siguiente sintaxis:

bartlett.test (fórmula, datos)

A continuación, se explica cómo utilizar esta función en nuestro ejemplo:

#realiza la prueba de
 Bartlett. prueba (puntuación ~ grupo, datos = gl)

	Prueba de Bartlett de homogeneidad de varianzas

datos: puntuación por grupo
K-cuadrado de Bartlett = 8,1066, gl = 2, valor p = 0,01737

El valor p ( .01737 ) de la prueba de Bartlett es menor que α = .05, lo que significa que podemos rechazar la hipótesis nula de que cada grupo tiene la misma varianza.

Por lo tanto, se viola el supuesto de varianzas iguales y podemos proceder a realizar el ANOVA de Welch.

Paso 3: Realice el ANOVA de Welch

Para realizar el ANOVA de Welch en R, podemos usar la función oneway.test () de la base R de la siguiente manera:

#Realice el ANOVA de Welch en un solo sentido
 . prueba (puntuación ~ grupo, datos = gl, var. igual = FALSO )

	Análisis unidireccional de medias (sin asumir varianzas iguales)

datos: puntuación y grupo
F = 5.3492, núm gl = 2.00, denom gl = 16.83, valor p = 0.01591

El valor p general ( .01591 ) de la tabla ANOVA es menor que α = .05, lo que significa que podemos rechazar la hipótesis nula de que los puntajes del examen son iguales entre las tres técnicas de estudio.

Luego podemos realizar una prueba post-hoc para determinar qué medias de grupo son diferentes. Consulte los siguientes tutoriales para ver cómo realizar varias pruebas post-hoc en R:

• Cómo realizar una corrección de Bonferroni en R
• Cómo realizar la prueba de Tukey en R
• Cómo realizar la prueba de Scheffe en R

Consulte este tutorial para determinar qué prueba post-hoc es mejor utilizar según su situación.

Recursos adicionales

Cómo realizar un ANOVA unidireccional en R
Cómo realizar un ANOVA bidireccional en R
Cómo realizar un ANOVA de medidas repetidas en R

• https://r-project.org
• https://www.python.org/
• https://www.stata.com/