«Leading минор порядка% не является положительно определенной» ошибка в MICE

Question

При запуске mice в R с достаточно большим набором данных (30 категориальных предикторов и предсказывал переменные,
п = 1000), я получаю следующее ошибка:

Error in chol.default(fit.sum$cov.unscaled) : the leading minor of order 17 is not positive definite

Я использую polyreg (полиномиальной регрессии) для большинства переменных как метод вменения. Поскольку я не знаю, почему эта ошибка возникает, я не могу опубликовать больше кода, который бы воспроизвести его.

Я нашел сообщение с аналогичной проблемой here при использовании линейной регрессии в MICE. Это сообщение не получило ответов на Cross Validated.

Возможно, это также связано с некоторой процедурой, которая становится необходимой в выборке Гиббса, как указано в сообщении this. ОП описал эту ошибку в самопрограммированной функции, используя MCMC, что связано с процедурами, которые использует MICE.

Я не знаю, что является ведущим несовершеннолетним - если бы кто-нибудь мог мне сказать, может быть, это было бы началом.

Answer 1

В моем случае модель, возможно, была уточнена. Если я уменьшу матрицу предиктора, проблема не возникнет. Один из вариантов - quickpred(), который автоматически генерирует более разреженную матрицу-предиктор на основе двумерных последовательных корреляций.