#fa <- factanal(x,4,covmat=Harman74.cor$cov)
#pc <- principal(Harman74.cor$cov,4)
#pcv <- varimax(pc$loading)
#round(factor.congruence(fa,pcv),2)
#
#round(factor.congruence(pcv,fa),2)
# Factor1 Factor2 Factor3 Factor4
#PC1 1.00 0.60 0.45 0.55
#PC2 0.44 0.49 1.00 0.56
#PC3 0.54 0.99 0.44 0.55
#PC4 0.47 0.52 0.48 0.99
#compare with
#round(cor(fa$loading,pcv$loading),2)
Run the code above in your browser using DataLab