0

一些可能有用的背景。我的位置是我有两个需要合并的地块。

  1. 我使用“生存”包从一组数据中生成 Kaplan-Meier 曲线。
  2. 我还创建并使用了一个自定义函数来生成生存曲线（来自相同的数据，由 WinBUGS 处理后）

我通过以下方式生成我的 KM 曲线。

library("survminer")
library("survival")


KM.est <- survfit(Surv(IPD.data[,1],IPD.data[,2])~1,
              data=IPD.data,type="kaplan-meier")
KM.plot <- ggsurvplot(KM.est, data = IPD.data)

当我现在调用 KM.plot 时，会生成我想要的图。我还可以通过 KM.plot$plot 仅提取绘图，检查后具有以下类：

> class(KM.plot$plot)
[1] "gg"     "ggplot"

https://i.stack.imgur.com/RnSIz.png

现在考虑我的自定义函数，（注意：nma 是来自 winbugs 的自定义对象）

prepare.survival <- prep_all_survivals(nma, "Study 1", "Treatment")
surv.plot <- survival_plot(prepare.survival)

它还输出必要的情节，与

> class(surv.plot)
[1] "gg"     "ggplot"

https://i.stack.imgur.com/yASja.png

我现在的问题是，是否可以将 surv.plot 和 KM.plot 组合在一个图中？

很难重新创建整个项目，因为您将需要 WinBUGS 模型和多个专有功能，但我认为如果您能建议一种将这些对象组合在一起的方法，那么我可以重新创建它。请注意，可以通过提到的包和我在下面提供的简短数据片段来复制 KM 对象。

"54" 25.06776181 24 0 1
"55" 25.10061602 23 0 1
"56" 25.62628337 22 0 1
"57" 25.98767967 21 0 1
"58" 26.21765914 20 0 1
"59" 26.41478439 19 1 1
"60" 26.51334702 17 0 1
"61" 26.54620123 16 0 1
"62" 26.67761807 15 0 1
"63" 27.72895277 14 0 1
"64" 28.28747433 13 0 1
"65" 29.24024641 12 0 1
"66" 30.35728953 11 1 1
"67" 31.47433265 9 0 1
"68" 31.96714579 8 0 1
"69" 32.13141684 7 0 1
"70" 33.57700205 6 0 1
"71" 34.98973306 5 1 1
"72" 35.44969199 3 0 1
"73" 36.13963039 2 0 1
"74" 38.47227926 1 0 1

感谢您的时间！

4

1 回答 1

0

编辑：我刚刚意识到您可能要求将曲线组合成一个面板，在这种情况下，我认为 markus 的评论适用。

您应该可以使用以下方法执行此操作grid.arrange：

library(gridExtra)
grid.arrange(KM.plot$plot, surv.plot)

或者，如果您想要更多控制（即添加标签、标题等），您可能需要查看cowplot。

cowplot::plot_grid(KM.plot$plot, surv.plot, ncol=1)
于 2018-07-20T12:49:30.063 回答