lmer系数是否标准化？

Question

我有一个非常基本的问题；也许有点太基础了，无法找到有用的应对措施。

我正在使用使用以下代码的lmer函数来计算多级模型：

lmer(H1_rirs, data= df_long_cl, REML = T)

我是否正确假设所检索的系数未标准化？如果可以，是否有一种简便的方法可以将它们标准化？

最好

卡罗琳

Answer 1

是的，默认情况下，不应用标准化。 如果您想获得标准化的系数，一种方法是在拟合模型之前对数据进行标准化。 effectsize-package中有一个可靠的实现。或者，您可以进行一些事后标准化（也effectsize-package）。后者产生不同的结果，最准确的是在模型拟合之前将数据标准化。

Answer 2

除了@daniel 对 effectsize 包的建议之外，还有“伪标准化”系数（Hoffman，2015 年），其中响应和预测变量根据预测水平进行标准化：

library(lme4)

m <- lmer(Reaction ~ Days + (Days | Subject), data = sleepstudy)

effectsize::standardize_parameters(m, method = "pseudo")
#> # Standardization method: pseudo
#> 
#> Parameter   | Coefficient (std.) |       95% CI
#> -----------------------------------------------
#> (Intercept) |               0.00 | [0.00, 0.00]
#> Days        |               0.68 | [0.48, 0.88]

^{由 reprex package (v2.0.0) 于 2021 年 6 月 7 日创建}