Question

我使用R的survival包并生成一个包含不同组的幸存对象。

我知道想要访问每个组的事件数量（records - right censored）。

当我使用print(survobject)时，会有一列events显示每个组的事件数。但是，我没有成功将此列提取为矢量。

Answer 1

示例使用来自库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

Answer 2

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

获取每个幸存对象组的事件总数

2 个答案: