When using the factanal() function from the stats package in R using the promax rotation, you are given factor correlations.
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"))
However, I can not seem to find how these correlations are calculated (no help in the R documentation). I also fail to recreate them myself using cor() on the factor scores.
I'm guessing it may be a polychoric correlation, but I can't be certain.
psychpackage to see if you get the same results. $\endgroup$
factor-rotation. You will find the formula how factor correlations are computed. And this answer (towards the end) explains that whether factor scores will or will not reproduce these correlations depends on the method of the scores estimation. $\endgroup$