Я использую R Studio и у меня есть данные опроса, в котором используются страты, веса и вложения.Соответственно, я использовал library(survey) для настройки svydesign, а затем использовал svyglm для завершения желаемой логистической регрессии, используя указанный svydesign (сохраненный как output2).
Когда я пытаюсь использовать library(epiDisplay) или любой другойиз других эпидемиологических пакетов на R я получаю следующую ошибку:
logistic.display (output2) Ошибка в logistic.display (output2): модель не из логистической регрессии
Пакет корректно работает для результатов регрессии, которые использовали glm, но не svyglm.Как мне включить мой план обследования и при этом получить эпидемиологические результаты?Доступен ли пакет, который может выполнять обе функции?
Редактировать
Пример со свободными данными (исключительно для целей примера): NHANES 2011-2012 Бромированные антипирены (BFR) - Сыворотка - Объединенные образцы(BFRPOL_G).
Доступно для скачивания здесь: https://wwwn.cdc.gov/Nchs/Nhanes/2011-2012/BFRPOL_G.XPT
После загрузки установите библиотеку (опрос) и создайте
BROM <- svydesign(id =new$SDMVPSU, strata = new$SDMVSTRA, weights =
new$WTSMSMPA, nest = TRUE, data = BFRPOL_G)
Сейчас (бессмысленно)регрессия, но, тем не менее, это демонстрирует проблему:
output2<-svyglm(RIAGENDR~RIGAGGRP+LBDBR66C, family = binomial, data=BFRPOL_G, design = BROM, maxit=100)
library(epiDisplay)
logistic.display(output2)
Результат?
> logistic.display(output2)
> Error in logistic.display(output2) : Model not from logistic regression
Здесь есть еще один вопрос / ответ по обмену стека (внизу страницы, предоставленной Эдвардом): https://stats.stackexchange.com/questions/8661/logistic-regression-in-r-odds-ratio показывает, насколько лаконично пакет epiDesign может отображать результаты.
С помощью логистической регрессии он может обеспечить грубое ИЛИ и 95% ДИ, скорректированное ИЛИ и 95% ДИ, тест Вальда и тест LR.Он также показывает количество наблюдений и AIC.
Я не знаю, как получить всю эту информацию из моей регрессии без этого пакета.Существуют и другие эпидемиологические библиотеки, но я не знаю ни одной, которая позволила бы мне включить вес / дизайн обследования.