与R竞争风险回归

Question

我正在尝试使用cmprsk软件包运行竞争风险，但不断收到错误。我无法解决的最后一个是这一个：

Error in solve.default(h, z[[2]]) : system is computationally singular: reciprocal condition number = 3.7676e-34

这是我使用的代码。

crr.matrix <- model.matrix(~ a + b + c + d + e -1, data=mydata) crr(HV_pT1$time,HV_pT1$status,crr.matrix,failcode=2)

Answer 1

我知道这是一个古老的问题，但这可能会对其他人有所帮助。

我怀疑问题出在e后的-1。

在协变量之后 需要[,-1]。也许这就是-1的目标？ [,-1]从model.matrix的输出中删除常数项，例如

crr.matrix <- model.matrix(~ a + b + c + d + e, data=mydata)[,-1]。

明智的做法是在cencode = y期末添加一个failcode = x。