手把手教你使用stata做竞争风险模型

在上一期我们已经介绍了怎么使用R语言做竞争风险模型,今天我们来演示一下怎么使用stata来做竞争风险模型,stata相对r语言来说的优势就是操作相对简单,大部分功能可以通过界面操作来实现,缺点就是界面相对代码不够灵活,但是适合初学者,我们继续使用上次膀胱癌bladder1的数据集来做竞争风险模型,这样可以对比一下stata做出来的结果和R语言有什么不同。
先打开stata导入数据,我们的数据不是很大,我就直接复制黏贴进去,下面就是数据的格式,
数据的名称解释为:
id: 患者编号;treatment(治疗方案):1. 安慰剂2. 维生素B6. 3.噻替帕 ;
number: 肿瘤的初始数量 size: 最大的初始肿瘤大小(厘米) recur: 复发次数
start,stop: 开始,停止:每个时间间隔的开始和结束时间 status(结局):1.存活2.复发3.死亡 rtumor: 复发时发现的肿瘤数量 rsize: 复发时最大的肿瘤大小
在这里插入图片描述
如下图箭头所示,依次点开红色箭头的菜单
在这里插入图片描述
打开下面菜单如下图所示,点击红色箭头
在这里插入图片描述
然后如下图所示填入时间变量和结局变量,点击确定
在这里插入图片描述
如下图所示,依次填入协变量和竞争时间,点击确定
在这里插入图片描述
结果就出来啦
在这里插入图片描述
还可以输入一下代码做图,图片还可以修改和美化,加入可信区间等

stcurve, cif at1(treatment=0) at2(treatment=1) at3(treatment=3)

在这里插入图片描述
从图上可以看出,3种治疗方法对患者肿瘤复发没有区别
多加几个协变量,点击确定
在这里插入图片描述
如下图所示结果和R语言做出来的一模一样,
在这里插入图片描述
Stata做出竞争风险模型非常简单,而且几乎不需要什么代码,非常适合初学者。
还可以通过stata做出列线图,这里就不多说啦,有兴趣的可以关注我的科研教程。
在这里插入图片描述

猜你喜欢

转载自blog.csdn.net/dege857/article/details/112706552