¿Cómo interpreto una curva de supervivencia del modelo de riesgo de Cox?


9

¿Cómo interpreta una curva de supervivencia del modelo de riesgo proporcional de Cox?

En este ejemplo de juguete, supongamos que tenemos un modelo de riesgo proporcional de Cox ageen kidneydatos variables y generamos la curva de supervivencia.

library(survival)
fit <- coxph(Surv(time, status)~age, data=kidney)
plot(conf.int="none", survfit(fit))
grid()

ingrese la descripción de la imagen aquí

Por ejemplo, en el momento , ¿qué afirmación es verdadera? o ambos están mal?200

  • Declaración 1: nos quedará un 20% de asignaturas (por ejemplo, si tenemos personas, para el día 200 , deberíamos tener aproximadamente 200 ), 1000200200

  • Declaración 2: Para una persona determinada, él / ella tiene un posibilidades de sobrevivir en el día 200 .20%200


βTx


Tenga en cuenta que su modelo asume independencia entre los tiempos del evento.
ocram

análisis de supervivencia puede tener suposiciones de independencia
Aksakal

Parece que la pregunta es realmente sobre la codificación R en lugar de las estadísticas puras. uno necesita saber la sintaxis y las características de funciones particulares usadas en el ejemplo. Si ese es el caso, ¿no está fuera de tema de alguna manera? de lo contrario, debe explicar lo que está sucediendo a aquellos que no usan R
Aksakal

Respuestas:


5

x

h(t;x)=h0(t)eβx.
H(t;x)=0th(u;x)du=0th0(u)eβxdu=H0(t)eβx,
H0(t)=0th0(u)dux
S(t;x)=eH(t;x)=eH0eβx=S0(t)eβx
S0(t)=eH0(t)

β^S^0(t)xS^(t;x)=S^0(t)eβ^x

El cálculo de esto en R especifica el valor de sus covariables en el newdataargumento. Por ejemplo, si desea la función de supervivencia para individuos de edad = 70, en R, haga

plot(survfit(fit, newdata=data.frame(age=70)))

newdata?survfit.coxphS0(t)eβx¯


Estoy de acuerdo contigo. Esta es una respuesta bien escrita. Pido disculpas al OP por mi error y aprecio la forma en que el OP lo corrigió.
Michael R. Chernick

@ hxd1101 Después de leer la página de ayuda con survfit.coxphmás cuidado, he corregido un error en mi respuesta, vea la actualización.
Jarle Tufto

2

Tendremos un 20% de asignaturas restantes (por ejemplo, si tenemos 1000 personas, para el día 200, deberíamos tener 200 restantes). o Para una persona determinada, ¿tiene un 20% de posibilidades de sobrevivir en el día 200?

En su forma más pura, la curva de Kaplan-Meier en su ejemplo no hace ninguna de las declaraciones anteriores.

La primera declaración tiene una proyección prospectiva tendrá . La curva de supervivencia básica solo describe el pasado, su muestra. Sí, el 20% de su muestra sobrevivió en el día 200. ¿Sobrevivirá el 20% en los próximos 200 días? No necesariamente.

Para hacer esa declaración, debe agregar más suposiciones, construir un modelo, etc. El modelo ni siquiera tiene que ser estadístico en un sentido como la regresión logística. Por ejemplo, podría PDE en epidemiología, etc.

Su segunda afirmación probablemente se basa en algún tipo de suposición de homogeneidad: todas las personas son iguales.


xβTx

@ hxd1011, depende de su modelo. Si estuvieras modelando piezas de automóviles, entonces podrías asumir que son iguales. por el contrario sus fallos podrían ser correlacionados por el número de lote, entonces ellos no son lo mismo, etc
Aksakal

Edité mi pregunta para que sea más específica en el modelo de Cox, ¿todavía se aplica su respuesta en la curva Kaplan_Meier?
Haitao Du

2

S0(t)S(t)

S0(t)S(t)x=0


0

S(t)=0.2t

t1h(t)h(t)

Con respecto a los supuestos: ¿pensé que las pruebas de coeficientes habituales en un entorno de regresión de Cox suponen independencia, condicional a las covariables observadas? Incluso la estimación de Kaplan-Meier parece requerir independencia entre el tiempo de supervivencia y la censura ( referencia ). Pero podría estar equivocado, entonces las correcciones son bienvenidas

Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.