Como plotar um modelo de risco Cox com splines

Eu tenho um seguinte modelo:

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

onde fulength é uma duração de acompanhamento (incluindo mortalidade), preditor é um preditor de mortalidade.

A saída do comando acima é esta:

                         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

Como plotar esse modelo para obter uma boa linha curvilínea com faixas de confiança de 95% e taxa de risco no eixo y? O que estou buscando é algo semelhante a isto:

questionAnswers(1)

yourAnswerToTheQuestion