Seltsame Inkonsistenz zwischen JMP (SAS) und R bei der Abschätzung der mittleren ÜberlebenszeitHallo allerseits,
Ich bin auf die folgende seltsame Inkonsistenz zwischen R und JMP (oder SAS) bei der mittleren Überlebenszeitschätzung gestoßen.
Ich verstehe, dass CI aufgrund unterschiedlicher Vorgehensweisen unterschiedlich sein kann, aber ich hätte nie gedacht, dass die Schätzung des mittleren Überlebenspunkts anders sein könnte. Hat hier jemand zufällig den tatsächlichen technischen Grund? Ich versuchte es herauszufinden, hatte aber kein Glück.
Für den folgenden einfachen Datensatz (Ereignis = 1 bedeutet Tod, 0 bedeutet Zensur),
Zeit | Veranstaltung | 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 kann keine Schätzung der mittleren Überlebenszeit für die Gruppe trt = 1 liefern, während R can kann.
JMP-Ausgabe: (JMP 12.1, Analyse -> Zuverlässigkeit und Überleben -> Überleben, Y = Zeit, Gruppierung = Trt, Cencor = Ereignis, Zensurcode = 0)
| Gruppe | Mittlere Zeit | ... |
| 0 | 116 | ... |
| 1 | . | ... |
R Ausgabe: (Code: Bibliothek (Überleben); Übersicht (survfit (Surv (Zeit, Ereignis) ~ trt, data = data)) $ table )
| | ... | Median | ... |
| trt = 0 | ... | 116 | ... |
| trt = 1 | ... | 150 | ... |
Vielen Dank im Voraus!
Re: Seltsame Inkonsistenz zwischen JMP (SAS) und R bei der Abschätzung der mittleren ÜberlebenszeitVielen Dank für die Erklärung! Das macht Sinn!
Jetzt muss ich sehen, ob ich eine Formel oder etwas finden kann, das eine klare Definition in SAS vs. R gibt, da ich R das SAS-Ergebnis ausgeben muss.
Danke noch einmal!
Beste,
Yan
Re: Seltsame Inkonsistenz zwischen JMP (SAS) und R bei der Abschätzung der mittleren ÜberlebenszeitIch glaube, es gibt zu viele rechtszensierte Beobachtungen, um den Median zu definieren. Für die 8 Beobachtungen in der Gruppe von trt = 1 sind 3 Fehler bei t = 138, 1 ist ein Fehler bei t = 150 und die Hälfte (4) von ihnen wird bei t = 150 rechtszensiert. Da der Median von 8 Werten der Durchschnitt der 4. und 5. Werte ist, müssten wir 150 und einige (unbekannte) Werte größer als 150 berechnen. Folglich meldet die Survival-Plattform einen fehlenden Wert für die Median Time.
Was R angibt, warum R einen Wert gibt, kann ich nur spekulieren, aber ich schätze, da der 4. geordnete Wert nicht zensiert wird und die vierte Beobachtung von 8 das 0,50-Quantil wäre, gibt R den Wert der 4. Beobachtung als mittlere Überlebenszeit. Wenn Sie eine weitere t = 150-Beobachtung (Zeile 10) ändern, um rechtszensiert zu werden, meldet R, dass der Median für trt = 1 fehlt (NA).
Hoffe das hilft,
Michael
