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

Question

我想为一些感兴趣的事件构建一个生命表。请参阅下面的示例：

我正在尝试使用survivalR 包执行此操作：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、生存、std.err、下限 95% CI 和上限 95% CI。但是，它没有按预期打印 n.exit.censored。我试图弄清楚为什么它不打印数字审查列。我下载了最新版的生存包。

packageVersion("survival")

印刷

[1] ‘2.43.3’

这是否发生在其他人身上？

标签： rsurvival

Answer 1

n.exit.censored 列只应该在输出中用于计算过程数据，这不是您提供的。见?summary.survfit。您可以通过查看n.events == 0. 您显然想要的格式（更典型的人口统计或精算类型的显示）按年度间隔分组，而summary.survfit输出未分组，而是报告事件和审查的确切时间。

查看代码，我发现即使使用计数过程 Surv 对象，该列仍可能不会像帮助页面所记录的那样输出，但 ratehr 似乎该列的名称为“censored”。