Prueba de los rangos con signo de Wilcoxon

La prueba de rangos con signo de Wilcoxon es la alternativa al t-test de muestras dependientes, en los casos en que las muestras no siguen una distribución normal (mostrando marcada asimetría o colas), o bien, cuando hay un tamaño muestral demasiado reducido para poder determinar si realmente proceden de poblaciones normales.

Condiciones para usar la prueba de los rangos con signo de Wilcoxon

Para poder utilizar esta prueba deben cumplirse las siguientes condiciones:

  • Los datos tienen que ser dependientes.
  • Los datos tienen que ser ordinales, se tienen que poder ordenar de menor a mayor o viceversa.
  • No es necesario asumir que las muestras se distribuyen de forma normal o que proceden de poblaciones normales. Pero sea cual sea el tipo de distribución de las diferencias, tiene que ser simétrica.
  • A pesar de considerarse el equivalente no paramétrico del t-test, el Wilcoxon signed-rank test trabaja con medianas, no con medias.
  • Preferible al t-test cuando hay valores atípicos, no hay normalidad de los datos o el tamaño de las muestras es pequeño.

Procedimiento (análisis manual)

A continuación, se explica la prueba, paso a paso:

  1. Planteamiento de las hipótesis con un ejemplo
  • H0: La mediana de las diferencias de cada par de datos es cero. Mediana(diferencias)=0Mediana(diferencias)=0.
  • H1: HaHa: La mediana de las diferencias entre cada par de datos es diferente de cero. Mediana(diferencias)≠0Mediana(diferencias)≠0.
  • 2. Se calculan las diferencias y la ordenación de las observaciones:
  • Se calcula la diferencia de cada par de observaciones.
diferencias <- c(antes - despues)  rbind(antes, despues, diferencias)
##             [,1] [,2] [,3] [,4] [,5] [,6] ## antes          2    5    4    6    1    3 ## despues        5    6    2    7    1    6 ## diferencias   -3   -1    2   -1    0   -3
Pre-testPost-testSigno (+/-)Diferencia pre-post testrangos
25-134.5
56-111.5
42123.0
67-111.5
11000.0
36-134.5
Toda la información necesaria en tabla

Cálculo del estadístico W

W=min(W+,W−)

  • W+W+ = suma de los rangos con signo positivo
  • W−W− = suma de los rangos con signo negativo
sumaPositivos <- sum(tabla[tabla$signo == 1,"rangos"]) sumaNegativos <- sum(tabla[ tabla$signo == -1,"rangos"] ) W <- min(c( sumaPositivos, sumaNegativos )) W
## [1] 3

4.Cálculo de p-value

Una vez ha obtenido el valor del estadístico W, se puede calcular cuál es la probabilidad de que adquiera valores igual o más extremos que el observado.

  • Si el tamaño n < 25, se compara el valor obtenido de W con los valores de una tabla Wilcoxon. Si W cae dentro del intervalo correspondiente en la tabla para ese n, la diferencia NO es significativa.
  • Si n > 25 se puede asumir que W se distribuye de forma aproximadamente Normal, rechazando H0H0 si Z calculado es menor que el valor de Z para el αα (esta es la que utiliza la función de R).

5. Cálculo del tamaño del efecto (side effect o TE)

Parece indicarse que la forma de calcular el tamaño del efecto de un test de Wilcoxon es mediante la correlación de rangos Rank correlation.

Tabla de valores críticos de la prueba T de Wilcoxon

Para calcular la significancia de un valor T, éste debe ser menor o igual que el valor especificado en la tabla dependiendo de si se trata de una prueba unilateral o bilateral

Prueba bilateralPrueba unilateral
Signif ⇒
n ⇓
0.050.010.050.01
5****0**
60**2**
72**30
83051
95183
1083103
11105137
12137179
131792112
1421122515
1525153019
1629193523
1734234127
1840274732
1946325337
2052376043
2158426749
2265487555
2373548362
2481619169
25896810076
26987511084
271078311992
2811691130101
29126100140110
30137109151120
Tabla de valores críticos de la prueba T de wilcoxon

Tabla de valores para obtener significancia de valor Z

Para que un valor z sea significativo, éste debe ser mayor que el del percentil correspondiente a un nivel de significancia determinado, dependiendo de si se trata de una prueba de una cola (unilateral) o de dos colas (bilateral)

Prueba UnilateralPrueba bilateral
Significancia0.050.010.0010.050.010.001
Percentil959999.997.599.599.95
z1.64492.32633.09021.96002.57583.2905
Puntúa esta entrada ❤️

Deja un comentario

Tu dirección de correo electrónico no será publicada. Los campos obligatorios están marcados con *

Scroll al inicio