Extraña inconsistencia entre JMP (SAS) y R en la estimación del tiempo de supervivencia promedioHola a todos,
Encontré la siguiente inconsistencia extraña entre R y JMP (o SAS) en la estimación del tiempo de supervivencia medio.
Entiendo que la IC puede ser diferente debido a diferentes formas de hacerlo, pero nunca esperé que la estimación del punto de supervivencia promedio pudiera ser diferente. ¿Alguien aquí por casualidad sabe la razón técnica real? Traté de averiguarlo, pero no tuve suerte.
Para el siguiente conjunto de datos simple (evento = 1 significa muerte, 0 significa censura),
hora | 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 no puede proporcionar una estimación del tiempo de supervivencia mediano para el grupo trt = 1 mientras que R puede.
Salida JMP: (JMP 12.1, uso analizar -> confiabilidad y supervivencia -> supervivencia, Y = tiempo, agrupación = trt, cencor = evento, código de censura = 0)
| Grupo | Tiempo medio | ... |
| 0 | 116 | ... |
| 1 | . | ... |
Salida R: (código: librería (supervivencia); resumen (survfit (Surv (tiempo, evento) ~ trt, datos = datos)) tabla $ )
| | ... | mediana | ... |
| trt = 0 | ... | 116 | ... |
| trt = 1 | ... | 150 | ... |
¡Gracias a todos de antemano!
Re: Extraña inconsistencia entre JMP (SAS) y R en la estimación del tiempo de supervivencia promedioMuchas gracias por la explicación! ¡Eso tiene sentido!
Ahora tengo que ver si puedo encontrar una fórmula o algo que ofrezca una definición clara en SAS vs R, ya que necesito dejar que R produzca el resultado de SAS.
¡Gracias de nuevo!
Mejor,
Yan
Re: Extraña inconsistencia entre JMP (SAS) y R en la estimación del tiempo de supervivencia promedioCreo que lo que está sucediendo es que hay demasiadas observaciones censuradas para definir la mediana. Para las 8 observaciones en el grupo trt = 1, 3 son fallas en t = 138, 1 es una falla en t = 150 y la mitad (4) de ellas están censuradas a la derecha en t = 150. Dado que la mediana de 8 valores es el promedio de los valores 4º y 5º, tendríamos que promediar 150 y algún valor (desconocido) mayor que 150. Como resultado, la plataforma Survival está reportando un valor faltante para el Tiempo medio.
En cuanto a por qué R da un valor, solo puedo especular, pero supongo que dado que el 4º valor ordenado no está censurado y la 4ª observación de 8 sería el cuantil 0,50, R informa que el valor de la 4ª observación es el valor tiempo medio de supervivencia. Tenga en cuenta que si cambia una observación t = 150 más (fila 10) para que esté censurada a la derecha, R informa que la mediana para que trt = 1 falte (NA).
Espero que esto ayude,
Miguel
