r - 使用 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

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

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

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

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

cencode = y在术语之后添加一个是明智的failcode = x。