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

Question

我想使用update()函数更新模型的随机部分，特别是添加一个随机效应。大多数例子（help("update")，help("update.formula")，lme4:用 R 进行混合效应建模）都关注于模型的固定部分。在下面的示例中，如何使用update()从fm0到fm1？

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

- Richard

Answer 1

我怀疑这对你没有用处，但是你必须去除随机效应，然后再加上需要的效应：

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

- Remko Duursma

---

这是完全正确的。很遗憾没有更紧凑的方法来做到这一点，但很难看出如何指定“在此指定的随机效应项内更新”（特别是因为内部机制只是将更新规范传递给R的内置update.formula()方法）。 - Ben Bolker