Перший оцінювач Каплана—Майєра

У цій вправі ми використаємо ті самі дані, що й у відео. Розглянемо функцію survfit() та об'єкт, який вона створює. Ця вправа допоможе вам дослідити об'єкт survfit.

Пакет survival уже підключено для вас у цій вправі.

Ця вправа є частиною курсу

Аналіз виживаності в R

Переглянути курс

Інструкції до вправи

Дослідіть, як працює функція survfit(), ввівши ?survfit у консолі.
Обчисліть оцінювач Каплана—Майєра за допомогою survfit().
Подивіться на структуру об'єкта survfit за допомогою str().
Створіть data.frame із чотирма часовими точками, відповідною кількістю під ризиком (n.risk), кількістю спостережень з подією (n.event), кількістю цензурованих спостережень (n.censor) і значенням кривої виживання (surv). Візьміть усю інформацію з об'єкта survfit.

Інтерактивна практична вправа

Спробуйте виконати цю вправу, доповнивши цей зразок коду.

# Create time and event data
time <- c(5, 6, 2, 4, 4)
event <- c(1, 0, 0, 1, 1)

# Compute Kaplan-Meier estimate
km <- survfit(___(___, ___) ~ ___)
km

# Take a look at the structure
str(___)

# Create data.frame
data.frame(time = km$time, n.risk = ____, n.event = ____,
  n.censor = ____, surv = ___)

Редагувати та запускати код