Как построить модель опасности Кокса с помощью сплайнов
У меня есть следующая модель:
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? Я стремлюсь к чему-то похожему на это: