| 查看: 2755 | 回复: 1 | ||
| 【悬赏金币】回答本帖问题，作者lintengfei88将赠送您 5 个金币 | ||
lintengfei88新虫 (初入文坛)
|
[求助]
如何使用R语言cmprsk包做以age为time scale的竞争风险生存分析
|
|
|
各位大神~ 我们希望用5年的生存分析数据来进行life time 风险的预测，需要求stroke的风险函数（以death为竞争风险事件），并年龄作为time scale（常规的分析使用生存时间来分析）。 我们希望将出组年龄fage作为发生事件的时间，将入组年龄age作为进入模型的时间，不知道在才cmprsk里面怎么实现？ 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有其他的包能够实现这个目的？  火狐截图_2018-01-09T12-33-25.009Z.png |
回复此楼» 猜你喜欢
领导好安排自己做杂事。在成果上给领导挂名，请其不要经常安排自己做杂事
已经有5人回复
-
校级优秀本科毕业论文，比例多少？看指导教师资历吗？
已经有7人回复
-
博士毕业找工作选择！
已经有13人回复
-
2024国家自然科学基金资助额度
已经有5人回复
-
Optics Express一直卡在Final Decision
已经有6人回复
-
审稿人一收一拒，编辑给了修
已经有18人回复
-
特聘教授，算高级职称吗？青年基金还没结题的情况下还能报面上吗？
已经有8人回复
-
求推荐博导
已经有4人回复
-
想请问下IEEE Trans最后一轮，你们等了多久？
已经有7人回复
-
等盲审
已经有16人回复
2楼2018-01-10 08:46:44
5