我使用 R 的生存包并生成了具有不同组的 survfit 对象。
我知道想要访问records - right censored每个组的事件数（即 ）。
当我使用时print(survobject)，有一列events显示每个组的事件数。但是，我没有成功将此列提取为向量。
问问题
2165 次
2 回答
2
作为示例使用来自库survival帮助页面的数据。print()在对象上使用函数survfit会打印所有summary.survfit()对象的一部分。
library(survival)
print(survfit( Surv(futime, fustat)~rx, data=ovarian))
Call: survfit(formula = Surv(futime, fustat) ~ rx, data = ovarian)
records n.max n.start events median 0.95LCL 0.95UCL
rx=1 13 13 13 7 638 268 NA
rx=2 13 13 13 5 NA 475 NA
如果我们保存summary()对象并查看它的结构，就会发现有一个元素table。
gg<-summary(survfit( Surv(futime, fustat)~rx, data=ovarian))
str(gg)
List of 15
$ n : int [1:2] 13 13
$ time : num [1:12] 59 115 156 268 329 431 638 353 365 464 ...
$ n.risk : num [1:12] 13 12 11 10 9 8 5 13 12 9 ...
$ n.event : num [1:12] 1 1 1 1 1 1 1 1 1 1 ...
$ n.censor : num [1:12] 0 0 0 0 0 0 0 0 0 0 ...
$ surv : num [1:12] 0.923 0.846 0.769 0.692 0.615 ...
$ type : chr "right"
$ strata : Factor w/ 2 levels "rx=1","rx=2": 1 1 1 1 1 1 1 2 2 2 ...
$ std.err : num [1:12] 0.0739 0.1001 0.1169 0.128 0.1349 ...
$ upper : num [1:12] 1 1 1 0.995 0.946 ...
$ lower : num [1:12] 0.789 0.671 0.571 0.482 0.4 ...
$ conf.type: chr "log"
$ conf.int : num 0.95
$ call : language survfit(formula = Surv(futime, fustat) ~ rx, data = ovarian)
$ table : num [1:2, 1:7] 13 13 13 13 13 13 7 5 638 NA ...
..- attr(*, "dimnames")=List of 2
.. ..$ : chr [1:2] "rx=1" "rx=2"
.. ..$ : chr [1:7] "records" "n.max" "n.start" "events" ...
- attr(*, "class")= chr "summary.survfit"
元素表包含与您获得的相同信息print()。因此event，只需选择您感兴趣的列即可访问值。
gg$table
records n.max n.start events median 0.95LCL 0.95UCL
rx=1 13 13 13 7 638 268 NA
rx=2 13 13 13 5 NA 475 NA
gg$table[,4]
rx=1 rx=2
7 5
于 2013-02-22T11:40:41.033 回答
1
str( Surv(heart$start, heart$stop, heart$event)[1:10, 1:3] )
#-----------------------------
labelled [1:10, 1:3] 0 0 0 1 0 36 0 0 0 51 ...
- attr(*, "dimnames")=List of 2
..$ : NULL
..$ : chr [1:3] "start" "stop" "status"
- attr(*, "units")= chr "Day"
这些是时代：
Surv(heart$start, heart$stop, heart$event)[1:10, "stop"]
# [Day]
# [1] 50 6 1 16 36 39 18 3 51 675
Surv(heart$start, heart$stop, heart$event)[1:10, "start"]
# [Day]
# [1] 0 0 0 1 0 36 0 0 0 51
因此，如果您想要与事件相关联的时间，您只需：
> Surv(heart$start, heart$stop, heart$event)[ , "stop"][ Surv(heart$start, heart$stop, heart$event)[ , "status"] == 1]
[Day]
[1] 50 6 16 39 18 3 675 40 85 58 153 8 81 1387 1 308 36
[18] 43 37 28 1032 51 733 219 263 72 35 852 16 77 12 100 66 5
[35] 53 3 2 40 45 996 72 9 980 285 102 188 3 61 149 343 68
[52] 2 69 584 78 32 285 68 30 90 17 2 21 96 80 334 5 110
[69] 207 186 340 165 16 21 6
或者要获得数字，只需将事件标记相加：
sum( Surv(heart$start, heart$stop, heart$event)[ , "status"] == 1)
[1] 75
并按潜在协变量制表：
> with(lung, table(Surv(time, status)[, "status"], ph.ecog ))
ph.ecog
0 1 2 3
0 26 31 6 0
1 37 82 44 1
于 2013-02-22T11:40:21.037 回答