Вправа з ігноруванням цензурування
Ви влаштовуєте вечірку, і о першій ночі гості раптово починають танцювати. Вам кортить проаналізувати, як довго вони танцюватимуть, тож ви починаєте збирати дані. Проблема в тому, що з часом ви втомлюєтеся і лягаєте спати.
Ви отримуєте такі праворуч цензуровані дані про тривалість танців у dancedat:
name— ім'я вашого друга.time— праворуч цензурована тривалість танцю.obs_endвказує, чи ви спостерігали кінець танцю друга (1), чи пішли спати до того, як він перестав танцювати (0).
Вранці ви починаєте аналіз, але ще втомлені й спершу ігноруєте те, що маєте цензуровані спостереження. Потім згадуєте цей курс на DataCamp і робите все правильно.
Пакет survival уже підключено для цієї вправи.
Ця вправа є частиною курсу
Аналіз виживаності в R
Інструкції до вправи
- Оцініть функцію виживання, удаючи, що всі цензуровані спостереження є повними спостереженнями.
- Оцініть функцію виживання для цього набору даних за методом Каплана — Меєра.
- Побудуйте правильну та хибну криві виживання разом за допомогою
ggsurvplot_combine()і порівняйте їх. Зверніть увагу, чим вони відрізняються.
Інтерактивна практична вправа
Спробуйте виконати цю вправу, доповнивши цей зразок коду.
# Create dancedat data
dancedat <- data.frame(
name = c("Chris", "Martin", "Conny", "Desi", "Reni", "Phil",
"Flo", "Andrea", "Isaac", "Dayra", "Caspar"),
time = c(20, 2, 14, 22, 3, 7, 4, 15, 25, 17, 12),
obs_end = c(1, 1, 0, 1, 1, 1, 1, 1, 0, 0, 0))
# Estimate the survivor function pretending that all censored observations are actual observations.
km_wrong <- survfit(___(time) ~ 1, data = dancedat)
# Estimate the survivor function from this dataset via kaplan-meier.
km <- survfit(___(___, ___) ~ ___, data = dancedat)
# Plot the two and compare
ggsurvplot_combine(list(correct = ___, wrong = ___))