问题标签 [survminer]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
0 回答
27 浏览

r - 如何在 ggsurvplot 中添加特定时间点?

我想在 7 天内打破我的生存图和风险表,但后续在 90 天结束,我还想在 90 天包括一个刻度。默认情况下,当中断时间为 7 天时,情节在 84 天结束。

添加自定义 x 轴是不够的,因为它还应该反映在风险表和累积事件表中。

我已经看到可以做到这一点(附图)。有谁知道怎么做?

示例代码:

期望的输出: ggsurvplot 添加了 7 天的时间间隔和第 90 天

0 投票
1 回答
36 浏览

list - 在同一个图中组合多条 (>2) 生存曲线(空模型)

我正在尝试使用包 survminer 中的函数 ggsurvplot_combine 在同一个图上组合多个 survfit 对象。当我列出 2 个 survfit 对象时,它非常有效。但是当我以不同的方式组合 3 个 survfit 对象时,我收到错误:

级别错误 - ( tmp value = as.character(levels)): 因子级别 3 重复

我读过关于组合生存图的类似帖子(https://cran.r-project.org/web/packages/survminer/survminer.pdfhttps://github.com/kassambara/survminer/issues/195R在同一个图中绘制多条生存曲线https://rpkgs.datanovia.com/survminer/reference/ggsurvplot_combine.html)并针对此特定错误,使用“唯一”提供解决方案。但是,我什至不明白这个错误说明了哪个因素变量。我无权分享我的数据或数字,所以我会尝试复制它:

  1. 数据:
  • 时间:直到事件或随访结束之间的随访
  • 端点:1=事件,0=无事件或审查员
  1. 空模型:
  • KM1 <- survfit(Surv(data$time1,data$endpoint1)~1,
    type="kaplan-meier", conf.type="log", data=data)
  • KM2 <- survfit(Surv(data$time2,data$endpoint2)~1, type="kaplan-meier",
    conf.type="log", data=data)
  • KM3 <- survfit(Surv(data$time3,data$endpoint3)~1, type="kaplan-meier",
    conf.type="log", data=data)
  1. 列出空模型:
  • list_that_works <- list(KM1,KM3)
  • list_that_fails <- list(KM1,KM2,KM3) 似乎列表只包含两个参数: list(PFS=, OS=)
  1. 在一个图中组合 >2 个空模型:

ggsurvplot_combine(list_that_works, data=data, conf.int=TRUE, fun="event", combine=TRUE) 这给出了我正在寻找的情节,但有 2 条累积发生率曲线

ggsurvplot_combine(list_that_fails, data=data, conf.int=TRUE, fun="event", combine=TRUE) 这给出了错误'级别错误 - ( tmp value = as.character(levels)): 因子级别 3 重复' .

当我尝试使用 ggsurvplot(c(KM1,KM2,KM3), data=data, conf.int=TRUE, fun="event", combine=TRUE)组合 3 个图时,会出现错误:

错误:mutate() 'column 'survsummary' survsummary = purrr::map2(grouped.d$fit, grouped.d$name, .surv_summary, data=data' 的问题。x $ 运算符对于原子向量无效。

非常感谢任何帮助!也非常欢迎另一种组合生存拟合的方法!我最好的选择是它与只包含两个参数的“list”函数有关:list(PFS=, OS=)