Цель
Я пытаюсь построить функцию риска для данных о выживании в R.
Покушение
Я использую набор данных рака легких NCCTG в пакете survival в R.
В моем понимании функция опасности имеет следующий вид:
h(t) = h0(t) * e^(b1 * x1 + ... + bn * xn)
Значения для x1, ..., xn являются независимыми переменными, полученными из данных self, и с помощью следующего кода я могу получить значения для коэффициентов b1, ..., bn.
# Sample code
library(survival)
lung$status <- lung$status - 1
coxph_mdl <- coxph( Surv(time, status ) ~ age, data = lung )
summary(coxph_mdl)
Вопрос
Как я могу построить функцию риска для группы пациентов с данными о выживаемости в R?