R: Задание случайных эффектов с помощью команды glmer

Question

Я анализирую категорические данные из вопросника, проведенного в разных школах, чтобы увидеть, какие факторы могли повлиять на ответы ученика.Поэтому я строю смешанную модель, используя команду glmer из пакета R lme4.Для каждого ответа на вопрос анкеты у меня есть шесть переменных-предикторов, и я хочу включить школу как случайный эффект таким образом, чтобы как перехват, так и наклон варьировались в зависимости от школы.Я искал долго и упорно и онлайн и в автономном режиме и обнаружил противоречивые сведения, касающиеся правильный способ кодирования для этого и, будучи R новичок, я не уверен, что правильно!Вот что я придумал (где Like - это переменная ответа):

LikeM1 <-glmer(Like ~ Treatment + Regularity + Learn + Age + Gender + Organisation_Membership_Summary + (1 + Like|School), data = MagpieData, na.action = "na.omit",family = binomial(logit))

Правильно ли я определил школу как случайный эффект, так что как перехват, так и наклон меняются в зависимости от школы или нет?Возможно, я должен упомянуть, что, будучи категориальными данными, все мои переменные являются факторами в R

Answer 1

Я думаю, что вы хотите сделать это:

LikeM1 <-glmer (как ~ Лечение + Регулярность + Обучение + Возраст + Пол + Организация_Члены сообщества + (1 | Школа) + (0 + Лечение +Регулярность + Обучение + Возраст + Пол + Organisation_Membership_Summary | School), data = MagpieData, na.action = "na.omit", family = binomial (logit)) </p>

Первая часть формулы вкруглые скобки - случайный перехват, а второй - случайный наклон. Эта ссылка дает действительно хорошее объяснение.

Answer 2

Если вы хотите, чтобы наклон и точка пересечения варьировались в зависимости от группы, общая форма: y ~ x + (1 + x | группа). В скобках 1 указывает, что перехват должен варьироваться в зависимости от группы, а x указывает, что коэффициент предиктора x должен варьироваться в зависимости от группы. У вас есть много предикторов в вашей модели. Я бы начал с одного предиктора за раз, чтобы сделать интерпретацию немного легче.