Как запустить регрессию GEE для повторных измерений с биномиальной переменной ответа?
Я пытаюсь выполнить логистическую регрессию для обобщенных оценочных уравнений (GEE). Я использую Ideo_Dich(состоит из 0 и 1) в качестве моей переменной ответа с:
sample2$Ideo_Dich <- ordered(factor(sample2$Ideo_Dich))
library (multgee)
nomLORgee(Ideo_Dich~ square+round, data=sample2,
id= Politician_ID,repeated=Country_ID)
но я получаю следующее сообщение об ошибке:
The response variable should have more than 2 categories
Мой примерный набор данных выглядит так:
Politician_ID Country_ID Ideo_Dich round square
<int> <int> <ord> <dbl> <dbl>
1 3917 1 0 0.374 -0.486
2 3921 1 0 0.682 -0.580
3 3931 1 0 0.463 -0.801
4 3932 1 0 0.00806 -0.296
5 3935 1 0 -0.250 -0.485
6 3936 1 0 0.814 -0.684
7 3937 1 0 -0.0876 -0.421
8 3942 1 0 0.630 -0.738
9 3944 1 0 0.0779 -0.499
10 3945 1 0 0.549 -1.30
Поскольку я новичок в методологиях регрессии, я хотел бы получить некоторые рекомендации по этому вопросу.
2 ответа
GEE - это модели среднего населения. Вам нужно указать только один идентификатор в модели. Если вы можете выбрать только политиков или округ в качестве идентификатора и структуры корреляции. sample2$Ideo_Dich<- factor(sample2$Ideo_Dich) библиотека (geepack) Ideo_Dich ~ geeglm(квадрат + круг, данные =sample2, id = Politician_ID, corstr= "exchengable") сводка (Ideo_Dich)
Проблема в том, что ваша переменная ответа является двоичной, а не категориальной. Эта конкретная модель (nomLORgee) предназначена для полиномиальной регрессии GEE. Таким образом, пакетов gee или geepack и биномиальной регрессии GEE должно хватить. Вышеупомянутое - хороший тому пример.