I have a question concerning the conditional probabilities in LCA. Usually, you assign a person to the latent class with the highest probability and then go on with this categorial variable. However, is it not allowed to go on with the probabilities themselves, meaning that you get and calculate with a score for each person and each class? I know that one requirement of LCA is that latent classes are disjunctive and exhaustive. I calculated my analyses with both the categorial variable (= a person belongs to one and only one latent class) and the (intervall scaled) probabilities (= a person has a value for each class). The results seem to me much more plausible using the scores... And in my case it's reasonnable to assume that a person "has" a bit of different classes.
Given that the LCA is estimated with individuals in all classes, using most likely class membership would introduce estimation errors. In Mplus one can add covariates or other relationships to the model and estimate it in one step.
Jon Heron posted on Wednesday, March 07, 2007 - 3:37 am
like you, I've recently had the same thoughts that a model using the posterior probabilities themselves is the way to go (e.g. an mlogit model with pweights in Stata).
Unfortunately, any model which starts off by you importing the data into another package will end up treating your derived class variable/probabilities as manifest whereas they are actually still latent.
Two issues here
1) You're going to underestimate SE's if you assume the class-variable is manifest. I guess it's the same as imputing some missing data and then assuming your dataset is complete.
2) To fit the model properly, the other variables which you bring into your model once you've derived the latent classes should at the very least be permitted to affect the posterior probabilities and with a model including covariates these should be allowed to affect the 'curves' themselves and hence potentially the interpretation of the classes.
This second point is currently causing me serious brain-ache.
Delucchi, Matzger, Weisner Drug and Alcohol Dependence 74 (2004) 235–244 Dependent and problem drinking over 5 years: a latent class growth analysis
For quite a clear stats-methods on the way to proceed
Just a comment on using posterior-probability-weighted multinomial regression for exploring "c ON x" relationships. There is no need for going outside Mplus here - you can do multinomial regression and instead of categorical observed dependent variable categories such as 0, 1, 2, 3 use fractional membership in all dependent variable categories.
Thank you for your answers. I have to explain in more detail what I did. My aim is to predict work outcomes with the clustered variables. That is, I do not predict membership of classes but instead use the classification variable as IV.
I did 2 steps: first I calculated an exploratory latent class analysis and found 4 classes (with PANMARK). In the next step (with SPSS) I used these class-variable as an IV and/or as a moderator in multiple regression analysis to predict work outcomes. Here, I did not use 4 binary IV (for each class) but 4 variables with the response probabilites of the LCA solution. This means that I have a (intervall scaled) number for each person and each class. This is against the idea that a person can only belong to one class, but in my case, I find it very plausible that a person does belong to different classes. Or is latent class analysis simply the wrong method if my assumption is that a person can be member of different classes? Or is it possible to introduce the condition that a person can be member of various classes in LCA? Is my approach described above justifiable?
for a latent class variable c and an observed outcome variable y, you can do this in one step in Mplus. A one-step analysis avoids errors of estimation for the parameters and their SEs. You don't say "y ON c", but the means/intercepts/thresholds will change by default across the c classes. Ex8.6 in the UG shows how to do this for a categorical outcome in the context of a growth mixture model.
If you prefer to think of a person as belonging to several classes, LCA is not the model to use. LCA assumes that a person belongs to only one class. Although the class probabilities that you get at the end gives a person a probability for each class, that is not a feature of the model. It should rather be viewed in line with factor score estimation where each person gets a posterior distribution (given the model and his/her data) and his/her estimated factor score is taken as the peak of that posterior distribution (like most likely class). If you want the model itself to acknowledge partial class membership you are into modeling such as "Grade of Membership". GoM is discussed in the Mplus framework in section 7 of Asparouhov-Muthen (2006) on our web site under Papers, Multilevel Mixture Modeling.
Dear Bengt, thanks a lot for your answer. "Grade of membership" seems to be what I need. Up to now, I used PANMARK and Mplus for LCA. Can I calculate GoM with them or which program instead is required?
Dear Bengt, Can you please tell me what the mplus command is to calculate a grade of membership model? I have read the Asparouhow/Muthen paper (2006) and went through the user guide, but was not able to find the command for this special type of lca. Thanks a lot, Marius
See Section 7 of the following paper which is available on the website:
Asparouhov, T. & Muthen, B. (2008). Multilevel mixture models. In Hancock, G. R., & Samuelsen, K. M. (Eds.), Advances in latent variable mixture models, pp. 27-51. Charlotte, NC: Information Age Publishing, Inc.
Hello Linda, Thanks for the reference. I must admit I'm struggling to convert the general GoM model specification from the paper into something concrete in the Mplus language. Is it something along the lines of example 10.7 in the user's guide? Can you also tell me how one should properly incorporate the dummy variables for items (Xqij)? Any advice appreciated, Nick
Sorry, not that I am aware of. Erashova has done Bayes work on it, but not using Mplus. Anybody else knows?
Chris Harper posted on Wednesday, January 14, 2015 - 10:38 am
Can a GoM model be estimated in Mplus with ML? I'm a little confused after reading through this thread. According to Asparouhov and Muthen (2006) GoM models can be estimated in Mplus, but this thread seems to indicate otherwise. If it is possible to estimate a GoM model in Mplus, how should the syntax differ from a LCA?