r - 提取 NLME 对象的范围

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

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

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

所以我想在这里获取相关性的参数。

最佳答案

费了好大劲才找到!我必须查看 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

关于r - 提取 NLME 对象的范围，我们在Stack Overflow上找到一个类似的问题：https://stackoverflow.com/questions/23963618/