替换R中的lmer系数

Question

在下面的岗位上，

Replace lm coefficients in [r]

我也有兴趣改变混合模型的系数与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