lmer 系数是否标准化？答案

【问题标题】：Are lmer-coefficients standardized or not?lmer 系数是否标准化？
【发布时间】：2020-03-29 00:41:20
【问题描述】：

我有一个非常基本的问题；可能有点太基本了，无法找到有用的回复。

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

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

我是否正确假设检索到的系数是非标准化的？如果是，有没有简单的方法来标准化它们？

最好的，

卡罗琳

【问题讨论】：

H1_rirs 是如何定义的？
嗨 Roland，它的定义是这样的：H1_rirs<-commresc ~ 1 + cvar1c + (1 + cvar1c | ID_new) 两个变量都是以人为中心的，尽管我想我不应该以我的结果为中心（它在更大的模型中充当中介)

标签： r lme4 coefficients standardized multilevel-analysis

【解决方案1】：

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

【讨论】：

【解决方案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 创建于 2021-06-07 (v2.0.0)}

【讨论】：