r - NLME对象的提取范围

标签 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/

相关文章:

r - 如何查找 R 的 Google Analytics 包的客户端 ID 和 key

r - 像通用线一样向现有的 ggplot 添加线

r - lme 模型的用户定义函数适合 : error

r - 从 nlme 对象中提取随机公式

r - 如何将受密码保护的 Excel 文件读入 r 中?

python - 我应该如何为 R 和 Octave 方法(可能使用 Python)进行快速 GUI 开发?

python - 通过 RPy 查看 R 数据示例(示例 : lmeSplines)

r - 如何使用lme4将没有随机效应的模型与具有随机效应的模型进行比较?

r - 随机效应变量是否自动被视为 R 中 lmer(或 lme)的因子?

r - 如果环境在创建对象的时间和计算对象的时间之间发生变化,R 中的延迟计算会导致意外结果