#r #regression #glm #multinomial
#r #регрессия #glm #многочленный
Вопрос:
Я использую большой набор данных и провел несколько логистических регрессий svyglm() . Сейчас я изучаю зависимую переменную с как минимум 5 возможными результатами. Я обнаружил svymultinom , но у меня проблемы с кодом. Я получаю следующую ошибку:
mmodel <- svymultinom(outcome3 ~ married, dataDesign)
Ошибка в UseMethod («withReplicates», design): нет применимого метода для ‘withReplicates’, применяемого к объекту класса «c (‘survey.design2’, ‘survey.design’)»
Может кто-нибудь привести пример, как использовать svymultinom ?
Ответ №1:
Похоже, что svymultinom (вы должны указать пакет, в котором он находится) работает для проектов с весом репликации и что ваш дизайн не является дизайном с весом репликации.
Я бы использовал новый svyVGAM пакет вместо,
library(svyVGAM)
mmodel <- svy_vglm(outcome3 ~ married, family=multinomial, design=dataDesign)
Вы также могли бы сделать
svymultinom(outcome3~married, design=as.svrepdesign(dataDesign))
но svyVGAM::svy_vglm решение чище
Комментарии:
1. Спасибо, доктор Ламли. Я только что просмотрел вашу документацию по svyVGAM и это именно то, что мне было нужно.
2. Я запустил свою модель, используя как svy_vglm, так и svymultinom, но результаты для них сильно отличаются. Я также провел серию отдельных регрессий svyglm, и эти результаты были похожи на svymultinom. Почему я получаю такие кардинально разные коэффициенты между svy_vglm и svymultinom?
3. Вы установили контрольные уровни для факторов одинаковыми? Возможно, вам понадобится
multinomial(refLevel=1). В противном случае вам нужно будет предоставить некоторую фактическую информацию о результатах4. Вот и все! Это было то, что я упустил из виду. Спасибо!