我试图将线性混合模型拟合到具有以下特殊特征的数据集:
有一个分类变量 A 取值 {1,2,...,K},它被认为是 {1,2} 的固定效应和 {3,...,K} 的随机效应。我应该如何在 R 中的lmer函数中实现它?
更具体地说,考虑最简单的模型y_i = b_0 + A + error。
(一个接近的解决方案是分别为 {1,2} 和 {3,...,K} 估计两个模型,但这会影响对误差项方差的估计。)
我试图将线性混合模型拟合到具有以下特殊特征的数据集:
有一个分类变量 A 取值 {1,2,...,K},它被认为是 {1,2} 的固定效应和 {3,...,K} 的随机效应。我应该如何在 R 中的lmer函数中实现它?
更具体地说,考虑最简单的模型y_i = b_0 + A + error。
(一个接近的解决方案是分别为 {1,2} 和 {3,...,K} 估计两个模型,但这会影响对误差项方差的估计。)