使用 GLMER 在混合效应模型中编码组级变量的正确语法

机器算法验证 r 混合模式 lme4-nlme 重复测量 多层次分析
2022-03-26 11:26:26

我在个人级别的数据集上运行多级对数赔率回归。每个人都属于一个地区,并且既有个人层面的变量(因人而异),也有地区层面的变量(每个地区都是常数)。让我们调用我们的解释变量 ind_var1、ind_var2、dis_var1、dis_var2、我们的区域 dis 和我们的输出变量 out

我一直在寻找有关如何编写代码的正确语法,这就是我所取得的成就:

mod1 = glmer(out ~ ind_var1 + ind_var2 + dis_var1 + dis_var2 + (1|dis), data = data, family = binomial(link = 'logit'))

我的问题是,上面的代码是构建模型的正确语法,还是需要编辑区级解释变量(dis_var1 和 dis_var2)的语法?

1个回答

当某人更多地来自多级背景和/或使用过诸如 HLM 或 MLwiN 之类的软件时,通常会混淆“变量变化的水平”与混合效应模型,其中(据我记得)有必要指定变量变化的水平。混合效应模型并非如此。软件“知道”变量变化的水平,所需要做的就是将变量包含在公式的固定效应部分中。

您使用的模型公式看起来非常好。