- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 0 个
- 通用积分
- 0
- 学术水平
- 0 点
- 热心指数
- 0 点
- 信用等级
- 0 点
- 经验
- 106 点
- 帖子
- 3
- 精华
- 0
- 在线时间
- 68 小时
- 注册时间
- 2013-9-24
- 最后登录
- 2022-8-19
高中生
还不是VIP/贵宾
- 威望
- 0 级
- 论坛币
 - 0 个
- 通用积分
- 0
- 学术水平
- 0 点
- 热心指数
- 0 点
- 信用等级
- 0 点
- 经验
- 106 点
- 帖子
- 3
- 精华
- 0
- 在线时间
- 68 小时
- 注册时间
- 2013-9-24
- 最后登录
- 2022-8-19
 | 开心 2018-1-13 23:55:37 |
|---|
签到天数: 2 天 连续签到: 1 天 [LV.1]初来乍到
|
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有其他的包能够实现这个目的?
|
|