Обчислення кривої виживання з моделі Cox
У цій вправі ми відтворимо приклад із відео, дотримуючись таких кроків:
- Обчислити модель Cox
- Визначити «уявних пацієнтів»
- Обчислити криві виживання
- Створити
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