Latent class analysis with nominal in...
Message/Author
 Christopher Hautmann posted on Wednesday, December 03, 2008 - 12:45 pm
Hello,

I want to compute a latent class analysis with nominal indicators (example 7.7 in the user's guide). I expected that for the values of the indicators conditional probabilities would be provided. Instead in the result section means are reported. How can a nominal variable have a mean? What mistake do I make? Does Mplus not make the same analysis like Latent Gold for example?

Regards

Christopher Hautmann
 Bengt O. Muthen posted on Thursday, December 04, 2008 - 9:25 am
The estimated item parameters that are reported are logits for the different categories of each nominal indicator. These logits can be translated into probabilities very easily using the multinomial logistic regression expression - this is discussed in Chapter 13 of the User's Guide.
 Ali posted on Tuesday, October 04, 2016 - 8:17 am
I run LCA with 4 nominal indicators and each indicator has 3 categories.
But, I am no sure if I get the conditional probabilities for each indicator correctly.
From the model results, there were means for each category in each latent class.
For example, in the latent class 1, means for ST53Q01#1 is 1.080 ,and ST53Q01#2 is 1.753, but ST53Q01#3 was fixed as reference group .
So, log odds (ST53Q01#1|C=1)=1.08, log odds (ST53Q01#2|C=1)=1.753,
and log odds (ST53Q01#3|C=1)=0 in Class 1. And,each log odds is exponentiated and summed. To get the conditional probability for each category in is that each expoentiated value is divide by the sum.

Is it the correct procedure?
 Bengt O. Muthen posted on Tuesday, October 04, 2016 - 2:31 pm
Check how it is done in the UG chapter 14. See the multinomial logistic regression example with covariates all = 0.
 Sam Crawley posted on Monday, October 22, 2018 - 5:46 pm
I have calculated the estimates for nominal indicators using the method described in chapter 14. However, is there any way to also get standard errors, as is provided for categorical indicators in probability scale?

Thanks.
 Bengt O. Muthen posted on Monday, October 22, 2018 - 5:56 pm
If you express them in Model Constraint, you also get their SEs.
 Sam Crawley posted on Monday, October 22, 2018 - 7:12 pm
Thanks. Is there an example of the syntax to do this somewhere? In particular, I'm not sure how to reference the estimates for the means for each combination of class/nominal category in the MODEL CONSTRAINT section.
 Bengt O. Muthen posted on Tuesday, October 23, 2018 - 11:14 am
The V8 UG pages 55-557 show an example with 4 nominal categories and in the first step only intercepts are used which corresponds to your situation. These are the intercepts you have for one of your observed nominal variables in one class. A nominal DV's intercepts are referred to as y#1, y#2, etc. Label these in the Model command and generalize from there in Model Constraint.
 Sam Crawley posted on Tuesday, October 23, 2018 - 3:16 pm
Thanks very much.

Just in case anyone in the future needs a more complete example (assuming 3 latent classes, and a nominal indicator variable y with 4 categories):

MODEL:
%c#1%
[y#1-y#3] (c1y1-c1y3);

%c#2%
[y#1-y#3] (c2y1-c2y3);

%c#3%
[y#1-y#3] (c3y1-c3y3);

MODEL CONSTRAINT:
NEW(c1yprob1-c1yprob4 c2yprob1-c2yprob4 c3yprob1-c3yprob4 c1y4 c2y4 c3y4);

DO(1,4) c#y4 = 0;
DO(\$,1,4) DO (%,1,3) c%yprob\$ = exp(c%y\$) /
(exp(c%y1) + exp(c%y2) + exp(c%y3) + exp(c%y4));
 Bengt O. Muthen posted on Tuesday, October 23, 2018 - 5:31 pm
That's a FAQ.