我有一个具有相关性的NLME对象。我想知道如何提取范围 模型的相关性。我正在运行模拟,因此我无法阅读摘要,并手动获取。
所以我的模型看起来像这样:
library(MASS)
library(nlme)
lme(fixed=temp ~ time, random=~1|day,correlation=corExp(form=~time),data=beav1)
所以我想在这里得到相关参数。
答案 0 :(得分:2)
花了很多钱去寻找!我必须查看nlme:::print.summary.lme
的代码才能找到$modelStruct$corStruct
,然后在nlme:::print.summary.corStruct
找到那里。
这应该有效
library(datasets)
library(nlme)
mod <- lme(fixed = temp ~ time, random = ~1|day,
correlation = corExp(form=~time), data=beaver1)
store_range <- coef(mod$modelStruct$corStruct, unconstrained = F)
屈服
> store_range
range
58.82908