替換 R 中的 lmer 系數

Question

從下面的帖子，

替換 [r] 中的 lm 系數

我也對更改裝有 lmer 的混合模型的系數感興趣。 例如在以下形式的模型中：

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

如何將 x1 和 x2 的系數更改為另一個數字（0 或 0.1 或 1 等）。

Answer 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