test <- congeneric.sim(short=FALSE)
round(cor(test$observed),2) # show a congeneric correlation matrix
f1=factor.pa(test$observed,1,scores=TRUE)
round(cor(f1$scores,test$latent),2) #factor score estimates are correlated with but not equal to the factor scores
Run the code above in your browser using DataLab