Comparación de funciones de supervivencia

Como los estimadores de Kaplan-Meier son estimadores de máxima verosimilitud y la teoría establece que estos estimadores son asintóticamente normales, una primera aproximación es usar dicha teoría, es decir, si de una población con función de supervivencia S(t) se extraen dos muestras independientes, y a partir de las mismas se calculan los estimadores y debido a su normalidad asintótica, las variables:

son, también asintóticamente, normales de media 0 y varianza 1. Estas variables se pueden usar para comparar para cada tiempo predeterminado, las dos estimaciones. Sin embargo, generalmente, no interesa comparar para un tiempo predeterminado (a veces sí, por ejemplo para tiempos de corte muy establecidos como supervivencia a los 5 años para el cáncer), sino comparar globalmente ambas curvas de supervivencia.

Evidentemente, se podrían comparar, por este procedimiento, punto a punto, pero este modo de proceder es poco eficiente ya que no se usan todos los datos en cada comparación y aparecen los problemas asociados a las comparaciones múltiples y, en consecuencia, se han desarrollado pruebas para realizar una única comparación global. La más popular es la conocida como prueba del log-rank y puede verse como una aplicación de la prueba de Mantel y Haenszel para tablas de contingencia estratificadas y que se puede generalizar a la comparación de r funciones de supervivencia, con un estadístico que se distribuye como una c2 con r - 1 grados de libertad.

Ejemplo 2

Compárense las funciones de supervivencia para los datos siguientes; donde, como en el ejemplo anterior, las pérdidas están indicadas con asteriscos; extraídos de dos muestras hipotéticas correspondientes a dos tratamientos distintos.

trat. A: 10, 13, 15*, 17, 19*, 20*

trat. B: 8, 10, 11, 11*, 12, 15*

Para hacerlo con el SPSS, se necesita crear un archivo con tres variables: tiempo, estado (codificado, por ejemplo, 0: pérdida, 1: muerte) y tratamiento (codificado, por ejemplo, 0: tratamiento A, 1: tratamiento B).

y el resultado

Survival Analysis for TIEMPO
Factor TRATA = A

Time Status Cumulative
Survival
Standard
Error
Cumulative
Events
Number
Remaining
10  muerte ,8333  ,1521 
13  muerte ,6667  ,1925 
15  pérdida    
17  muerte ,4444  ,2222 
19  pérdida    
20  pérdida    

Number of Cases: 6 Censored: 3 ( 50,00%) Events: 3

Survival
Time
Standard
Error
95% Confidence
Interval
Mean:
(Limited to 20) 
17  ( 13; 20 ) 
Median:  17  ( 9; 25 ) 



Survival Analysis for TIEMPO

Factor TRATA = B

Time Status Cumulative
Survival
Standard
Error
Cumulative
Events
Number
Remaining
muerte ,8333  ,1521 
10  muerte ,6667  ,1925 
11  muerte ,5000  ,2041 
11  pérdida    
12  muerte ,2500  ,2041 
15  pérdida    

Number of Cases: 6 Censored: 2 ( 50,00%) Events: 4

Survival
Time
Standard
Error
95% Confidence
Interval
Mean:
(Limited to 15) 
12  ( 10; 14 ) 
Median:  11  ( 9; 13 ) 


Survival Analysis for TIEMPO


Total Number
Events
Number
Censored
Percent
Censored
TRATA A  50,00 
TRATA B  33,00 
OVERALL  12  41,67 


Test Statistics for Equality of Survival Distributions for TRATA


Statistic df Significance
Log Rank   2,23  ,1352 


O con AlcEst para ver los intervalos de confianza

Ejemplo 3

Se quiere comparar la supervivencia a una cierta intervención quirúrgica en 3 hospitales distintos. Para ello se sigue en cada hospital a una muestra aleatoria de pacientes intervenidos (habitualmente, y para aumentar el tamaño muestral, se sigue a todos los intervenidos, considerándolos entonces una muestra de todos los potenciales enfermos) y se obtienen los siguientes tiempos, en meses:

Hosp A: 1, 3, 7, 8*, 12, 12, 15*

Hosp B: 2, 2, 3, 8, 10, 10, 12*, 15

Hosp C: 1, 1, 3, 7, 10*, 12, 12, 14, 15*

El resultado es:

Survival Analysis for TIEMPO
Factor HOSPITAL = A

Time Status Cumulative
Survival
Standard
Error
Cumulative
Events
Number
Remaining
muerte ,8571  ,1323 
muerte ,7143  ,1707 
muerte ,5714  ,1870 
pérdida    
12  muerte    
12  muerte ,1905  ,1676 
15  pérdida    

Number of Cases: 7 Censored: 2 ( 28,57%) Events: 5

Survival
Time
Standard
Error
95% Confidence
Interval
Mean:
(Limited to 15) 
( 5; 13 ) 
Median:  12  ( 8; 16 ) 


Survival Analysis for TIEMPO

Factor HOSPITAL = B

Time Status Cumulative
Survival
Standard
Error
Cumulative
Events
Number
Remaining
muerte    
muerte ,7500  ,1531 
muerte ,6250  ,1712 
muerte ,5000  ,1768 
10  muerte    
10  muerte ,2500  ,1531 
12  pérdida    
15  muerte ,0000  ,0000 

Number of Cases: 8 Censored: 1 ( 12,50%) Events: 7

Survival
Time
Standard
Error
95% Confidence
Interval
Mean:  ( 4; 12 ) 
Median:  ( 2; 14 ) 


Survival Analysis for TIEMPO

Factor HOSPITAL = C

Time Status Cumulative
Survival
Standard
Error
Cumulative
Events
Number
Remaining
muerte    
muerte ,7778  ,1386 
muerte ,6667  ,1571 
muerte ,5556  ,1656 
10  pérdida    
12  muerte    
12  muerte ,2778  ,1617 
14  muerte ,1389  ,1272 
15  pérdida    

Number of Cases: 9 Censored: 2 ( 22,22%) Events: 7

Survival
Time
Standard
Error
95% Confidence
Interval
Mean:
(Limited to 15) 
( 5; 12 ) 
Median:  12  ( 6; 18 ) 


Survival Analysis for TIEMPO


Total Number
Events
Number
Censored
Percent
Censored
HOSPITAL A  28,57 
HOSPITAL B  12,50 
HOSPITAL C   22,22 
OVERALL  24  19  20,83 


Test Statistics for Equality of Survival Distributions for TRATA


Statistic df Significance
Log Rank   ,41  ,8164 



INDICE CAPÍTULO ANTERIOR SIGUIENTE CAPÍTULO VOLVER A BIOESTADÍSTICA