Анализ выживания: количество цензурированных событий в таблице жизни

Question

Я хочу построить жизненный стол для какого-нибудь интересного события. Смотрите ниже пример:

Я пытаюсь сделать это с survival Пакет R: https://cran.r-project.org/web/packages/survival/survival.pdf.

library(survival)
surv <- Surv(veteran$time, veteran$status)
km <- survfit(surv~trt, data = veteran)
summary(km, censored = TRUE)

Обратите внимание, что я указал censored = TRUE. Указание этой опции должно включить время цензуры в выходные данные.

Он печатает время, риск, риск, выживание, стандартную ошибку, более низкий 95% CI и верхний 95% CI. Тем не менее, он не печатает n.exit.censored, как ожидалось. Я пытаюсь понять, почему не печатается столбец с цензурой. Я скачал последнюю версию пакета для выживания.

packageVersion("survival")

Печать

[1] ‘2.43.3’

Это происходит для кого-то еще?

Answer 1

Предполагается, что столбец n.exit.censored должен присутствовать только в выходных данных для подсчета данных процесса, а это не то, что вы предлагали. Увидеть ?summary.survfit, Вы можете восстановить время цензуры, посмотрев на места, где n.events == 0, Формат, который вы, очевидно, желаете (который более типичен для демографии или актуарного типа отображения), группируется по годовому интервалу, тогда как summary.survfit вывод не группируется, а скорее сообщает точное время событий и цензуры.

Глядя на код, я вижу, что даже при подсчете объекта Surv, столбец может все еще не выводиться, как описано на странице справки, но, скорее всего, столбец будет иметь имя "censored".