NLME对象的提取范围

Question

我有一个有关联的NLME对象。我想知道如何从模型中提取关联的范围。我正在运行模拟，所以我不能只是阅读摘要，并得到它的手动。

我的模型看起来是这样的：

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

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

Answer 1

花了很长时间才找到的！我必须看一下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