0

私はカラカルの素晴らしい例から二分データの因子分析を行っていますが、psychパッケージのfa.poly関数によって生成されたオブジェクトから因子を抽出するのに苦労しています。

誰かが私がfa.polyオブジェクトから要因を抽出するのを手伝ってくれる（そして相関関係を見る）ことができますか？

動作例については、カラカルの例を参照してください。

4

1 に答える 1

2

この例では、次のオブジェクトを作成します。

faPCdirect <- fa.poly(XdiNum, nfactors=2, rotate="varimax")    # polychoric FA

したがって、faPCdirectのどこかに、必要なものがあります。str()の構造を調べるために使用することをお勧めしますfaPCdirect

> str(faPCdirect)

List of 5
 $ fa   :List of 34
  ..$ residual    : num [1:6, 1:6] 4.79e-01 7.78e-02 -2.97e-0...
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:6] "X1" "X2" "X3" "X4" ...
  .. .. ..$ : chr [1:6] "X1" "X2" "X3" "X4" ...
  ..$ dof         : num 4
  ..$ fit      
      ...skip stuff....
  ..$ BIC         : num 4.11
  ..$ r.scores    : num [1:2, 1:2] 1 0.0508 0.0508 1
  .. ..- attr(*, "dimnames")=List of 2
  .. .. ..$ : chr [1:2] "MR2" "MR1"
  .. .. ..$ : chr [1:2] "MR2" "MR1"
  ..$ R2          : Named num [1:2] 0.709 0.989
  .. ..- attr(*, "names")= chr [1:2] "MR2" "MR1"
  ..$ valid       : num [1:2] 0.819 0.987
  ..$ score.cor   : num [1:2, 1:2] 1 0.212 0.212 1

つまり、このオブジェクトは5つのリストであり、最初の要素が呼び出され、2x2行列であるfaという要素が含まれています。score.corあなたが欲しいのはオフダイアゴナルだと思います。

> faPCdirect$fa$score.cor
          [,1]      [,2]
[1,] 1.0000000 0.2117457
[2,] 0.2117457 1.0000000
于 2012-07-19T02:06:01.580 に答える