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

Question

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

抱歉，如果这个问题在某个地方得到了回答...我无法在这里 google 或 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() 可以得到组的总和，然后使用 [ ]

进行索引

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

R Survival Analysis: Show total Patients (n) for Each Group

r

survival-analysis