Я анализирую данные Европейского социального опроса.Из-за недостатка данных я использовал пакет amelia для вменения.Зависимое значение является порядковым с 4 значениями, и поэтому я планировал выполнить упорядоченную логистическую регрессию с помощью функции «ologit» в пакете Zelig:
z.out <- zelig(as.factor(Y) ~ X1 + X2, model = "ologit", data = ameliadata)
Этот код будет выполняться, но когда я спрашиваюдля результатов показан следующий код ошибки:
z.out:
Модель: комбинированная ошибка вменений в se [i,] <- sqrt (diag (vcovlist)[[i]])): количество заменяемых элементов не кратно длине замены </strong>
У меня есть пять отдельных вмененных наборов данных.Анализируя отдельно, я могу использовать функцию Zelig и функцию «ologit» для каждого из этих пяти элементов.Проблема возникает только тогда, когда я использую свой комбинированный объект данных amelia.Я пытался оценить разные модели с одинаковым выходом амелии, и у меня, похоже, есть проблема только с теми, которые связаны с упорядоченной регрессией.Например, «ls» -модель работает нормально, и если я изменяю зависимую переменную на дихотомическую, я также могу без проблем запустить «logit» -модель.
Поэтому мне интересно, был ли кто-нибудьв состоянии запустить "ologit" с zelig на данных Амелии ранее, или если у кого-нибудь есть идеи о том, в чем может быть проблема?Я буду очень признателен за любые идеи и предложения.Большое спасибо за ваше время и помощь.
PS: Я извиняюсь за то, что не привел пример.Я подумал, что это может не относиться к делу, поскольку проблема заключается в совместимости вмененного набора данных и пакета zelig, но если кто-то посчитает нужным соответствующий пример, я приведу его в ближайшее время.