0

我使用 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)

请问有什么建议吗?

4

0 回答 0