急性白血病患者生存时间比较：6-MP 治疗与安慰剂效果分析

42 位急性白血病患者参与一项临床跟踪研究，每位患者随机地采用 6-MP 治疗或安慰剂 (placebo)，一年后停止研究，以下是缓解时间（单位是周）：

使用 6-MP 者 (21 人): 6，6，6，7，10，13，16，22，23，6+，9+，10+，11+，17+，19+，20+，25+，32+，32+, 34+，35+ 使用安慰剂者 (21 人): 1，1，2，2，3，4，4，5，5，8，8，8，8，11，11，12，12，15，17，22，23

(1) 将两个治疗组的生存函数的乘积限估计图像绘制在一起，直观比较生存分布的差异;

使用 R 语言内容：进行绘图，代码如下：

library(survival)
library(survminer)

# 创建数据框
time <- c(6,6,6,7,10,13,16,22,23,6,9,10,11,17,19,20,25,32,32,34,35,1,1,2,2,3,4,4,5,5,8,8,8,8,11,11,12,12,15,17,22,23)
status <- c(rep(1,21),rep(0,21))
treat <- c(rep('6-MP',21),rep('Placebo',21))
data <- data.frame(time, status, treat)

# 计算生存函数的乘积限估计
fit <- survfit(Surv(time, status) ~ treat, data = data)
ggsurvplot(fit, data = data, pval = TRUE, legend.title = 'Treatment',
           risk.table = TRUE, risk.table.col = 'strata', ggtheme = theme_minimal())

运行代码后，得到如下图像：

从图像中可以看出，使用 6-MP 治疗的患者生存时间更长，生存曲线高于安慰剂组，两组之间的差异显著 (p=0.004)。