Message/Author 


Hi Folks, This is my first post, and so please excuse any normative slipups on my part. We would like to estimate a SecondOrder Latent Class Model, with a mix of categorical/nominal observed indicators and categorical/nominal latent variables on the first and the second order. Graphically, the model looks exactly like the model in example 5.6 of the Mplus manual, but with all categorical/nominal variables instead of the continuous observed and latent variables given in the example. Our question is, what's the most efficient Mplus code to achieve this type of latent class model? Any and all help is greatly appreciated. Thanks! 


Secondorder LCA has the flavor of UG ex 8.14 but without c2 on c1 and without the observed indicator (threshold) restrictions. For the second level, you would have the Overall model with c1 on c and c2 on c. For the first level you would have a Model c1 and a Model c2. Look at Tech1 to see that you get what you want. 


Hi Folks, Is this example (8.14) the place I should start if I want to run a multiple indicator GMM, similar to example 6.14? I've got longitudinal data that doesn't quite provide the complete picture if its reduced to one factor, so I'm wondering if I can run (or if I should even consider running) something akin to a multiple indicator GMM. Also, my data is cluster randomized, so I'm going to be analyzing it in a multilevel context after I've completed the first step of the single level analysis. Thanks. 


You should start with the one class model in Example 6.14. Then estimate it for two classes etc. by adding TYPE=MIXTURE; to the ANALYSIS command and the CLASSES option to the VARIABLE command. 


Execellent. I have another question regarding LCA and a different data set. I have data of 7 binary drug and sexual risk behaviors at baseline, immediate postintervention, 6 months, and 12 months. The LCA at baseline found a 3class solution. I would like to look at the three subpopulations found at baseline and see how they responded to the intervention. I guess one approach would be to simply run growth models in each subpopulation on each risk behavior one after the other. Or I could do the multiple indicator growth model as a multiple group growth analysis. Is there another approach you can recommend where I could do all of this in one step? 


If your 7 binary indicators are well summarized by 3 latent classes, perhaps you would be interested in seeing how subjects transition between those classes over time and how those transitions are affected by the intervention. If so, Latent Transition Analysis is suitable. We can send more information about Mplus input for intervention analysis if this is where you want to go. You mention growth modeling in each of the 3 classes, which sounds like you think there is variation across people within classes. If that's the case, perhaps LCA is not the best measurement model, but instead Factor Mixture Analysis (not imposing measurement invariance across classes). As an advanced approach, an FMA measurement model can be used for a combined growth and transition analysis. As yet another alternative, you can see the 7 binary items as indicators of a (say) single continuous factor, where this factor has 3 classes and you use that factor as your dependent variable in a GMM  that is then a multipleindicator GMM. Its measurement model can be seen a Mixture Factor Analysis (classinvariant measurement parameters). 

Mplususer posted on Wednesday, January 30, 2008  4:56 pm



Hi there, Following your response to Matthew Cole on July 11, I have a question about multipleindicator GMM. Suppose we have scores for 7 items across 3 time points and we don't know if they are all loaded on the same factor across time. Do we need to run factor analyses at each time point to preliminarily examine if the items are loaded on this factor, and then use the factor scores as the values of a dependent variable for GMM analysis? Is Example 6.14 still the place to start the analysis? Thanks. 


You need to check for measurement invariance across time. Example 6.14 without the growth model is the place to start. I would then continue with a multiple indicator growth model. I would not use factor scores with less than 15 factor indicators. 

Mplususer posted on Thursday, January 31, 2008  11:50 am



Thanks, Linda! 


With respect to the secondorder LCA described in Bengt's 10/1/06 response (above), can this be extended to the case where c1 and c2 are both GMMs? 


I think you mean that c1 and c2 are latent class variables influencing growth factor means for GMMs. Yes, that kind of secondorder LCA is possible to do in Mplus. 


No, actually the other way around, where the second order latent class variable (i.e., the c in example 8.14) is a function of two latent categorical variables, c1 (externalizing trajectory class) and c2 (internalizing trajectory class). If this is possible, I seem to be having a problem figuring out how to specify it in the model statements. Specifically, I keep getting error messages indicating that all growth factors need to be specified in the overall statement. However, in the model I'm describing, there are growth factors specific to the c1 and c2 models, but not second order latent class variable. I'd appreciate any help you can offer. Thanks. 


Seems like you can give the two growth models in the Overall part and then in the cspecific parts mention the classspecific growth factor means for the respective growth model. So like: Model c1: %c1#1% [i1]; %c1#2% [i1]; and the same for Model c2 using [i2] (and other growth factor means. This implies that [i1] is held equal across the c2 (and c) classes, and so is only a function of c1 classes. And analogously for [i2]. Check in Tech1 that you get the equalities you need. That's quite a complex model. 


Thanks so much! The model runs beautifully. 


Hi All  Two general questions about LCA with a secondorder factor. I am conducting an analyses similar to example 7.18, except with 4 and 5 categorical indicators for c1 and c2, respectively. How do you determine the number of classes for c1 and c2? Through previous LCA, or through this secondorder analyses? Secondly, can you use the auxiliary(e) function for a second order LCA in the same way you do for standard LCA? If not, how would compare expected means across secondorder classes? Thanks again! 


Q1. I would do it in separate regular LCAs for each c. Q2. We don't recommend aux(e) any longer; see the tables at the end of the paper on our website: Asparouhov, T. & Muthén, B. (2014). Auxiliary variables in mixture modeling: Using the BCH method in Mplus to estimate a distal outcome model and an arbitrary second model. Web note 21. Try aux(bch) and see what you get. 


Thank you very much for the suggestion and reference. When trying to run with BCH, I get the error: *** ERROR in VARIABLE command Auxiliary variables with E, R, R3STEP, DU3STEP, DE3STEP, DCATEGORICAL, DCONTINUOUS, or BCH are not available with TYPE=MIXTURE with more than one categorical latent variable. Is there an error in my coding, or if not, any suggestions to get around this restriction? Thanks again, Ryan 


This cannot be done with more than one categorical latent variable. 


Thank you for clarifying and the continued help. If I have a model based off of two categorical latent variables, is there any another other way to estimate and test differences in auxiliary variables across groups other than the auxiliary functions? 


Not that I know of. 

Back to top 