r - 创建具有响应百分比和时间的逆 Kaplan Meier 曲线

Question

我正在尝试创建患者对药物治疗做出反应所需时间的逆 KM 图。

Time    response
3   57
4   35
4   85
4   90
5   55
6   65
6   89
6   72
9   97
9   89
9   98
10  99
10  92
13  99
14  50
15  97
18  60
21  70
25  76
28  77
40  82
48  86

时间以天为单位，响应以百分比为单位。起初我以为我可以使用生存分析来尝试这个，但认为危险图会更好。我不知道该怎么做。

这是一篇已发表文章的链接，第三张图显示了这一点。我还不是 KMplots 方面的专家，但我们将不胜感激任何帮助和批评！

https://www.researchgate.net/publication/7789803_Bortezomib_therapy_alone_and_in_combination_with_dexamethasone_for_previously_untreatment_症状_multiple_myeloma

Answer 1

为了解决您的问题，我首先将您的数据重组为我习惯的生存数据。这是每个事件/审查员一行。然后我拟合生存模型并绘制 KM。

dt <- as.data.frame(matrix(c(3,57
,4,35
,4,85
,4,90
,5,55
,6,65
,6,89
,6,72
,9,97
,9,89
,9,98
,10,99
,10,92
,13,99
,14,50
,15,97
,18,60
,21,70
,25,76
,28,77
,40,82
,48,86),ncol=2,byrow = TRUE))

colnames(dt) <- c("time","response")

#translate percentage of responders at each time to number of responders if we start with a population of 10000
dt$individuals <- round(10000*sapply((1:nrow(dt)),function(x){prod(dt[1:x,"response"]/100)}))

s <- data.frame(time = with(dt,rep(time, individuals))
                ,event = 1)

library(survival)

sobj <- Surv(s$time, s$event)
fit <- survfit(sobj ~ 1)
plot(fit, fun="event")