如何在lmer（）中使用update（）作为随机部分？

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

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

2017-07-19 Richard

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

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

2017-07-19 10:41:53

这是完全正确的。这太糟糕了，没有一种更简洁的方法来做到这一点，但很难看出如何指定“在*这个指定的随机效应项中更新*”（尤其是当内部机器只是将更新规范传递给R的内置' update.formula（）'方法） –

回答