R:fit混合效应模型

假设我们有以下线性混合效应模型：

我们如何在

R

中拟合此嵌套模型

目前，我尝试了两件事：

将模型重写为：

然后使用

lme4

包中的

lmer

函数拟合混合效应模型，并将

Xi

作为随机和固定效应协变量，如下所示：

lmer（y~X-1+（0+X |主题））

但是，当我将结果传递给BIC进行模型选择时，它总是选择最简单的模型，这是不正确的

我首先尝试在

X\u I

上回归

y\u I

，并将

X\u I

视为固定效应，然后我将得到斜率的估计值，即

phi\u I

向量。然后将

phi_i

作为新的观察值，并再次对

C_i

进行回归，得到

beta

。但这似乎不正确，因为我们不知道实际问题中的

C_i

，而且似乎

C_i

和

beta

共同决定了系数

那么，在

R

中是否还有其他方法适合这种模型？我的错误在哪里

---

谢谢你的帮助

看起来更像是CrossValidated.com类的问题看起来更像是CrossValidated.com类的问题