Message/Author 


I am modelling pubertal growth using LTA. Two time points with a 4 category latent variable at each time point. The LTA is really nice. But I want to model the effect of the transition from T1 to T2 on a distal outcome. How might this be accomplished? Do I need a second order latent categorical variable, as illustrated in the moverstayer model? With thanks. 


I think you are correct that you want a model like Example 8.14 with c pointing to the distal outcome 


I have a model like 8.14 running. But if "c" (the second order latent categorical variable) is categorical with, say, two classes, and the distal outcome is a continuous variable (such as a latent factor), won't there be a problem with the regression estimates? By that I mean "c" is not dummy (0,1) or effects (1,1) coded, but is coded as 1,2. How might that be handled? 


The regression of a continuous distal outcome on the secondorder catgorical latent variable is reflected by the change of means of the continuous distal outcome across the classes of the secondorder categorical latent variable. 

Kaigang Li posted on Thursday, May 29, 2008  12:12 am



Hello Linda, If I have a THREE category categorical distal outcome, U, and three covariates, ie. x1, x2, and x3. How can I write in syntax for looking at the effect of covariates on the U? Is the following commands correct? %OVERALL% U on x1 x2 x3; %c#1% [U$1 U$2]; %c#2% [U1$1 U2$2]; Thanks, Kaigang 


It looks correct. 

Kaigang Li posted on Thursday, May 29, 2008  11:59 pm



Thanks. I have another question. I am reading the chapter "Secondgeneration structural equation modeling with a combination of categorical and continuous latent variables: New opportunities for latent class/latent growth modeling." but I found the examples penn18.inp at http://statmodel.com/examples/penn.shtml are not consistent with the ones used in the book. Could you please instruct me how I can use those examples? Thanks, Kaigang 

Kaigang Li posted on Friday, May 30, 2008  12:22 am



Linda, One more question. In the article "Jung, T. & Wickrama, K.A.S. (2008). An introduction to latent class growth analysis and growth mixture modeling. Social and Personality Psychology Compass, 2, 302317." posted at http://statmodel.com/papers.shtml, the author Specified a singleclass latent growth curve model using the following Model: i s  t1@0 t2@1 t3@2; But in the penn1.inp, the command % OVERALL% was used. I run the program with and without % OVERALL% , the results are different. I am not sure which is correct if I specify a single class model. Thanks for clarification. Kaigang 


If you do not obtain the same results for a singlelevel growth model and a oneclass growth model, you should send your two output files and your license number to support@statmodel.com. The results should be identical. 

Kaigang Li posted on Friday, May 30, 2008  8:24 am



Thanks Linda, I will double check the results. Do you have any comments on the question posted on Thursday, May 29, 2008  11:59 pm right above the question you answered? Thanks, Kaigang 


Colleagues and I have run a LTA with gender as a known class. The next step in our analysis is to include distal outcomes. We have looked at K. Nylund’s dissertation as a guide of how to include distal outcomes in a latent transition model. Her example is helpful in regards to regressing the distal outcomes on a second order mover stayer variable and on estimating the means for a given class within a wave (e.g., class 3 wave 3). However, we are interested in estimating differences in distal outcomes for different transition configurations, separately for boys and girls. For example, we would like to know if girls that transition from class 1 to class 3 have different means on a distal outcome than girls that transition from class 2 to class 3 or girls that remain stationary across the two time points. Is there a way to estimate this in MPLUS? One important thing to note is that we have previously tried to run a mover stayer model with gender as a known class and that model did not replicate (with 2000 starts). Additionally, we tried to run a mover stayer model with gender as a covariate and that model did not replicate either. Any suggestions on how to examine distal outcomes for different transition configurations would be appreciated. Thank you. 


I think you would have to use the dot specification for two latent class variables. See UG. For instance: %c1#1.c2#1% [distal] (p11); [u1$1u5$1] (15); [u6$1u10$1] (15); %c1#1.c2#2% [distal] (p12); [u1$1u5$1] (15); [u6$1u10$1] (610); %c1#2.c2#1% [distal] (p21); [u1$1u5$1] (610); [u6$1u10$1] (15); %c1#2.c2#2% [distal] (p22); [u1$1u5$1] (610); [u6$1u10$1] (610); Model Constraint: new(diff); diff = p12p21; 


I setup code based upon your suggestion. However, when I run this I get the following warning: "*** ERROR in MODEL command Unknown class model name ASB1.ASB2 specified in Cspecific MODEL command." Am I missing a step before setting up the Model asb1.asb2 statement? A portion of the code is below to show how I begin: MODEL: %overall% asb1#1 asb1#2 asb1#3 on H1GI20 race_cat poverty_sum_r; asb2#1 asb2#2 asb2#3 on depress alcw2 anydrugw2; asb2#1asb2#3 on asb1#1asb1#3; Model asb1.asb2: %asb1#1.asb2#1% ... Then sets up the item response parameter constraints and the distal outcome parameters (as well as the model constraint statement). 


Remove Model asb1.asb2: 


I did as you indicated, but now get the following error: *** ERROR The following MODEL statements are ignored: * Statements in Class %ASB1#1.ASB2#1% of MODEL: [ ARRW3_DIC ] !distal outcome variable This error is repeated for each of the asb1.asb2 combinations (4 levels of each) 


Please send the full output and your license number to support@statmodel.com. 


Related to the above described LTA models  I have simplified the model and have everything working in terms of assessing distal outcomes by class consistent with K. Nylund's dissertation approach. To compare distal outcomes across classes I am using the model test approach. However, I cannot seem to locate the actual statistical tests in the output. Can you clarify where these are located (happy to send the output, if that would help). 


It is with the other fit statistics before the parameter estimates. If you can't see it, please send your full output and license number to support@statmodel.com. 


I am using the manual 3step approach to estimate a multiplegroup LTA with a continuous distal outcome (two classes at two time points assuming measurement invariance, so four latent statuses). In the third step, I successfully obtain distal means for each of the four latent statuses separately for each of the two observed groups. Of the eight means estimated, six are as expected, but two are quite different from the others (two of the four within a single group). Is there a way to manually estimate these latent status means from the data file created in step 2? I'm trying to get a handle on why the means for those two statuses are so different from the rest. I appreciate any guidance you can provide. 


Appendix I of our 3step paper shows that the "c2.dat" data set contains n1 and n2, which are the most likely class variables for the two latent class variables. You can get the distal outcome mean for each crossclassification of n1 and n2. 


Thanks for your quick reply. The means estimated in step 3 do no match the means I get from averaging the distal outcomes within each status. For example, for group 2, the step 3 output displays means of 8.83, 43.55, 40.82, 3.21 for statuses 14, respectively. Manually estimated means (exporting status and summarizing in Stata) are 6.51, 8.52, 6.25, and 5.47. The 43.55 and 40.82 are the "quite different" means in question. Thoughts? 


You may want to check your steps in detail against the Nylund et al (2014) 3step LTA article on our website. Also check that the class formations stay the same. If that doesn't help, please send data and outputs from your steps to Support so we can see what's happening. 


To follow up on this thread, the user did not use the proposed steps: To study the influence of a covariate on the latent class variables in an LTA with measurement invariance, one should go by the Appendix KN setups, which takes the measurement invariance approach. See the Appendices of the paper on our website: Asparouhov, T. & Muthén, B. (2014). Auxiliary variables in mixture modeling: Threestep approaches using Mplus. Structural Equation Modeling: A Multidisciplinary Journal, 21:3, 329341. The posted version corrects several typos in the published version. An earlier version of this paper was posted as web note 15. Appendices with Mplus scripts are available here. 


Hello, I have a question regarding looking at different "directions" of effects in one model. I have an LTA with four parenting profiles at four time points, and I'm modeling children's emotion regulation as both a covariate and a distal outcome. I'm trying to get as close as I can to a something that conceptually resembles an autoregressive crosslagged panel model. I understand it's not wholly possible, as I cannot ask for x/y on c. I ran one model ("parentdriven") where I used emotion regulation as a covariate of profile membership and as an influence on transition probabilities. I ran another model ("childdriven") where I allowed the means of emotion regulation as an outcome to vary across the profiles and used Wald's tests to examine differences. My question: Can I combine these two sets of syntax? Would I gain different information than by running them separately? I am struggling to wrap my mind around it. If you know of any examples that examine bidirectional/transactional relationships with latent classes, I would be very appreciative of your recommendations. Thank you very much for your time. 


It sounds like you have an LTA with a 4class latent variable c_t determined by parent outcomes and you want to combine that with a child outcome at the different time points, say z_t. And you want z_t to both influence c_t to c_t+1 transitions and be influenced by c_t (or c_t1 perhaps). If that is a correct understanding, I think it is doable. Even though you don't say z_t ON c_t1 the z_t means can change over the c_t1 classes and therefore represent the hypothesis. But it is a complex model; I haven't tried it. I would recommend starting with only 2 time points. 


hello, I'm trying to write a moverstayer model with covariates of the ms class. the classes at t1 and t2 are supposed to be sequential steps, that is people in c1.1 can only transition upwards to c2.2 at t2 or remain in c2.1; but people in c1.2 cannot go in c2.1. I was able to do this with Parameterization = Probability; however to insert covariates I can only use Parameterization = logit; and then I am not sure how I can constrain the last class to not to go downwards thank you for your reply 


You can accomplish the same using Param = logit. See the V7Part2 handout on the Mplus short course video page under the heading: Dutch Mplus Users Group and Mplus Version 7 workshop, Utrecht, August 2012 Videos and handouts from 3day Version 7 workshop. See in particular slides 57 and 6971. 


thank you Bengt, it's quite clear now. I have another question on the UG 8.15 it's specified %OVERALL% c1 ON c; with c = moverstayer while in the slide you suggested it is c1 c2 c3 ON c !Relating c1, c2, c3 to c: (Movers) Why is that different? and what is the best way to specify the moverstayer model? 


I have also another question: I would need to differences in some distal variables between mover and stayer for only certain classes:e.g., for people who remained in class 1 (c#1.c#2) and people in class1 that moved to the next class (c#1.c#2) does it make sense to compare the coefficients like this: %c#1.c#1% [distal] (d1); %c#1.c#2% [distal] (d2); model test: d1 = d2; and above all, are those coefficients the logit of what exactly? the probability of the distal variable to be high in each pair of classes? thanks again 


Q1. That's a difference between the logit and the probability parameterization. Q2. You can do it that way. The logit for a binary distal is related to the probability of 1 instead of 0 for the distal: P(distal=1) = 1/(1+exp(L)) where L is the logit. 


me again with another question: I'm requesting the means of a distal variables for each transition pattern with: %overall% [distal]; However I noticed that the output show means even for patterns for which the most likely size is zero. Does it make sense to keep those in the model (I guess they shouldn't be interpreted though) or it would be more correct to constrain those to zero? 


I assume you don't hold the distal means equal across patterns but they are free. If so, small class/pattern frequencies (which would give zero most likely class freq) give unstable estimates  if the pattern freq is zero it is not even identified. I think you can fix them to anything. 


Hi, I am trying to run an LTA model (two time points) with covariates and distal outcomes. I was able to use example 8.14 to model the effect of covariates on the transition from c1 to c2. However, I am finding it tricky to model the effect of the c1c2 transition on a distal outcome, modeled as a continuous latent variable, f2. I am also controlling for the outcome at time 1, f1, by regressing it on f2. To regress the c1c2 transition on f2, I would need to vary the means of the continuous latent variable across the classes of model c1.c2. However, in order to do this I would need an extra categorical latent variable in my model. Since I would like to regress each individual transition on the outcome, is there a way to perform this regression without having to use a MoverStayer model? Thanks for your help. 


I think you mean that you want to regress f2 on c1 x c2 classes and not the other way around. I don't see why you need an extra latent class variable to let the f2 mean vary  why not just use the "dot" approach for the c1 x c2 combination classes: %c1#1.c2#1% [f2]; etc 


Dr. Muthen, I read in the manual that you need more than two categorical latent variables to specify MODEL c1.c2. However, I realized that I could use the %c1#1.c2#1% approach under the %OVERALL% label without any problems. Also, I do want to regress f2 on the c1c2 transition, not the other way. Sorry for the confusion and thanks for helping me resolve my issue. Thanks, Raghav 


You got it. 


In modeling an LTA with a distal outcome that varies only over the last process (i.e., C3 in my case), how does one obtain the means for the C3 classes only? In the Nylund diss p. 127 she reports these means yet when I estimate them in model the output gives the means of the distal for every latent class pattern (rather than only the last process)... What am I missing? With thanks, J 


If you have specified the model to allow the distal to vary across only the C3 classes, then the distal means should vary over only those classes. So if C3 has say 4 classes, there should be only 4 distinct distal means. 


Many thanks. I see now my two means for the relevant C3 process are the only two means provided for all LT patterns as should be. 


Hello again: I am estimating a 4time point LTA with a mover stayer variable. I am also estimating a distal at T5 for the stayers. All of the above is in line with the examples provided in the Nylund et al. series of documents. When I estimate the distal for the stayers, I also get different means for the mover patterns. This is problematic because for some patterns it seems to cause the nonpositive issue or needs fixed to avoid singularity (I'm assuming because that pattern is too sparse to generate an estimate?). Is there any way to fix means for the movers at 0 on the distal? I am not interested in them for the moment. 


That should be possible  you just have to state for which class combination (perhaps using the dot command) you want which mean. I don't know that you want them to be zero, but rather work with equality constraints. 


This worked nicely, thank you. I set the means for the movers equivalent. 


Hello, I am estimating a 2time point LTA with predictor of class membership at both points in time (using the manual 3step procedure). I am now assessing distal outcomes of class membership at T2 (while controlling for predictors). I have added [outcome] to each T2 class. However, the class proportions and latent transition probabilities changed quite a bit when I ran the LTA with the distal outcomes. I tried constraining some of the proportions/probabilities using MODEL CONSTRAINT and that seemed to work. Is this an appropriate route to take? Next I used MODEL TEST to compare outcome means across T2 classes. I also requested cintervals. I am puzzled about the results. the cintervals for the outcome means across T2 classes overlap heavily even when comparing class outcome mean parings that, according to the corresponding Wald test are statistically significantly different (most cases with p < .001). Would you have any insight as to why this might be? I have tried using MODEL CONSTRAINT to create diff variables, but I get a number of error messages (these dissapear if I remove the probability/proportion constraints, but then the outcome means across classes are quite different than those in the model with MODEL TEST). Is there a way for me to calculate the diff variables, their standard errors and cintervals manually? Thank you for your help, Dayuma 


I think it is more expedient if you send the output you use for what you describe in your second paragraph and we take it from there. Send to Support along with your license number and this message. 

Daniel Lee posted on Thursday, February 21, 2019  2:34 pm



Hi, I am trying to look at whether certain transitions (e.g., low to high, moderate to low) report different levels of an outcome (e.g., depression). Would you have documentation and syntax that demonstrates how to implement this type of analysis in full (how to set up the analysis and model command)? For example, I saw some helpful syntax above, but wasn't sure if I should still input C2 on C1 under %OVERALL% or input the original commands from the LTAs under Model C1 and Model C2 (e.g., Model C1: %C1#1%...%C1#2%...MODEL C2: %C2#1%). 


Just compare the depression means/thresholds for the combinations of classes that you are interested in  using Model Constraint. 

Daniel Lee posted on Saturday, February 23, 2019  9:22 am



Thank you, but I am particularly interesting in setting up the LTA correctly using the dot method (e.g., %c1#1.c2#1%) you mention in this discussion...I have not seen the dot method in the UG manual...my question is how I would set up the dot method if I conducted my LTA using the method below: MODEL: %OVERALL% c2 ON c1; MODEL c1: %c1#1% c2 ON x; [u11$1] (1); [u12$1] (2); [u13$1] (3); %c1#2% c2 ON x; [u11$1] (6); [u12$1] (7); [u13$1] (8); MODEL c2: %c2#1% [u21$1] (1); [u22$1] (2); [u23$1] (3); %c2#2% [u21$1] (6); [u22$1] (7); [u23$1] (8); 


Say %c1#1.c2#1% c2 on x (10); [u11$1] (1); [u12$1] (2); [u13$1] (3); [u21$1] (1); [u22$1] (2); [u23$1] (3); %c1#1.c2#2% c2 on x (10); [u11$1] (1); [u12$1] (2); [u13$1] (3); [u21$1] (6); [u22$1] (7); [u23$1] (8); %c1#2.c2#1% c2 on x (11); [u11$1] (6); [u12$1] (7); [u13$1] (8); [u21$1] (1); [u22$1] (2); [u23$1] (3); %c1#2.c2#2% c2 on x (11); [u11$1] (6); [u12$1] (7); [u13$1] (8); [u21$1] (6); [u22$1] (7); [u23$1] (8); 

Takeo Kato posted on Wednesday, May 08, 2019  6:20 am



I am estimating a 2time point LTA with three steps. I have been done LTA with covariates to examine the latent transition probabilities allowing for covariates. Now I'm in the process of LTA with the distal outcomes which are continuous variables. I could get the mean of each distal outcome variable on latent class assignment at c2. However, I couldn't get the mean of each distal outcome variable on each transition pattern. How to solve it? My syntax is below. MODEL: %OVERALL% c2 ON c1; c1 ON male ga moed01 mobm01 faag02 ptin100 pwbwratio10 anbmi_whosds18; MODEL c1: %c1#1% [c18#1@13.690]; [c18#2@11.434]; [c18#3@10.155]; [c18#4@0.0001]; c2 ON moed01 mobm01 anbmi_whosds18; %c1#2% . . MODEL c2: %c2#1% [c32#1@6.726]; [c32#2@4.722]; [c32#3@0.367]; [c32#4@1.461]; [vlstd40 vlcv40 vlmv40 vlsv40 vldv40]; %c2#2% [c32#1@4.451]; [c32#2@7.814]; [c32#3@6.713]; [c32#4@1.461]; [vlstd40 vlcv40 vlmv40 vlsv40 vldv40]; . . 


You can use the dot approach instead of Model c1, Model c2. So %c1#1.c2#1% etc for all combinations. 

Takeo Kato posted on Wednesday, May 08, 2019  8:19 pm



Thank you. I tried it based on your instruction. However, I got the error messages. MODEL: %OVERALL% c2 ON c1; c1 ON male ga moed01 mobm01 faag02 ptin100 pwbwratio10 anbmi_whosds18; %c1#1.c2#1% [c18#1@13.690]; [c18#2@11.434]; [c18#3@10.155]; [c18#4@0.0001]; c2 ON moed01 mobm01 anbmi_whosds18; %c1#1.c2#1% [c18#1@2.902]; [c18#2@1.908]; [c18#3@3.978]; [c18#4@5.281]; c2 ON moed01 mobm01 anbmi_whosds18; . . %c1#1.c2#1% [c32#1@6.726]; [c32#2@4.722]; [c32#3@0.367]; [c32#4@1.461]; [vlstd40 vlcv40 vlmv40 vlsv40 vldv40]; %c1#1.c2#1% [c32#1@4.451]; [c32#2@7.814]; [c32#3@6.713]; [c32#4@1.461]; [vlstd40 vlcv40 vlmv40 vlsv40 vldv40]; . . 

Takeo Kato posted on Wednesday, May 08, 2019  8:21 pm



sorry, separated post The error message here *** ERROR The following MODEL statements are ignored: * Statements in Class %C1#1.C2#1% of MODEL: C2#1 ON MOED01 C2#1 ON MOBM01 C2#1 ON ANBMI_WHOSDS18 C2#1 ON MOED01 C2#1 ON MOBM01 C2#1 ON ANBMI_WHOSDS18 C2#1 ON MOED01 C2#1 ON MOBM01 C2#1 ON ANBMI_WHOSDS18 C2#1 ON MOED01 C2#1 ON MOBM01 C2#1 ON ANBMI_WHOSDS18 C2#1 ON MOED01 C2#1 ON MOBM01 C2#1 ON ANBMI_WHOSDS18 C2#2 ON MOED01 C2#2 ON MOBM01 C2#2 ON ANBMI_WHOSDS18 . . *** ERROR One or more MODEL statements were ignored. Note that ON statements must appear in the OVERALL class before they can be modified in classspecific models. Some statements are only supported by ALGORITHM=INTEGRATION. 


The error message is clear, saying "Note that ON statements must appear in the OVERALL class before they can be modified in classspecific models." You can add such a statement in Overall with slopes fixed @0. We ask that postings be limited to one window. For longer messages, send to Support. 

Back to top 