从nlme对象中提取随机公式

Question

5

我正在尝试从使用lme构建的模型中提取随机结构，但似乎除了固定公式外无法得到其他内容。例如：

library(nlme)
fm1 <- lme(distance ~ age, Orthodont, random = ~ age | Subject)
deparse(terms(fm1))
# "distance ~ age"

使用findbars()，lmer可以实现这一点：

library(lmerTest)
fm2 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)
findbars(formula(fm2))
# [[1]]
# Days | Subject

我希望能够提取以下内容：

我想获得以下信息：

# ~ age | Subject
# (Days | Subject)

我可以使用regexpr来实现这个，但我也希望它适用于更复杂的结构（多个随机斜率、嵌套随机变量等），可能包括加性或随机斜率。谢谢！

- jslefche

1个回答

网页内容由stack overflow 提供, 点击上面的

可以查看英文原文，
原文链接

- konvas · Accepted Answer

您可以通过以下方式访问这些内容：

fm1$call$fixed
# distance ~ age
fm1$call$random
# ~age | Subject