目录

R语言tidycmprsk包分析竞争风险模型

竞争风险模型就是指在临床事件中出现和它竞争的结局事件,这是事件会导致原有结局的改变,因此叫做竞争风险模型。比如我们想观察患者肿瘤的复发情况,但是患者在观察期突然车祸死亡,或者因其他疾病死亡,这样我们就观察不到复发情况了,这种情况下不能把缺失数据仅仅当做右删失处理,这样的话会造成数据的估值错误。这是我们应该优先选择竞争风险模型来做数据分析,而不是COX回归。竞争风险模型在数据挖掘中经常使用到。既往我们在文章《手把手教你使用R语言做竞争风险模型并绘制列线图》和《一步到位:手把手教你R语言竞争风险模型建模-列线图-校准曲线-K折验证-外部验证- 决策曲线》已经介绍了使用cmprsk包绘制和QHScrnomo包进行分析,今天咱们来介绍一下tidycmprsk包,它也是基于cmprsk包进行分析,但是功能更多,绘制图形更加精美。

下面立即开始,先导入数据和R包,使用的是tidycmprsk 包自带的trial数据

r 复制代码
library(tidycmprsk)
data(trial)

这是一个模拟两种化疗药物的模拟治疗肿瘤的数据,数据集包含200名接受药物A或药物B的患者的基线特征。数据集还包含肿瘤对治疗的反应结果。

Trt:化疗的药物,age年龄,marker标记物水平,stage阶段,grade等级,response物流反应,death_cr死亡状态,这个是结局变量,ttdeath时间变量

假设咱们想研究trt不同化疗药物对结局的影响,咱们先来个单因素的分析

r 复制代码
out<-cuminc(Surv(ttdeath, death_cr) ~ trt, trial)
out

咱们可以进一步绘图,如果没有指定,一般是默认第一个结局

r 复制代码
out %>%
  ggcuminc()

咱们也可以进一步指定,这样清楚一点

r 复制代码
out %>%
  ggcuminc(outcome = c("death from cancer", "death other causes"))

加上生存分析人数的表格

r 复制代码
out %>%
  ggcuminc(outcome = c("death from cancer")) +
  add_risktable()

也可以换个格式显示

r 复制代码
out %>%
  ggcuminc(outcome = c("death from cancer")) +
  add_risktable(
    risktable_stats = "n.risk",
    stats_label = list(n.risk = "Number at Risk"),
  )

还可以进一步绘制可信区间

r 复制代码
out %>%
  ggcuminc(outcome = c("death from cancer")) +
  add_risktable(
    risktable_stats = "n.risk",
    stats_label = list(n.risk = "Number at Risk"),
  )  +
  add_confidence_interval()

接下来咱们来演示多因素分析,如果你不指定结局的话,会默认的把第一个结局当成存活,最后一个结局当成竞争风险指标,就是咱们常见的0,1,2形式

r 复制代码
crr_mod <- crr(Surv(ttdeath, death_cr) ~ age + trt, trial)
crr_mod

如上图结果就出来了,结果还表明了选择了癌症为结局。我们也可以指定,我觉得这样比较好点,使用failcode这个参数指定

r 复制代码
crr_mod <- crr(Surv(ttdeath, death_cr) ~ age + trt, failcode="death from cancer",trial)
crr_mod

结果表明药物B相比药物A并没有明显优势。

本文是转载文章,点击查看原文
如有侵权,请联系 xyy@jishuzhan.net 删除
相关推荐
储悠然1 分钟前
Lisp语言的物联网数据分析
开发语言·后端·golang
东方珵蕴28 分钟前
COBOL语言的折线图
开发语言·后端·golang
知识中的海王28 分钟前
js逆向入门图灵爬虫练习平台 第四题学习
开发语言·前端·javascript
光算科技28 分钟前
无限滚动(Infinite Scroll)页面谷歌不收录!必须改回分页吗?
java·开发语言
重生之我要成为代码大佬30 分钟前
从零讲透DFS-深度优先搜索-2(排序与组合)
开发语言·python·算法·深度优先遍历
双叶83637 分钟前
(51单片机)独立按键控制流水灯LED流向(独立按键教程)(LED使用教程)
c语言·开发语言·数据结构·单片机·嵌入式硬件·游戏·51单片机
褚眠莘1 小时前
Lua语言的边缘计算
开发语言·后端·golang
云之兕1 小时前
Java内存模型详解:堆、栈、方法区
java·开发语言·jvm
@蓝莓果粒茶1 小时前
LeetCode第131题_分割回文串
开发语言·前端·算法·leetcode·职场和发展·c#·.net
nlog3n1 小时前
Java 代理模式 详解
java·开发语言·代理模式