Ich habe folgendes Modell:

coxph(Surv(fulength, mortality == 1) ~ pspline(predictor))

wobei Fülle eine Follow-up-Dauer (einschließlich Mortalität) ist, ist Prädiktor ein Prädiktor für Mortalität.

Die Ausgabe des obigen Befehls lautet:

                         coef  se(coef) se2    Chisq DF   p    
pspline(predictor), line 0.174 0.0563   0.0562 9.52  1.00 0.002
pspline(predictor), nonl                       4.74  3.09 0.200

Wie kann ich dieses Modell so zeichnen, dass ich die schöne Kurve mit 95% Konfidenzband und Hazard Ratio auf der y-Achse erhalte? Was ich anstrebe, ist etwas Ähnliches: