r - NLME 对象的提取范围

Question

我有一个具有相关性的 NLME 对象。我想知道如何从模型中提取相关范围。我正在运行模拟，所以我不能只阅读摘要，然后手动获取它。

所以我的模型看起来像这样：

library(MASS)
library(nlme)
lme(fixed=temp ~ time, random=~1|day,correlation=corExp(form=~time),data=beav1)

所以我想在这里得到相关性的参数。

Answer 1

这需要大量的挖掘才能找到！我必须查看代码 for nlme:::print.summary.lmeto find $modelStruct$corStruct，然后 atnlme:::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