楼主: 腾飞flynn
4216 2

[问答] 如何使用R语言cmprsk包做以年龄为time scale的竞争风险生存分析 [推广有奖]

  • 0关注
  • 0粉丝

高中生

72%

还不是VIP/贵宾

-

威望
0
论坛币
0 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
106 点
帖子
3
精华
0
在线时间
68 小时
注册时间
2013-9-24
最后登录
2022-8-19

楼主
腾飞flynn 发表于 2018-1-6 14:40:35 |AI写论文
5论坛币
各位大神~
我们希望用5年的生存分析数据来进行life time 风险的预测,需要求stroke的风险函数(以death为竞争风险事件),并年龄作为time scale(常规的分析使用生存时间来分析)。
我们希望将出组年龄fage作为发生事件的时间,将入组年龄age作为进入模型的时间,不知道在才cmprsk里面怎么实现?


ID

age

py

rpy

fstatus

treatment group

fage

1

66.58726899

5.054072553

5.054072553

0

2

71.64134155

7

59.8384668

3.02532512

3.02532512

0

2

62.86379192

8

65.99589322

5.056810404

5.056810404

1

2

71.05270363

9

70.68309377

5.054072553

5.054072553

0

2

75.73716632

10

60.03832991

5.054072553

5.054072553

1

1

65.09240246

11

72.60780287

5.054072553

5.054072553

2

2

77.66187543

12

70.275154

4.783025325

4.783025325

0

1

75.05817933

……

……


……

……



ID

入组年龄

death终点生存时间

stroke终点生存时间

stroke=1
  death=2
  删失=0

用药组别:
  1=安慰剂
  2=干预组

出组年龄



use (a)
mcrr  <-  crr(fage,fstatus, as.factor(treatment group))
在cmprsk里面只能设置出组年龄为ftime,不能考虑入组年龄age的设置?

mcoxph <- coxph(Surv(age, fage, fstatus==1) ~ as.factor(treatment group) , data=a)
用coxph 函数虽然能够设置入组年龄、出组年龄为fstart和fstop,但模型只考虑了fstatus==1即stroke的风险,其他的1和0编码都默认为删失事件,算出来的不是竞争风险。

我的问题是:
想问问各位大神,用crr函数的时候,有没有办法像coxph一样,在模型里面将入组年龄age设置为进入模型的时间?
或者说在coxph函数基础上能不能进行竞争事件的设置,来做竞争风险分析?
或者说survial有其他的包能够实现这个目的?



关键词:R语言 生存分析 cmprsk 竞争风险

沙发
Arkxiao 发表于 2019-3-25 15:55:07
是否已解决?求解决方案

藤椅
Arkxiao 发表于 2019-3-25 15:55:09
是否已解决?求解决方案

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2026-2-7 21:37