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/