I have run LCGA and identified 4 class trajectories. How can I find out if their latent class means differ from each other (class 1 from class 2, class 2 from class 3 etc.). Is there some kind of a post-hoc test available and how to request it? Thank you!
Michael Lap posted on Monday, March 12, 2018 - 11:33 am
To clarify- in a 3-trajectory class solution, in the output, there is a section: Parameterization using Reference Class 1
You are probably referring to the means of the intercept and slope growth factors [i] and [s]. If so, simply label them for each class in the Model command and do test of their differences in Model Test.
As for c#2 and c#3, they have to do with class sizes and not trajectory differences.
Thank you, Bengt. 1) Where can I find an example of the Model test? 2) "As for c#2 and c#3, they have to do with class sizes and not trajectory differences." Does it mean that in my case, the size of Class 3 (number of cases) is not different from the size of Class 1? But Class 2 and Class 1 are different in terms of their sizes? Thank you, I appreciate your help.
I would not go about the analysis that way. I would first focus on whether C or C-1 classes is required by looking at BIC. If say 2 classes has a much better BIC than 1 class, I would report the 2-class solution. I would not bother with testing if [i1], [s1] are different from [i2], [s2] - referring the the means of the intercept and slope growth factors in the case of no covariates. But if you want to, simply say
With covariates, these parameters would refer to intercepts in the regressions on the covariates and would therefore be less informative about trajectory differences.
Michael Lap posted on Wednesday, March 14, 2018 - 12:34 pm
Thank you very much, Bengt. Actually, the number of class trajectories was finalized based on the BIC, entropy, and Lo-Mendel-Rubin statistic. By using the Model test, I was hoping to make sure that all detected class trajectories do indeed differ from each other (intercepts and slopes). However, as you said, the Model test is less informative in the situation with covariates in the model.
Dear Drs. Muthen, I am running a 3-class GMM CI with predictors and covariates. I am using the 3-step approach with adjustment for classification errors.
Referring to Micheal's post (the 4th from the top) and your related reply "Yes (assuming you have no C on X).", I was wondering if you could help me interpreting the size of classes when covariates are included in the model.
In my model, when Class 1 is the reference one, and when age, education, sex, smoking habit and alcohol habit are included, I obtain these values:
If you want to estimate the class counts in the situation with covariate it would be best to use the output
FINAL CLASS COUNTS AND PROPORTIONS FOR THE LATENT CLASSES BASED ON THE ESTIMATED MODEL
Unfortunately these don't come with confidence interval yet.
Looking at the intercepts alone is not very useful because in evaluating class counts you would need to include the entire multinomial regression with the covariates. In addition the posterior class probabilities vary across observations because the covariates vary.
For a particular set of covariates you can obtain CI for the posterior class probablities using the multinomial regression formula in model constraints. For the entire population, as an approximation, one can use the average covariate value (see output:tech7) instead of actual covariate values. This can yield good confidence interval estimates for the class counts.