如何在lmer（）中使用update（）进行随机部分？

Question

我想使用update()函数来更新模型的随机部分，特别是添加随机效果。大多数示例（help("update")，help("update.formula")，lme4:mixed effects modeling with R）都侧重于模型的固定部分。在下面的示例中，如何使用fm0从fm1转到update()？

library(lme4)
(fm0 <- lmer(Reaction ~ Days + (1 | Subject), sleepstudy))
(fm1 <- lmer(Reaction ~ Days + (1 + Days | Subject), sleepstudy))

Answer 1

我怀疑这对你的情况有用，但你必须删除随机效果，然后在背面添加所需的：

update(fm0, . ~ . -(1|Subject) + (1 + Days | Subject))