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

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

0 голосов
/ 25 февраля 2019

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

enter image description here

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

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

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

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

packageVersion("survival")

Отпечатки 

[1] ‘2.43.3’

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

1 Ответ

0 голосов
/ 26 февраля 2019

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

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

...