使用“icenReg”包在 R 中进行间隔删失生存分析

如何解决使用“icenReg”包在 R 中进行间隔删失生存分析

我正在尝试使用带有函数“ic_sp”的 R 包“icenReg”来执行区间删失生存分析。在我们的研究中，我们想比较阿尔茨海默病多基因风险高与低的组，就病理性淀粉样蛋白发现 (A+) 而言。该研究的参与者被跟踪了 0 到 5.4 年。我让审查个人在纳入时为 A+，在纳入时 l=0 和 u=age。随访结束时仍为 A- 的个体进行右删失，l=随访结束时的年龄，u=Inf。对于在研究期间成为 A+ 的个人，我在评估时设置了 l=r=age。

我的代码是这样的：

res.CoxAD <- ic_sp(cbind(l,u) ~ AD.HiLo + apoe4,data=myData,model=ph,bs_samples=500)

变量 AD.HiLO 为“HI”或“LO”，apoe4 为 0 或 1。

运行此程序，我收到一条错误消息： “if (nrow(x)

我无法理解导致此错误的原因。大多数参与者要么是从右到左要么被审查，这会是一个问题吗？ 任何反馈都非常感谢。