与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]从模型的输出中删除了恒定项。

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

明智的是在failcode = x项之后添加cencode = y。