我正在尝试从使用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来解决这个问题,但我也希望将其应用于更复杂的结构(多个随机斜率,嵌套的随机变量等),并且可能包括加法或随机斜率。谢谢!

最佳答案

您可以通过以下方式访问它们

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

关于r - 从nlme对象中提取随机公式,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/26912618/

10-12 17:30
查看更多