Ошибка "ведущий младший из% заказа не является положительно определенной" в MICE

Question

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

Ошибка в chol.default(fit.sum$cov.unscaled): ведущий несовершеннолетний порядка 17 не является положительно определенным

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

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

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

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

Answer 1

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