我有两个从constsurvplot包的plot_surv_at_t生成的曲线图.一项是研究5年无复发存活率与肿瘤大小的关系.另一个是癌症特有的存活率与肿瘤大小的关系.

model1 <- coxph(Surv(Time_To_Recurrence,Recurrence) ~ Tumor_Size, data = data, x=TRUE)
model2 <- coxph(Surv(CSS,Cancer_Death) ~ Tumor_Size, data = data, x=TRUE)

plot_surv_at_t(time = "Time_To_Recurrence",
               status = "Recurrence", 
               variable = "Tumor_Size",
               data = data1, 
               model = model1, 
               t=60)

plot_surv_at_t(time = "CSS",
               status = "Cancer_Death", 
               variable = "Tumor_Size",
               data = data1, 
               model = model2, 
               t=60)

每个单独的曲线图看起来都很棒,但我想在同一曲线图中展示大小对无复发和癌症特定生存率的影响.我试过patchwork的套餐,能把它们挨着放在一起,但我还没能把它们搭在一起.

谢谢

推荐答案

一种可能的 Select 是从第二个ggplot2对象拉取数据并将其添加到第一个ggplot2对象,例如

library(survival)
library(survminer)
#> Loading required package: ggplot2
#> Loading required package: ggpubr
#> 
#> Attaching package: 'survminer'
#> The following object is masked from 'package:survival':
#> 
#>     myeloma
#options(timeout = 6000)
#install.packages("contsurvplot", dependencies = TRUE)
library(contsurvplot)

fit <- coxph(Surv(time, status) ~ age, data = lung, x = TRUE)
plot1 <- plot_surv_at_t(time = "time",
               status = "status", 
               variable = "age",
               data = lung, 
               model = fit, 
               t=60)

plot1


lung2 <- lung %>%
  mutate(CSS = sample(1:1000, 228))
fit2 <- coxph(Surv(CSS, status) ~ age, data = lung2, x = TRUE)
plot2 <- plot_surv_at_t(time = "CSS",
                        status = "status", 
                        variable = "age",
                        data = lung2, 
                        model = fit2,
                        t=60)
plot2


plot1 +
  geom_line(data = plot2$data,
            aes(x = cont, y = est),
            color = "red")

创建于2023-03-19年第reprex v2.0.2

R相关问答推荐

使用case_when和Mutate搜索多个列以寻找条件

在ubuntu 22.04上更新到R4.4后包安装出现编译错误

如果窗口在CLARME或集团之外,则有条件领先/滞后滚动总和返回NA

列出用m n个值替换来绘制n个数字的所有方法(i.o.w.:R中大小为n的集合的所有划分为m个不同子集)

更新合适的R mgcv::bam模型报告无效类型(关闭).'';错误

如何在R中合并和合并多个rabrame?

根据模式将一列拆分为多列,并在R中进行拆分

如何在格子中添加双曲曲线

`lazy_dt`不支持`dplyr/across`?

当我添加美学时,geom_point未对齐

R -在先前group_by级别汇总时获取最大大小子组的计数

根据列表中项目的名称合并数据框和列表

计算直线上点到参考点的总距离

R如何将列名转换为更好的年和月格式

使用gt_summary是否有一种方法来限制每个变量集进行配对比较?

无法将条件case_when()应用于使用!!创建的新变量Mutations

如何移动点以使它们的打印不重叠

将数据从一列转换为按组累计计数的单个虚拟变量

如何在访问之前下载的输入时同时上传和处理所有指定的shiny 输入?

从字符串01JAN2021创建日期