Mixture modeling with multilevel data PreviousNext
Mplus Discussion > Latent Variable Mixture Modeling >
 Anonymous posted on Tuesday, February 19, 2002 - 3:06 pm
I would like to know if a multilevel or multisample mixture model is possible with Mplus. If so, can you please point me to further readings?

Many thanks!
 Linda K. Muthen posted on Tuesday, February 19, 2002 - 3:55 pm
Mplus can currently do Mixture Complex with and without missing data. It can also do Mixture mutisample using training data to specify the groups. Mixture multilevel is currently under development.
 Scott Grey posted on Wednesday, June 14, 2006 - 11:47 am

I'm trying to run a LPA using 'TYPE IS MIXTURE COMPLEX' but I do not get standard errors in my output. Here's the code:

FILE IS "C:\Documents and Settings\insthealthsa4\My Documents\DARE\
External prevention programming\CFA3.dat";

NAMES ARE n1q44a_7 n1q44a_8 n1q44a_9 n1q44b_7 n1q44b_8 n1q44b_9
n1q44c_7 n1q44c_8 n1q44c_9 n1q44d_7 n1q44d_8 n1q44d_9 v10q52 v10q55
n2q49 utq39 crsswlk ms7dist treatms inclass outclass dare other
t7 t8 t9 t10 msdist;
USEVARIABLES ARE inclass outclass dare other;
AUXILIARY = crsswlk;
CLUSTER IS ms7dist;
CLASSES = class(5);

STARTS = 50 5;
LOGHIGH = +25;
LOGLOW = -25;
LOGCRITERION = 0.0000001;
CONVERGENCE = 0.000001;
MCONVERGENCE = 0.000001;

inclass-other WITH outclass-other;



 Linda K. Muthen posted on Wednesday, June 14, 2006 - 1:47 pm
I'm afraid this information does not tell me why you don't get standard errors. Please send your input, data, output, and license number to support@statmodel.com.
 Hao Duong posted on Saturday, October 18, 2008 - 11:02 am
Dr. Muthen,
I am confused about interpretations for ib sb on w, and c on w in example 10.9 in Mplus user's guide in 2007.
Would you please explain them for me?
Thank you
I appreciate all your help!
Hao Duong
 Linda K. Muthen posted on Monday, October 20, 2008 - 8:43 am
All regressions are linear regressions. ib and sb are continuous growth factors. c#1 is a random intercept, also a continuous latent variable.
 Miguel Villodas posted on Monday, August 22, 2011 - 10:22 am

I have two questions. I am attempting to establish a GMM with four time points. However, the time points are not equally spaced. I do have an indicator of time (age at each assessment).

1. Would this preclude me from using a traditional GMM? I was under the impression that I would have to use a two level model, rather than LGCM because of the unequal time between interviews, and regress my outcome on my time variable.

2. Also, if multilevel modeling is required, would it still be possible to examine classes of trajectories?
 Bengt O. Muthen posted on Monday, August 22, 2011 - 10:39 am
1. No. A single-level model is sufficient. Just use time scores that reflect the non-equidistance.

2. Yes, you can use between-level latent class variables.
 Miguel Villodas posted on Monday, August 22, 2011 - 11:52 am
Thank you very much Dr. Muthen. The problem that I have run into is that the spacing of each assessment varies individually. Is this still possible?
 Miguel Villodas posted on Monday, August 22, 2011 - 2:24 pm
I should clarify as this is unclear. I have been able to run these models using the TYPE = RANDOM command for the individually varying time scores. I am wondering if there is a way, using this option, to examine classes of trajectories? If so, how would I specify this model in MPlus?
 Bengt O. Muthen posted on Monday, August 22, 2011 - 3:49 pm
Individually-varying times of observations is handled by the TSCORES option (see UG) and it can be combined with Type=Mixtures, although it takes a little longer.
 FANG BEE LING posted on Thursday, November 01, 2012 - 2:33 am

I'm doing multilevel latent class analysis with 12 items which have 4 categorical options. I tried few times but keep pop up some error message. And I'm not sure what's gone wrong. Hope someone can advise me on this. many thanks.

my input instruction:
NAMES ARE wel mor dif enj str qui bor lik hel oth uni job schid ;
USEVARIABLES = wel mor dif enj str qui bor lik hel oth uni job ;
CATEGORICAL ARE wel mor dif enj str qui bor lik hel oth uni job ;
MISSING ARE all (9);
WITHIN=wel mor dif enj str qui bor lik hel oth uni job ;
STARTS = 20 10;
C#1; C#2; C#1 WITH C#2;

Categorical variable WEL contains 158 categories.
This exceeds the maximum allowed of 10.
 Linda K. Muthen posted on Thursday, November 01, 2012 - 7:46 am
It sounds like you are reading your data incorrectly. Perhaps the number of variable names is not the same as the number of columns in your data set or you have blanks in your data. If you can't figure this out, send your files and license number to support@statmodel.com.
 jilke@fsw.eur.nl posted on Thursday, April 18, 2013 - 3:03 am
I am implementing a IRT multilevel mixture model with item bias effects in Mplus. My structure has items nested in individuals, nested in countries.

I would like to code the mixtures nominal, yielding discrete random effects. For identification I use effect-coding for the mixtures. But when using the following syntax, there are no differences for the latent means across mixtures; both are zero.

Am I applying the correct syntax for effect-coding of mixtures?

NAMES = id item1 item2 item3 country;
USEVARIABLES = item1 item2 item3;
CATEGORICAL = item1 item2 item3;
CLUSTER= country;
CLASSES= eta3 (2);
BETWEEN= eta3;


item1* (a)
item2 (b)
item3 (c);

item1* (a)
item2 (b)
item3 (c);

%eta3#1% !Mixture #1
[ETA2b] (d); !Mean of latent var for mixture #1

item2; !item bias effects

%eta3#2% !Mixture #2
[ETA2b] (e); !Mean of latent var for mixture #2

item2; !item bias effects

model constraint:
0= d+e; !effect coding for mixtures
 Linda K. Muthen posted on Thursday, April 18, 2013 - 9:33 am
The latent variable mean must be fixed at zero in one class for model identification. So d and e in MODEL CONSTRAINT are not both identified.
 Artur Pokropek posted on Thursday, February 13, 2014 - 8:51 am
I'm running Multilevel Latent Class model similar to model presented in Henry&Muthén(2010). The core of the model is multilevel logistic model: resp on negd1 negd2 negd3; and 2 latent classes are specified. Estimation of the model works when I've two regression in each of the latent class:

resp on negd1 negd2 negd3;
cw#1 on time ;
resp on negd1 negd2 negd3 ;
resp on negd1 negd2 negd3 ;
resp ;
[resp$1] ;
resp ;
[resp$1] ;
resp ;
[resp$1] ;

But the problem appears when I want to have empty logistic model in class 2. i.e model without explanatory variables. I was trying several specifications. None of them worked. While I'm declaring ON statement in %OVERALL% part I get ON statements in all classes. When ON statement is not declared in %OVERALL% part I'm not allowed to specify it for class 1. Is there a way to specify logistic model with explanatory variables in class 1 and empty model in class 2?

Thank You!
 Bengt O. Muthen posted on Friday, February 14, 2014 - 12:02 pm
Please send to Support the data and output for the case where you specified the ON statement in the Overall part and fix the slopes in class 2.
 Lisa M. Yarnell posted on Wednesday, April 09, 2014 - 6:25 pm
Hello, I plan to use UG Example 10.12 (two-level LTA with a covariate) for my analyses. I have students nested in schools. I understand the code in the example, except I want to clarify one aspect of it.

In the code below, why are the indicators for the latent classes modeled at the between level? If these reflect individual responses (such as from individual students), wouldn't those be on the within level?

Or is it that because we are estimating probabilities (or mean responses) for the items for persons, conditional on class, this becomes an average across persons--no longer on the within level? Thank you.

c2 ON c1 x;
c1 ON x;
c1#1 ON w;
c2#1 ON c1#1 w;
c1#1 c2#1;

[u11$1-u14$1] (1-4);
[u11$1-u14$1] (5-8);
[u21$1-u24$1] (1-4);
[u21$1-u24$1] (5-8);
 Linda K. Muthen posted on Thursday, April 10, 2014 - 9:45 am
In multilevel modeling, all mean/threshold/intercept parameters are on the highest level.
 Lisa M. Yarnell posted on Thursday, April 10, 2014 - 5:18 pm
OK, may I request your input on the following questions regarding this section of the code:
c1#1 ON w;
c2#1 ON c1#1 w;
c1#1 c2#1;

Understanding the code above:
1) Why is c2#2 ON c1#2 (and other combinations such as c2#2 ON c1#1) not above, similar to the second line? I believe that the code above regresses the cluster-level (average or intercept) latent status for class 1 of c2 on that for class 1 of c1. This is part of the random intercept setup. But wouldn't regressing the second class of each latent class variable make sense to do as well?

2) By the same token, why does the above code not show "c1#2 c2#2;" as well? Is it that by allowing intercepts for the first class for each latent class variable to vary across clusters, these are already free to differ from those for the second latent classes?

Altering the code for my data:
3) If I have 3 latent class variables rather than 2, and want to have random intercepts, I would also model c3#1 ON c2#1, right?

4) Finally, if I try a model without random intercepts, would I remove both the "c2#1 ON c1#1", and "c1#1 c2#1;" sections of code?

Thank you sincerely.
 Bengt O. Muthen posted on Friday, April 11, 2014 - 11:36 am
There are only 2 classes in which case there is no between counterpart for the second class just like in multinomial logistic regression. With 3 classes you have between counterparts for the first 2.
 Jiwon Shannon Choi posted on Thursday, May 01, 2014 - 5:02 am
Dear Muthen
I am writing this post to ask MLCA analysis.
is it possible to run "multilevel LCA" with covariates and distal outcome simultaneously?
the number of individual cases are 580 nested with 30 organizations.

Thank you in advance.
 Linda K. Muthen posted on Thursday, May 01, 2014 - 10:34 am
Thirty organizations is the minimum you should have. Yes, this model is possible.
 Yasumasa Otsuka posted on Monday, October 06, 2014 - 6:44 pm
I'm now considering to calculate the MOR proposed by Larsen & Merlo (2005). But due to my poor statistical ability, I cannot understand the mathematical expression shown in page 83, 1st line.
Could anyone give me some example with actual number to calculate MOR in this article?
 Bengt O. Muthen posted on Tuesday, October 07, 2014 - 9:23 am
Try contacting the authors, or post on Multilevelnet.
  Chris Kenaszchuk posted on Friday, June 03, 2016 - 2:33 pm
Some segments of Example 10.1, 'Two-level mixture regression for a continuous dependent variable,' are below.

I'd like to modify the program for two objectives:
(1) Incorporate the measurement error table output from step 1/step 2 of manual 3-step estimation, i.e., "Logits for the Classification Probabilities for the Most Likely..."
(2) Temporarily omit the regressions from the overall model in the between part of the model.

The example does not mention %C#2% in the within part of the model. Is it necessary to omit %C#2%? If so, then how could the parameter N#1@ be used for %C#2%?

WITHIN = X1 X2 N; ! N is my addition
! BETWEEN = W; ! Not using for now


Y ON X1 X2;
C ON X1;
[N#1@1.901]; ! As per Webnote 15, appendix E, step 3 of
! manual 3-step estimation
Y ON X2;

! Y ON W; ! No between-level regression
! C#1 ON W; ! No between-level regression
! C#1*1; ! A starting value is not needed
 Bengt O. Muthen posted on Sunday, June 05, 2016 - 11:54 am
It is not necessary to omit %C#2%. We did because we didn't need to say anything about c#2.

So go ahead and try it.
Back to top
Add Your Message Here
Username: Posting Information:
This is a private posting area. Only registered users and moderators may post messages here.
Options: Enable HTML code in message
Automatically activate URLs in message