Я применил линейную модель смешанного эффекта, используя lme4, однако это иерархическая структура данных, и у меня есть два разных набора данных с уровня 1 и уровня 2, но я изо всех сил пытался включить переменную из обоих наборов данных в линейная модель смешанного эффекта.
Вот пример: dt1 - это набор данных для учащихся мужского и женского пола в школах a, b и c (в моих данных они расположены по порядку, например a, a, a, б, б, б ...). Результат y - это окончательный результат теста, который является непрерывной переменной.
> dt1
# A tibble: 9 x 3
School gender y
<chr> <chr> <dbl>
1 a m 1
2 b F 3
3 c m 5
4 a F 4
5 b m 2
6 c F 1
7 a m 4
8 b F 3
9 c m 1
В dt2, q и w - переменные на уровне школы
> dt2
# A tibble: 3 x 3
School q w
<chr> <dbl> <dbl>
1 a 2 8
2 b 4 2.5
3 c 4 5
Мне удалось запустить MLM из dt1 следующим образом:
lmer(y~gender +(1|School), data= dt1)
Но как я могу включить переменную из dt2 в предыдущую модель?
Я пробовал это, но это не сработало:
lmer(y~ gender+dt2$q[dt2$q==q] +(1|School), data = dt1)
Подскажите, пожалуйста?