全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
4217 2
2018-01-06
悬赏 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有其他的包能够实现这个目的?



二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2019-3-25 15:55:07
是否已解决?求解决方案
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2019-3-25 15:55:09
是否已解决?求解决方案
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群