生存分析：生命表中被检查事件的数量

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)

请注意，我指定审查为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对象进行计数，该列仍可能不会按照帮助页面的记录进行输出，但是似乎该列的名称为“被审查”。