Ich habe ein proportionales Cox-Gefährdungsmodell eingerichtet, das den folgenden Code in R verwendet, der die Mortalität vorhersagt. Die Kovariaten A, B und C werden hinzugefügt, um Verwechslungen zu vermeiden (d. H. Alter, Geschlecht, Rasse), aber wir sind wirklich an dem Prädiktor X interessiert. X ist eine kontinuierliche Variabl

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

Jetzt habe ich Probleme, eine Kaplan-Meier-Kurve dafür zu zeichnen. Ich habe gesucht, wie ich diese Figur erschaffe, aber ich hatte nicht viel Glück. Ich bin nicht sicher, ob ein Kaplan-Meier-Diagramm für ein Cox-Modell möglich ist. Passt sich der Kaplan-Meier meinen Kovariaten an oder braucht er sie nicht?

Was ich versucht habe ist unten, aber mir wurde gesagt, das ist nicht richtig.

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

Ich habe auch versucht, eine Zahl zu zeichnen, die die kumulative Sterblichkeitsgefahr zeigt. Ich weiß nicht, ob ich es richtig mache, da ich es auf verschiedene Arten ausprobiert habe und unterschiedliche Ergebnisse erzielt habe. Idealerweise möchte ich zwei Linien zeichnen, eine, die das Sterberisiko für das 75. Perzentil von X und eine, die das 25. Perzentil von X zeigt. Wie kann ich das tun?

Ich könnte alles auflisten, was ich versucht habe, aber ich möchte niemanden verwirren!

Danke vielmals