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: