R NLME对象的提取范围
我有一个具有相关性的NLME对象。我想知道怎样才能提取出 模型的相关性。我正在运行模拟,所以我不能只是阅读摘要,然后手动获取 所以我的模型看起来像这样:R NLME对象的提取范围,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的代码才能找
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)
除非你也有数据,否则我无法运行该模型来查看弹出的内容。你能使用内置的数据集吗?您能显示摘要以指示您要提取的部分吗?我现在已将其更改为标准数据集。@shujaa比我快,但您现在有了答案。举个简单的例子,事情就简单多了。@Flick先生希望我不是在偷猎。。。虽然我很确定一周前我的名声比你高——现在你比我高出了50%:)别担心@shujaa。我只是喜欢问题得到回答。它不必在我身边。我承认最近我有点沉迷于在这个网站上回答问题。希望我的兴趣很快会降到一个更合理的水平。我也在做同样的事情。我只是想看看
print.summary.corStruct
中的类是如何改变coef
的行为的。泛型函数的使用非常有趣。是的,我仍然对类更改的必要性感到困惑。在lme4
中迁移到S4课程,无疑获得了赞赏。您可以将其缩短一点——您不需要summary()
coef(mod$modelStruct$corStruct,unconstraint=FALSE)
@BenBolker谢谢!采纳了你的建议。也感谢您维护lme4,回答这个问题让我更加感激。
> store_range
range
58.82908