Как построить модель опасности Кокса с помощью сплайнов

У меня есть следующая модель:

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

где длина - это длительность наблюдения (включая смертность), предиктор - предиктор смертности.

Вывод команды выше такой:

                         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

Как я могу построить эту модель так, чтобы получить красивую кривую линию с 95% доверительными полосами и коэффициентом опасности на оси Y? Я стремлюсь к чему-то похожему на это:

Ответы на вопрос(1)

Ваш ответ на вопрос