找了好久,中文的竞争风险模型的学习资料好少哦,再加上帮粉丝做了一个竞争模型的分析,今天顺带就给大家写一个竞争风险回归的例子。也是接着上一篇文章的续R数据分析:竞争风险模型的做法和解释
实例描述
我们有177例干细胞移植的急性白血病患者,我们关心这些病人的白血病复发情况,但是记住,因为白细胞移植,其移植相关死亡风险也加大。那么死亡和复发存在竞争关系,毕竟病人有可能还没等到复发就死了嘛,所以做数据分析的时候要考虑。对于这样的数据我们应该用风险竞争模型。
同时我们更关心协变量比如说病人性别,白血病类型,移植时期,干细胞来源等等这些变量是如何影响病人白血病复发风险的。
我们的数据长这样:

其中ftime为时间变量,status为删失或者竞争事件,我们要用到的预测变量包括Age, Sex, D, Phase, and Source共4个,分别代表病人的年龄,性别,疾病类型,阶段,干细胞来源。
首先我们必须将预测变量中的因子类型转化为哑变量,我们需要构建一个哑变量转化函数:
factor2ind <- function(x, baseline)
{
xname
本文介绍了如何使用R语言进行竞争风险模型分析,以研究177例干细胞移植急性白血病患者的复发风险。通过factor2ind函数处理分类变量,构建预测变量矩阵,并使用crr函数拟合模型。结果显示,疾病阶段对复发风险有显著影响,而其他变量如性别、年龄和干细胞来源的影响不显著。该模型可用于理解不同因素如何影响白血病复发的相对风险。
订阅专栏 解锁全文

5699

被折叠的 条评论
为什么被折叠?



