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