#### 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
# Run qgraph:
qgraph.efa(cor(Y),5)
# Show crossloadings:
qgraph.efa(cor(Y),5,crossloadings=T,cut=0)
Run the code above in your browser using DataLab