如何在 lmer() 中对随机部分使用 update()？

Question

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

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))

如何在 lmer() 中对随机部分使用 update()？

how to use update() for random part in lmer()?

r

lme4