Analisi di sopravvivenza di Kaplan-Meier
| Comando: | Statistiche |
|---|
Questa procedura esegue l'analisi di sopravvivenza e genera un diagramma di sopravvivenza di Kaplan-Meier.
Negli studi clinici lo sperimentatore è spesso interessato al tempo che trascorre prima che i partecipanti a uno studio presentino un determinato evento o endpoint. Questo evento è di solito un esito clinico, come il decesso, la scomparsa di un tumore, ecc.
I partecipanti vengono seguiti a partire da un certo punto di inizio e viene registrato il tempo necessario perché si verifichi l'evento di interesse.
Di solito la fine dello studio viene raggiunta prima che tutti i partecipanti abbiano presentato questo evento, e l'esito dei pazienti rimanenti è sconosciuto. È sconosciuto anche l'esito dei partecipanti che si sono ritirati dallo studio. Per tutti questi casi viene registrato il tempo di follow-up (dati censurati).
In MedCalc questi dati possono essere analizzati mediante una tavola di sopravvivenza (life table) o una curva di Kaplan-Meier, che è il metodo più comune per descrivere le caratteristiche di sopravvivenza.
Come inserire i dati
Per poter analizzare i dati, devi inserirli nel foglio di lavoro come segue:
- in una colonna si può inserire un codice per assegnare il caso a un determinato gruppo (gruppo di studio - gruppo di controllo).
- in una seconda colonna va registrato il tempo di sopravvivenza
- in una terza colonna va registrato se il caso ha raggiunto l'endpoint (inserendo il codice 1) oppure se il tempo è censurato, cioè l'esito è sconosciuto (inserendo il codice 0);
L'ordine di queste colonne naturalmente non ha importanza. Anche le righe non devono essere ordinate in alcun modo.
| A | B | C | D | |
|---|---|---|---|---|
| GROUP | TIME | REMISSION | ||
| 1 | 1 | 10 | 1 | |
| 2 | 2 | 1 | 1 | |
| 3 | 1 | 9 | 0 | |
| 4 | 2 | 4 | 1 | |
| 5 | 2 | 2 | 1 | |
| 6 | 1 | 34 | 0 | |
| 7 | 1 | 35 | 0 | |
| 8 | 2 | 2 | 1 | |
| 9 | 1 | 6 | 1 | |
| 10 | 2 | 11 | 1 | |
| 11 | 1 | 13 | 1 |
Il caso nella riga 1 apparteneva al gruppo 1 e ha raggiunto l'endpoint dopo 10 unità di tempo. Anche il caso nella riga 3 apparteneva al gruppo 1 ed è stato seguito per 9 unità di tempo. L'esito di questo caso è sconosciuto (ritirato dallo studio o fine dello studio) (dati di Freireich et al., Blood 1963; 21:699-716).
Da questi dati MedCalc può facilmente calcolare e costruire la curva di Kaplan-Meier.
Input richiesto
In questa finestra di dialogo vanno inseriti i seguenti dati:
- Tempo di sopravvivenza
Il nome della variabile che contiene il tempo per raggiungere l'evento di interesse, o il tempo di follow-up.
- Endpoint
Il nome di una variabile che contiene il codice 1 per i casi che hanno raggiunto l'endpoint, o il codice 0 per i casi che non hanno raggiunto l'endpoint, perché si sono ritirati dallo studio o perché è stata raggiunta la fine dello studio. Se i tuoi dati sono codificati in modo diverso, puoi usare lo strumento Definisci stato per ricodificarli.
- Fattore
Per Fattore seleziona una variabile qualitativa o discreta (variabile di raggruppamento, ad es. GROUP). Questo fattore qualitativo può essere costituito da codici alfanumerici o numerici. Questi codici vengono usati per suddividere i dati in diversi sottogruppi. Se vuoi studiare l'effetto di una variabile continua sul tempo di sopravvivenza, puoi convertire questa variabile continua in una variabile discreta con lo strumento Crea gruppi.MedCalc permette il confronto delle curve di sopravvivenza di un massimo di 6 sottogruppi.Se non è selezionata alcuna variabile Fattore, MedCalc mostra una sola curva di sopravvivenza (tutti i dati sono considerati appartenenti a un unico gruppo).
- Filtro
Un filtro per includere nel grafico solo un sottogruppo selezionato di casi.- Opzioni
- Tendenza lineare per i livelli del fattore: permette di verificare la presenza di una tendenza lineare tra i livelli del fattore (Altman, 1991). È appropriato se i livelli del fattore hanno un ordine naturale (ad esempio, se i codici del fattore rappresentano le dosi somministrate a gruppi diversi). Kaplan-Meier assume che i livelli del fattore siano equidistanti.
- Tempo medio di sopravvivenza ristretto: opzione per calcolare il tempo medio di sopravvivenza ristretto (RMST) ed eseguire il confronto dell'RMST tra i gruppi. Il tempo medio di sopravvivenza ristretto è il tempo medio di sopravvivenza libero da eventi fino a un punto temporale prestabilito. Equivale all'area sotto la curva di Kaplan-Meier dall'inizio dello studio fino a quel punto temporale. Il punto temporale selezionato deve trovarsi tra il primo e l'ultimo evento di ogni gruppo; in caso contrario il software lo reimposta al punto temporale più basso dell'ultimo evento tra i diversi gruppi.
- Grafico:
- Probabilità di sopravvivenza (%): traccia la probabilità di sopravvivenza (%) in funzione del tempo (curve decrescenti)
- 100 - Probabilità di sopravvivenza (%): traccia 100 - probabilità di sopravvivenza (%) in funzione del tempo (curve crescenti)
- Includi IC al 95% nel grafico: traccia l'intervallo di confidenza al 95% delle curve di sopravvivenza.
- Contrassegna dati censurati nel grafico: contrassegna i dati censurati nel grafico con una piccola linea verticale.
- Tabella del numero a rischio sotto al grafico: mostra sotto al grafico una tabella con il numero di soggetti a rischio.
Dopo aver inserito tutti i dati, fai clic su OK. MedCalc apre 2 finestre: una con i grafici di sopravvivenza e una con i risultati statistici.
Grafico
Le curve di sopravvivenza sono tracciate come una funzione a gradini, come nell'esempio seguente:
Con l'opzione "Includi IC al 95% nel grafico" selezionata, il grafico si presenta così:
Quando è selezionata l'opzione "Tabella del numero a rischio sotto al grafico", il risultato è:
Risultati
Riepilogo dei casi
Questa tabella mostra il numero di casi che hanno raggiunto l'endpoint (numero di eventi), il numero di casi che non hanno raggiunto l'endpoint (numero di casi censurati) e il numero totale di casi.
Tempo di sopravvivenza
TIME
Endpoint
REMISSION
Fattore
GROUP
Riepilogo dei casi
Numero di eventi a
Numero di casi censurati b
Fattore
N
%
N
%
Dimensione totale campione
1
9
42.86
12
57.14
21
2
21
100.00
0
0.00
21
Generale
30
71.43
12
28.57
42
a REMISSION = 1
b REMISSION = 0Sopravvivenza media e mediana
Il tempo di sopravvivenza medio e mediano sono riportati con il relativo intervallo di confidenza (IC) al 95%.
Il tempo di sopravvivenza medio è stimato come l'area sotto la curva di sopravvivenza nell'intervallo da 0 a tmax (Klein & Moeschberger, 2003).
La sopravvivenza mediana è il tempo più breve in corrispondenza del quale la probabilità di sopravvivenza scende a 0,5 (50%) o al di sotto. Se la curva di sopravvivenza non scende a 0,5 o al di sotto, il tempo mediano non può essere calcolato. Il tempo di sopravvivenza mediano e il suo IC al 95% sono calcolati secondo Brookmeyer & Crowley, 1982.
Sopravvivenza media e mediana
Fattore
Media
SE
IC al 95% per la media
Mediana
IC al 95% per la mediana
1
23.287
2.827
17.746 a 28.829
23.000
13.000 a 23.000
2
8.667
1.411
5.900 a 11.433
8.000
4.000 a 11.000
Generale
15.339
1.860
11.693 a 18.985
12.000
8.000 a 17.000
Tempo medio di sopravvivenza ristretto
Il tempo medio di sopravvivenza ristretto (RMST) è riportato con il relativo intervallo di confidenza al 95%. Se sono definiti dei gruppi, viene visualizzata una tabella con le differenze di RMST tra i gruppi, l'IC al 95% della differenza e il valore P associato (Royston & Parmar, 2013).
Tempo medio di sopravvivenza ristretto al punto temporale 18
Fattore
Media
SE
IC al 95% per la media
1
14.862
0.956
12.988 a 16.735
2
8.238
1.153
5.978 a 10.499
Generale
11.436
0.882
9.708 a 13.164
Confronto del tempo medio di sopravvivenza ristretto
Fattore
1
2
1
-
-6.6235
-9.5596 a -3.6874
P<0.00012
6.6235
3.6874 a 9.5596
P<0.0001-
Tabella di sopravvivenza
Per ogni punto temporale osservato sono elencate le proporzioni di sopravvivenza (con errore standard) di tutti i gruppi, oltre alla proporzione di sopravvivenza complessiva.
Tabella di sopravvivenza
[Mostra] Tabella di sopravvivenza [Nascondi]
Fattore 1
2
Generale
Sopravvivenza
tempoProporzione
di sopravvivenzaErrore
standardProporzione
di sopravvivenzaErrore
standardProporzione
di sopravvivenzaErrore
standard1
-
-
0.905
0.0641
0.952
0.0329
2
-
-
0.810
0.0857
0.905
0.0453
3
-
-
0.762
0.0929
0.881
0.0500
4
-
-
0.667
0.103
0.833
0.0575
5
-
-
0.571
0.108
0.786
0.0633
6
0.857
0.0764
-
-
0.714
0.0697
7
0.807
0.0869
-
-
0.690
0.0715
8
-
-
0.381
0.106
0.591
0.0764
9
-
-
-
-
-
-
10
0.753
0.0963
-
-
0.565
0.0773
11
-
-
0.286
0.0986
0.512
0.0788
12
-
-
0.190
0.0857
0.455
0.0796
13
0.690
0.107
-
-
0.426
0.0795
15
-
-
0.143
0.0764
0.398
0.0791
16
0.627
0.114
-
-
0.369
0.0784
17
-
-
0.0952
0.0641
0.341
0.0774
19
-
-
-
-
-
-
20
-
-
-
-
-
-
22
0.538
0.128
0.0476
0.0465
0.265
0.0765
23
0.448
0.135
0.000
0.000
0.189
0.0710
25
-
-
-
-
-
-
32
-
-
-
-
-
-
34
-
-
-
-
-
-
35
-
-
-
-
-
-
Endpoint: n osservato
9.0
21.0
n atteso
19.3
10.7
Osservato/Atteso
0.4675
1.9536
Confronto delle curve di sopravvivenza (test log-rank)
Scorrendo verso il basso, vedi il risultato del test log-rank per il confronto tra le due curve di sopravvivenza:
In questo esempio, 9 casi nel gruppo 1 e 21 casi nel gruppo 2 hanno presentato l'esito di interesse. La statistica chi-quadrato era 16,79 con un valore P associato inferiore a 0,0001. Si conclude quindi che, statisticamente, le due curve di sopravvivenza differiscono significativamente, ovvero che la variabile di raggruppamento ha un'influenza significativa sul tempo di sopravvivenza.
Confronto delle curve di sopravvivenza (test log-rank)
Chi-quadrato
16.7929
DF
1
Significatività
P < 0.0001
Hazard ratio con intervallo di confidenza al 95%
Se hai specificato un fattore, MedCalc calcola anche gli hazard ratio con l'intervallo di confidenza (IC) al 95%. L'hazard (rischio istantaneo) è una misura della rapidità con cui si verifica l'evento di interesse. L'hazard ratio confronta gli hazard di due gruppi.
Nell'esempio l'hazard ratio è 5,1462, quindi il rischio relativo stimato che l'evento di interesse si verifichi nel gruppo 2 è 5,1462 volte più alto che nel gruppo 1. Questo hazard ratio è significativamente diverso dal valore 1 (che corrisponde a hazard uguali), poiché l'intervallo di confidenza da 2,3506 a 11,2663 non include il valore 1.
Gli hazard ratio e gli intervalli di confidenza sono calcolati secondo Altman et al., 2000.
Nota che il calcolo dell'hazard ratio assume che il rapporto sia costante nel tempo: se le curve di sopravvivenza si incrociano, la statistica hazard ratio va quindi ignorata.
Rapporti di rischio a con un intervallo di confidenza al 95%
Fattore
1
2
1
-
5.1462
2.3506 a 11.26632
0.1943
0.08876 a 0.4254-
a Colonna/Riga
Test log-rank per la tendenza
Se si confrontano più di due curve di sopravvivenza e i gruppi hanno un ordine naturale, MedCalc può eseguire anche il test log-rank per la tendenza. Questo test verifica la probabilità che esista una tendenza nei punteggi di sopravvivenza tra i gruppi.
Letteratura
- Altman DG (1991) Practical statistics for medical research. London: Chapman and Hall.
- Altman DG, Machin D, Bryant TN, Gardner MJ (Eds) (2000) Statistics with confidence, 2nd ed. BMJ Books.
- Brookmeyer R, Crowley JA (1982) A confidence interval for the median survival time. Biometrics 38:29-41.
- Klein JP, Moeschberger ML (2003) Survival Analysis. Techniques for censored and truncated data, 2nd ed. New York: Springer Publishers.
- Royston P, Parmar MK (2013) Restricted mean survival time: an alternative to the hazard ratio for the design and analysis of randomized trials with a time-to-event
outcome. BMC Medical Research Methodology 13:152.
Vedi anche
- Filtro