r - 从 nlme 对象中提取随机公式

Question

我试图从使用 lme 构建的模型中提取随机结构，但除了固定公式之外，我似乎什么也得不到。例如。，

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

这对于 lmer 来说是可能的使用 findbars() :

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

我希望能够提取:

# ~ age | Subject
# (Days | Subject)

我可以使用 regexpr 来解决这个问题但我也希望这适用于更复杂的结构(多个随机斜率、嵌套随机变量等)，这可能包括加性或随机斜率。谢谢!

Answer 1

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

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