提取NLME对象的范围

时间:2014-05-30 21:23:44

标签: r nlme

我有一个具有相关性的NLME对象。我想知道如何提取范围 模型的相关性。我正在运行模拟,因此我无法阅读摘要,并手动获取。

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

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

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

1 个答案:

答案 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