我正在使用R库lme4的glmer函数来拟合使用拉普拉斯近似的一般线性混合(GLM)模型和二进制响应变量BINARY_r。我有一个二级固定效果变量('FIXED',比方说)和两个二级交叉分类随机效应变量(比如'RANDOM1'和'RANDOM2')。一级二进制响应变量BINARY_r分别嵌套在两个二级变量的每一个中。 logit函数用作链接函数,用于表示响应变量的非高斯特性。两个随机效应变量之间的交互效果由'RANDOM1:RANDOM2'表示。所有三个独立变量都是分类的。该模型采用
形式BINARY_r~FIXED + RANDOM1 + RANDOM2 + RANDOM1:RANDOM2。
'FIXED'和'BINARY_r'缺少数据,我希望通过对这两个变量中的每一个应用多次插补来探索模型的改进。
然而,我很不清楚如何使用MI在R中使用glmer生成一个新函数,该函数与原始函数相同,但现在包括FIXED和BINARY_r的插补数据。你能帮忙吗?
非常感谢提前