Message/Author 


Is it possible to do an EFA with a threelevel model? If yes, is there an example dataset and MPlus code anywhere? I couldn't find an example. 


No this is not yet available. 


Since there is no threelevel EFA option, can you recommend an analytical approach that I could use with MPlus 7? I looked at the ICC's for level 2 and level 3 in my data set, and they are usually between .15  .21 for all the outcomes. So, it indicates I should really have a threelevel EFA. So, would you recommend doing a twolevel EFA (ignoring either the second level or the third level)? Then, should I do a threelevel CFA based on the results for the twolevel EFA? Alternatively, I could move directly to a threelevel CFA and try several different models. Thanks for any help! 


I would do the EFA using Type=Complex Twolevel; and use level 2 for the twolevel clustering and level 3 for complex. Note also that you want a large number of units on both the 2nd and 3rd levels (at least 20, preferably 30100). 


So, I have 14 consecutive days of data nested within 98 individuals nested within 49 dyads. Based on the manual, I should use a twolevel model. By using Type = Complex Twolevel, am I then accounting for both the individual and dyad level clustering in my EFA? If not, what is that command doing? My ultimate goal is to find the factor structure for several dailylevel variables that takes into account the nested nature of the data. So, can the output of a multilevel EFA potentially produce a different factor structure at each level (i.e. time, individual, and dyad)? In addition, I read thru the user's guide and examples for CH 4 & 9, but I am still not sure if I am coding this multilevel EFA correctly. First, am I coding the three levels correctly (or two levels in your program)? Second, did I choose the correct EFA options to account for these three levels? See below for code. TITLE: EFA with continuous and count factor indicators DATA: FILE = daily_data.dat; VARIABLE: NAMES = y1y6 individual dyad MISSING = BLANK; COUNT = y2 y4; WITHIN = y1y6; BETWEEN = individual; CLUSTER = dyad; ANALYSIS: TYPE = COMPLEX TWOLEVEL EFA 1 3 UW 1 3 UB; ESTIMATION = FIML; 


I would treat this as a multivariate rather than a multilevel model. This means the data would be in wide format rather than long format. The nonindependence of observations is then handled by multivariate analysis. The modeling is more flexible allowing more parameters for each dyad member. 


Can you point me to an example of a multivariate model with threelevel data? Thanks so much! 


See Examples 9.21, 9.22, and 9.23. 


I have a followup question about this. We are evaluating a measure of therapist feelings after psychotherapy sessions, and the level 1 data consists of repeated measurements over sessions with some patients having only a few sessions while some have over one hundred. Level 2 consists of 200 patients, and level 3 is about 100 therapists. I thought the idea of using a two level EFA while controlling for therapist differences using TYPE=COMPLEX TWOLEVEL that Bengt mentioned would be perfect since we will use the measure on a patient level in later analyses (and thus we are not interested in the factor solution for the therapist level). However, I wonder if there is a problem that there are so many therapists who see only one patient (the mean number of patients per therapist is about 1.8)? Also, I wonder how to interpret the level 1 factor solution? Level 2 is reasonably clear, it should be the way different therapist feelings cluster together across patients, after adjusting for therapist differences. But I'm not sure about how to interpret the factor solution for the repeated measurements? Does it have something to do with factor stability over time, after adjusting for differences between patients? Any thoughts on this would be great! Best wishes, Fredrik Falkenström 


Regarding paragraph 1, I don't see that this would cause a problem. Regarding paragraph 2, your interpretation seems reasonable. 


Great, thanks! Fredrik 


A quick followup question: Is it possible to estimate an EFA with Type = Complex Twolevel using categorical dependent variables (i.e. ordinal level with a scale from 04)? It doesn't seem to work when I try WLSMV? Fredrik 


Sorry, the scale is ordinal with range 03 (not 04 as I wrote). Fredrik 


TYPE=EFA COMPLEX TWOLEVEL is not available. See pages601602 of the user's guide where the available combinations of types are shown. 


Ok, thanks. Fredrik 


I'm thinking of alternatives where I can use WLSMV. I thought I'd use only one session (thereby skipping the repeated measures level) and estimate a two level model with the therapist level (level two) estimated as an unstructured covariance matrix (since I don't think the data is enough to estimate a separate factor model for level 2, since many therapists only see one patient). But then I realized that model fit criteria probably look too good when I have an unstructured level 2 model? Especially the CFI, but I guess the Chisquare test will also look too good? I also thought of using a onelevel factor model with level two as Complex. Perhaps this is best, but since I chose only one session I got a warning that standard errors may be biased if I used the Useobservations command. When I instead used the Subpopulation command (as recommended in this warning), I got an error that in one of the categories the response was zero for all observations (and the model was not estimated). Do you have any recommendations in this situation? Thank you very much for your help! Best, Fredrik 


Thinking a bit further, I may be able to specify the above into two questions: 1. Am I right that when estimating a saturated betweenlevel model in addition to the withinpart of the model (the one I am interested in), the inclusion of a part of the model with perfect fit between modelimplied and observed covariance matrices will bias the Chisquare test of overall model fit toward nonsignificance compared to if the within part of the model had been somehow tested "on its own"? If not, this may be our preferred approach. 2. Does it matter that the standard errors may be biased when estimating a model using Type = Complex for a subpopulation using the "Useobservations" command, if we are not really interested in the significance of the parameters but in their values (i.e. size of factor loadings) only? Best, Fredrik 


I can't comment on all the issues since some touch on more general analysis strategies better suited for a general discussion list like Multilevelnet, but here are a couple of quick Mplusrelated comments. Note that 2level analysis can be done using WLSMV. Don't worry about using the Subpopulation command  it would seldom change results in a noticeable way. It is fine to use a saturated betweenlevel model but it also uses more parameters than are typically needed  often a 1factor model is sufficient on between. Please try to keep your postings within one window. 


Ok, thanks! Sorry for the long message, I was thinking along the way. Fredrik 


Thank you in advance for the help. I'm trying to do a threelevel EFA but the scale only has 3 items. I want to assess the psychometric utility of the three items loading on to one factor. My results indicate I have 0 degrees of freedom, since I only have three items. I'm trying to figure out if I can equate or impose other constraints on my model to make it interpretible. The research question I'm trying to address is that I have data consisting of supervisor ratings on a 3item measure for 56 clinicians nested within 17 supervision groups. The number of ratings per clinician is variable: 2  16. So, time nested within clinicians and clinicians nested within supervision group. Here is what I was trying: TITLE: "3Level, 3Item Multilevel EFA" ; DATA: FILE IS "EFA01.txt" ; VARIABLE: NAMES ARE study PtID SsnNum DaysBL Clncn ConGroup ConTmPrd DaysCon Enthu Skills Partic ; USEVARIABLES ARE Enthu Skills Partic ; CLUSTER = ConGroup clncn ; MISSING ARE ALL (99999); MODEL: %WITHIN% W_ESP by enthu@1 skills partic ; %BETWEEN% B_ESP by enthu@1 skills partic ; ANALYSIS: TYPE = TWOLEVEL COMPLEX EFA 1 1 UW 1 1 UB ; OUTPUT: TECH1 SAMPSTAT ; 


With 3 items you can fit only 1 factor so EFA is not relevant. 

Back to top 