Вычисление кривой выживаемости на основе модели Кокса
В этом упражнении мы воспроизведём пример из видео, следуя шагам:
- Вычислить модель Кокса
- Определить «воображаемых пациентов»
- Вычислить кривые выживаемости
- Создать
data.frameс данными кривых выживаемости - Построить график
Сейчас мы сосредоточимся на первых трёх шагах, а два оставшихся выполним в следующем упражнении.
Пакеты survival и survminer, а также данные GBSG2 уже загружены для вас.
Это упражнение является частью курса
Анализ выживаемости в R
Интерактивное практическое упражнение
Попробуйте выполнить это упражнение, дополнив этот пример кода.
# Cox model
cxmod <- ___(Surv(time, cens) ~ horTh + tsize, data = GBSG2)
# Imaginary patients
newdat <- expand.grid(
horTh = levels(GBSG2$___),
tsize = quantile(GBSG2$___, probs = c(0.25, 0.5, 0.75)))
___(newdat) <- letters[1:6]
# Inspect newdat
newdat