提取 gls 的估计值

Question

我正在尝试从 ML 拟合的模型中提取参数，该模型具有从 corPagel 获得的相关结构（在包 ape 中定义）。感兴趣的参数位于摘要中名为 apVar 的元素内。

library(ape); library(nlme)
cor_lambdaML <- corPagel( 1, avestree, fixed=FALSE )
m <- gls(AUTUMN ~ GREGARIOUSNESS + SEASON + TARSUS + MATURATION + SPRING + MIGRATION + DICHROMATISM + HABITAT + POSTJUVENILE + CONSPICUOUSNESS, correlation = cor_lambdaML, data = avesdata, method = 'ML' )
names(m)
m$apVar

此代码产生以下输出（请下载模型为 m here 的 dput）：

           corStruct     lSigma
corStruct 0.03848807 0.03090663
lSigma    0.03090663 0.04403168
attr(,"Pars")
  corStruct      lSigma 
0.858126838 0.008131981 
attr(,"natural")
[1] TRUE

提取第一个矩阵的值很简单，但是如何提取元素 attr(, "Pars") 中的 corStruct 的值？

Answer 1

不确定您需要挖掘多远，但m$Pars$corStruct 是存储数据的位置。也就是说：m 是一个有很多组件的 list，其中一个是 Pars 。这本身就是一个list，带有一个组件corStruct。

?glsObject 将引导您了解所有列表元素的详细信息。