问题描述
由于我想比较观察和治疗组与未治疗组之间的生存时间,因此我使用倾向评分加权法获得IPTW。在此之后,我想获得竞争风险的累积发生率曲线(CIF)。我知道如何获得一个非加权的CIF。但是,我找不到加权CIF的R函数。谢谢大家。
非加权CIF crr(ftime,fstatus,cov1,cov2,tf,cengroup,failcode = 1,cencode = 0, 子集,na.action = na.omit,gtol = 1e-06,maxiter = 10,init,variance = TRUE)
关于, 斯坦利
解决方法
到目前为止,还没有实现此方法的R-Package。我目前正在R-Package上进行这种工作(包括置信区间,假设检验等),但是它仍在开发中,因此我还不能共享我的代码。
您可以执行以下操作: 1.)估计参数化竞争风险模型(精细和灰色/特定因果关系...) 2.)使用模型预测每个人在所有事件时间的CIF 3.)在您感兴趣的每个组内的每个时间点(使用IPTW)进行加权均值
某些论文中介绍了此过程。例如参见Neumann (2016)