Конкурирующий риск регрессии с 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 срок.

Другие вопросы по тегам