Una prueba de rango logarítmico es la forma más común de comparar las curvas de supervivencia entre dos grupos.
Esta prueba utiliza las siguientes hipótesis:
H 0 : No hay diferencia en la supervivencia entre los dos grupos.
H A : No es una diferencia en la supervivencia entre los dos grupos.
Si el valor p de la prueba es menor que algún nivel de significancia (por ejemplo, α = .05), entonces podemos rechazar la hipótesis nula y concluir que hay evidencia suficiente para decir que hay una diferencia en la supervivencia entre los dos grupos.
Para realizar una prueba de rango logarítmico en R, podemos usar la función survdiff () del paquete de supervivencia , que usa la siguiente sintaxis:
sobreviviente (Surv (tiempo, estado) ~ predictores, datos)
Esta función devuelve un estadístico de prueba de chi-cuadrado y un valor p correspondiente.
El siguiente ejemplo muestra cómo utilizar esta función para realizar una prueba de rango logarítmico en R.
Ejemplo: prueba de rango logarítmico en R
Para este ejemplo, usaremos el conjunto de datos ováricos del paquete de supervivencia . Este conjunto de datos contiene la siguiente información sobre 26 pacientes:
- Tiempo de supervivencia (en meses) después del diagnóstico de cáncer de ovario
- Si se censuró o no el tiempo de supervivencia
- Tipo de tratamiento recibido (rx = 1 o rx = 2)
El siguiente código muestra cómo ver las primeras seis filas de este conjunto de datos:
biblioteca (supervivencia) #ver las primeras seis filas del encabezado del conjunto de datos (ovario) futime fustat age resid.ds rx ecog.ps 1 59 1 72,3315 2 1 1 2115 1 74,4932 2 1 1 3156 1 66,4658 2 1 2 4421 0 53,3644 2 2 1 5431 1 50,3397 2 1 1 6448 0 56,4301 1 1 2
El siguiente código muestra cómo realizar una prueba de rango logarítmico para determinar si existe una diferencia en la supervivencia entre los pacientes que recibieron diferentes tratamientos:
# realizar la prueba de rango logarítmico
survdiff (Surv (futime, fustat) ~ rx, data = ovarian)
Llamada:
survdiff (fórmula = Surv (futuro, fustat) ~ rx, data = ovárico)
N Observado esperado (OE) ^ 2 / E (OE) ^ 2 / V
rx = 1 13 7 5,23 0,596 1,06
rx = 2 13 5 6,77 0,461 1,06
Chisq = 1,1 con 1 grado de libertad, p = 0,3
El estadístico de la prueba de chi-cuadrado es 1,1 con 1 grado de libertad y el valor p correspondiente es 0,3 . Dado que este valor p no es menor que .05, no rechazamos la hipótesis nula.
En otras palabras, no tenemos evidencia suficiente para decir que existe una diferencia estadísticamente significativa en la supervivencia entre los dos tratamientos.
También podemos trazar las curvas de supervivencia para cada grupo usando la siguiente sintaxis:
#plote las curvas de supervivencia para cada grupo de tratamiento
(survfit (Surv (futime, fustat) ~ rx, data = ovarian),
xlab = " Hora ",
ylab = " Probabilidad de supervivencia global ")
Podemos ver que las curvas de supervivencia son ligeramente diferentes, pero la prueba de rango logarítmico nos dijo que la diferencia no es estadísticamente significativa.
- https://r-project.org
- https://www.python.org/
- https://www.stata.com/
