R生存分析：显示每组的患者总数（n）

Question

编辑：添加您可以使用的公共数据集。

很抱歉，如果在某个地方得到了解答......我无法在这里，google或youtube找到这个例子或解释。

我试图找到一种方法来显示每组患者（n）。

如果使用公共数据源：

>data(leukemia)
>leukemia

      time  status      x
      9       1     Maintained
      13      1     Maintained
      13      0     Maintained
      .       .         .
      .       .         .
      22   43      1 Nonmaintained
      23   45      1 Nonmaintained

一个简单的脚本来进行分析：

library(survival)

group <- leukemia$x

surv=Surv(leukemia$time,leukemia$status)
surv.data.group <- survfit(surv~group,type="kaplan-meier",conf.type="none")

随后是带有图例的简单绘图，以显示每个组的颜色：

plot(surv.data.group, col=c(1,2), xlab='Time (Months)', 
xmax=120,mark.time=FALSE, main="Survival Curve Example")

legend("topright", legend=unique(group), col=(1:2), lwd=0.5, bty='n')

我想要的结果是在图例中或在情节本身上显示每组患者的总数（n）。例如，图例可能如下所示：

Maintained (n=11)
Nonmaintained (n=12)

Answer 1

如果您有两个小组，则可以采用以下解决方法：

group1 <- paste("Group 1 ", "(n=", table(data$group)[1], ")", sep="")
group2 <- paste("Group 2 ", "(n=", table(data$group)[2], ")", sep="")

group <- data$group
legend("topright", legend= c(group1, group2), col=(1:2), lwd=1.0, bty='n')

使用table()，您可获得各组的总和，然后使用[ ]

进行索引