替换 R 中的 lmer 系数

【问题标题】：Replace lmer coefficients in R替换 R 中的 lmer 系数
【发布时间】：2020-02-07 14:32:48
【问题描述】：

从下面的帖子中，

我还对更改装有 lmer 的混合模型的系数感兴趣。例如在以下形式的模型中：

mod <- lmer(y ~ x1 + x2 + x3 + (1|class/subjects), data=data)

如何将 x1 和 x2 的系数更改为另一个数字（0 或 0.1 或 1 等）。

【问题讨论】：

为什么？只需传递零作为预测变量。
我想看看这是否可能。我以0为例

标签： r lme4 mixed-models

【解决方案1】：

我没有看到这样做的正当理由，您可能会遇到一些强大的龙，但很容易做到。

library(lme4)
fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)
summary(fm1)$coef
#             Estimate Std. Error   t value
#(Intercept) 251.40510   6.823773 36.842535
#Days         10.46729   1.545958  6.770744

fm1@beta[names(fixef(fm1)) == "Days"] <- 0
summary(fm1)$coef
#            Estimate Std. Error  t value
#(Intercept) 251.4051   6.823773 36.84253
#Days          0.0000   1.545958  0.00000

【讨论】：