いつの間にか(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")
で,答えが違う。