当使用promax旋转在R中使用stats包中的factanal()函数时,将得到因子相关性。
tmp = cbind(rnorm(200,2,1.5),rnorm(200,2,1.5),rnorm(200,2,1.5),rnorm(200,2,1.5),rnorm(200,2,1.5),rnorm(200,2,1.5))
print(factanal(tmp, 3, rotation="promax"))
但是,我似乎找不到如何计算这些相关性(R文档中没有帮助)。我也无法使用系数得分的cor()自己重新创建它们。
我猜测可能是多因素相关,但我不确定。
答案 0 :(得分:0)
使用以下代码计算得出,可以在https://github.com/SurajGupta/r-source/blob/master/src/library/stats/R/factanal.R中找到:
tmat <- solve(tmp$rotmat)
R <- tmat %*% t(tmat)