r - 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