Tengo un modelo de riesgos proporcionales de Cox configurado con el siguiente código en R que predice la mortalidad. Las covariables A, B y C se agregan simplemente para evitar confusiones (es decir, edad, sexo, raza), pero estamos realmente interesados en el predictor X. X es una variable continua.

cox.model <- coxph(Surv(time, dead) ~ A + B + C + X, data = df)

Ahora, tengo problemas para trazar una curva de Kaplan-Meier para esto. He estado buscando cómo crear esta figura, pero no he tenido mucha suerte. No estoy seguro de si es posible trazar un Kaplan-Meier para un modelo de Cox. ¿El Kaplan-Meier se ajusta a mis covariables o no las necesita?

Lo que probé está a continuación, pero me han dicho que esto no está bien.

plot(survfit(cox.model), xlab = 'Time (years)', ylab = 'Survival Probabilities')

También intenté trazar una figura que muestra el riesgo acumulado de mortalidad. No sé si lo estoy haciendo bien, ya que lo intenté de diferentes maneras y obtuve resultados diferentes. Idealmente, me gustaría trazar dos líneas, una que muestre el riesgo de mortalidad para el percentil 75 de X y otra que muestre el percentil 25 de X. ¿Cómo puedo hacer esto?

¡Podría enumerar todo lo que he intentado, pero no quiero confundir a nadie!

Muchas gracias.