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

时间:2016-10-13 19:14:00

标签: r survival-analysis inverse hazard

我正在尝试创建一个逆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_untreated_symptomatic_multiple_myeloma

1 个答案:

答案 0 :(得分:-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")

enter image description here