Tengo un siguiente modelo:

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

donde la longitud es una duración del seguimiento (incluida la mortalidad), el predictor es un predictor de mortalidad.

El resultado del comando anterior es este:

                         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

¿Cómo puedo trazar este modelo para obtener una bonita línea curva con bandas de confianza del 95% y una relación de riesgo en el eje y? Lo que busco es algo similar a esto: