我研究竞争风险并使用R. 我想在Fine和Gray(1999)中使用该模型,比例风险模型用于竞争风险的子分配,JASA,94:496-509。
我找到了cmprsk包。
但是,我有一个“间隔数据”配置,每个间隔的开始时间为t0
,结束时间为t1
,t1
为退出或右侧审查在给定实体的最后一个间隔中。以下是数据集的摘录
entity t0 t1 cov
1 0 3 12
1 3 7 4
1 7 9 1
2 2 3 2
2 3 10 9
3 0 10 11
4 0 1 0
4 1 6 21
4 6 7 12
...
我没有找到如何使用cmprsk
实现它,而它是在生存包(Surv(time,time2,…)
)中实现FOR EXAMPLE。
是否可以使用cmprsk
执行此操作,还是应该转到其他程序包?
我知道有一个Stata包(stcrreg
)正在做这件事,但我更喜欢使用R。