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

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。指定此选项应使审查时间包含在输出中。

它会打印时间、n.risk、n.event、survival、std.err、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 对象，该列仍可能不会像帮助页面记录的那样输出，但 ratehr 似乎该列的名称为“censored”。