stata如何显示累积发生率图的number at risk

可以用曲线下面积代表危险程度。

stata不能直接导出含有number at risk的KM图，需要自己把表格贴在KM图下面；用R就可以，survival和ggkm包就可以，代码如下：

library(survival)

library(ggkm)

data(colon)

attach(colon)

fit <- survfit(Surv(time,status)~rx, data=colon)

ggkm(fit,table = FALSE)

ggkm(fit,table = TRUE) #不生成p-value

ggkm(fit,table = TRUE,main = "Survival curve for those that do adhere",pval = TRUE,ci=TRUE) #生成p-value