Мои данные являются порядковыми, поэтому отсутствующие значения вменяются методом polr из пакета MICE. Теперь у меня есть несколько наборов данных, на которых я могу запустить Порядковую логистическую регрессию. Но, как упоминалось в заголовке: я хочу выполнить тест Бранта, чтобы проверить предположение о параллельной регрессии. Как я могу выполнить такой тест для моих вмененных наборов данных?
olr <- with(imputed, polr(target ~ var1+var2))
olrsummary <- summary(pool(olr))
> brant(olr)
Error in formula.default(model) : invalid formula
> brant(olrsummary)
Error in temp.data[, name] : incorrect number of dimensions
Я знаю, что могу взять первый набор данных с полным (вмененный, 1) и использовать его для моего теста Бранта. Но это просто не так.
Заранее спасибо