R生存分析: Show total Patients (n) for Each Group

Question

编辑:添加您可以使用的公共(public)数据集。

抱歉，如果这个问题在某个地方得到了回答...我在谷歌或 YouTube 上找不到关于我的生活的示例或解释。

我正在尝试找到一种方法来显示每组的患者总数 (n)。

如果使用公共(public)数据源:

>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()，您可以获得组的总和，然后使用[]进行索引