0

corPagel (パッケージ ape で定義) から取得した相関構造を持つ ML によって適合するモデルからパラメーターを抽出しようとしています。対象のパラメーターは、概要の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 の 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 の値を抽出するにはどうすればよいですか?

4

1 に答える 1

0

どこまで掘り下げる必要があるのか​​ わかりませんm$Pars$corStructが、データが保存されている場所です。つまり m、listには多くのコンポーネントがあり、そのうちの 1 つが ですPars。それ自体はlist、コンポーネントを持つ ですcorStruct。

?glsObjectすべてのリスト要素の詳細への道を歩み始めます。

于 2014-02-03T15:17:40.743 に答える