裏 RjpWiki

Julia ときどき R, Python によるコンピュータプログラム,コンピュータ・サイエンス,統計学

因子間相関係数行列がへん?

2011年06月10日 | ブログラミング

いつの間にか(R 2.13.0 ですでに),factanal(   rotation="promax")  で,因子間相関係数行列が出力されるようになっているが,因子との対応がおかしいのではないか?
Loadings は SS loadings の大きい順に並べ替えられているのに,因子間相関係数行列の行と列はその並べ替えができていないのじゃないか?
つまり,
temp <- factanal(x, factors=n, rotation="none")
ans <- promax(loadings(temp))
solve(t(ans$rotmat) %*% ans$rotmat)

のようにしたときと,
factanal(x, factors=n, rotation="promax")
で,答えが違う。

コメント (1)
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする

PVアクセスランキング にほんブログ村

PVアクセスランキング にほんブログ村