#### VISUALIZE CORRELATION MATRIX ###
eta=matrix(rnorm(200*5),ncol=5)
lam=matrix(0,nrow=100,ncol=5)
for (i in 1:5) lam[(20*i-19):(20*i),i]=rnorm(20,0.7,0.3)
eps=matrix(rnorm(200*100),ncol=100)
Y=eta%*%t(lam)+eps
loadings=loadings(factanal(covmat=cor(Y),factors=5))[1:100,1:5]
qgraph.loadings(loadings)
Run the code above in your browser using DataLab