我使用 lme4 拟合了一个线性混合效应模型,但是,它是数据的层次结构,我有两个来自级别 1 和级别 2 的不同数据集,但我正在努力将两个数据集中的变量包含在线性混合效果中模型。
这是一个示例:dt1 是 a、b 和 c 学校中男女学生的数据集(它们在我的数据中按顺序排列,如 a、a、a、b、b、b ...)。结果 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
我能够从 dt1 运行 MLM,如下所示:
lmer(y~gender +(1|School), data= dt1)
但是如何在以前的模型中包含来自 dt2 的变量?
我已经尝试过了,但没有奏效:
lmer(y~ gender+dt2$q[dt2$q==q] +(1|School), data = dt1)
请问有什么建议吗?