Strana incoerenza tra JMP (SAS) e R nella stima del tempo di sopravvivenza medianoCiao a tutti,
Mi sono imbattuto nella seguente strana incoerenza tra R e JMP (o SAS) nella stima del tempo di sopravvivenza mediano.
Comprendo che l'IC può essere diverso a causa di diversi modi di farlo, ma non mi sarei mai aspettato che la stima del punto di sopravvivenza mediana potesse essere diversa. Qualcuno qui capita di conoscere la vera ragione tecnica? Ho provato a capirlo, ma non ho avuto fortuna.
Per il seguente set di dati semplice (evento = 1 significa morte, 0 significa censura),
tempo | evento | trt |
116 | 1 | 0 |
116 | 0 | 0 |
110 | 1 | 0 |
116 | 1 | 0 |
116 | 1 | 0 |
110 | 1 | 0 |
116 | 1 | 0 |
116 | 1 | 0 |
150 | 0 | 1 |
150 | 1 | 1 |
138 | 1 | 1 |
138 | 1 | 1 |
150 | 0 | 1 |
150 | 0 | 1 |
150 | 0 | 1 |
138 | 1 | 1 |
JMP non può fornire una stima del tempo di sopravvivenza mediano per il gruppo trt = 1 mentre R può.
Output JMP: (JMP 12.1, usa analizza -> affidabilità e sopravvivenza -> sopravvivenza, Y = tempo, raggruppamento = trt, cencor = evento, codice di censura = 0)
| Gruppo | Tempo medio | ... |
| 0 | 116 | ... |
| 1 | . | ... |
R output: (codice: libreria (sopravvivenza); sommario (tabella di sopravvivenza (Surv (tempo, evento) ~ trt, data = data)) $ )
| | ... | mediano | ... |
| trt = 0 | ... | 116 | ... |
| trt = 1 | ... | 150 | ... |
Grazie a tutti in anticipo!
Ri: Strana incongruenza tra JMP (SAS) e R nella stima del tempo di sopravvivenza medianoGrazie mille per la spiegazione! Ciò ha senso!
Ora devo vedere se riesco a trovare una formula o qualcosa che dia una definizione chiara in SAS vs R, dal momento che ho bisogno di lasciare che R produca il risultato SAS.
Grazie ancora!
Migliore,
Yan
Ri: Strana incongruenza tra JMP (SAS) e R nella stima del tempo di sopravvivenza medianoPenso che quello che sta succedendo è che ci sono troppe osservazioni troncate per definire la mediana. Per le 8 osservazioni nel gruppo trt = 1, 3 sono guasti a t = 138, 1 è un fallimento a t = 150, e metà (4) di esse sono a destra con t = 150. Poiché la mediana di 8 valori è la media del 4 ° e del 5 ° valore, avremmo bisogno di una media di 150 e di un valore (sconosciuto) maggiore di 150. Di conseguenza, la piattaforma Survival riporta un valore mancante per il Tempo medio.
Riguardo al perché R dà un valore, posso solo speculare, ma suppongo che dal momento che il 4 ° valore ordinato non è censurato e la 4a osservazione su 8 sia il quantile 0,50, R sta riportando il valore della quarta osservazione come tempo di sopravvivenza mediano. Notare che se si cambia un'altra t = 150 osservazione (riga 10) per essere correttamente censurato, R riporta la mediana per trt = 1 mancante (NA).
Spero che questo ti aiuti,
Michael
