Конкурирующий риск регрессии с R
Я пытаюсь рискнуть с пакетом 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)
0 ответов
Я знаю, что это старый вопрос, но он может помочь другим.
Я подозреваю, что проблема в -1
после е.
А [,-1]
требуется после ковариат. Возможно, это то, к чему вы стремились с помощью-1
? В[,-1]
удаляет постоянный член из вывода model.matrix, например
crr.matrix <- model.matrix(~ a + b + c + d + e, data=mydata)[,-1]
.
Целесообразно также добавить cencode = y
после failcode = x
срок.