Extracting EFA factor loadings and factor correlations

[Tracey Chau]

Dear all,

I am analysing data where the participants are clustered by family. It is easy enough to run an EFA on Mplus with TYPE = COMPLEX to account for this, but I have struggled to find an equivalent on R. The analyses that I'm running on R require multiple EFAs (i.e., 1 factor, 2 factors, ... all the way up to n factors) as well as correlations between all of the factors from all n EFAs. I need the factor loadings from each of the n EFAs too. 

The syntax for all these analyses have already been written on R but accounting for the clustered nature of the data has been the biggest challenge as there's nothing in the psych or lavaan packages as far as I can see. 

I am wondering whether I can utilise the MplusAutomation package to run these EFAs on Mplus with TYPE = COMPLEX and somehow extract the equivalent of the psych fa $loadings and $weights outputs in order to run the rest of the syntax?

Thank you so much in advance for your help, any direction or advice is greatly appreciated.

Kind regards,

Tracey