Two-part growth model with continuous... PreviousNext
Mplus Discussion > Growth Modeling of Longitudinal Data >
 V X posted on Tuesday, February 06, 2007 - 12:01 pm
Hi, Linda,

I used to use Mplus v3 for the two-part growth model and I am now using v4.

I noticed that new commands have been added for the model (DATA TWOPART) in the latest version of Mplus and the original data are log transformed automatically for the continuous part of the model.

I am wondering, what if I ignore the DATA TWOPART command in Mplus v4 enviroment, will the log tranformation automatically apply to the analysis?

Thank you.
 Linda K. Muthen posted on Tuesday, February 06, 2007 - 12:07 pm
No, if you don't use DATA TWOPART, the log transformation will not be done. But you can do in using DEFINE.
 V X posted on Friday, February 09, 2007 - 5:15 pm
Thank you,Linda. This is of great help.

I have another question with regard to Monte Carlo simulation study for a two-part growth model for a continuous outcome.

I plan to study the impact of alternative distribution assumption to parameter estimation. So, can I use Mplus v4 to generate a longitudinal data set with alternative assumption, say, gamma distribution for the continuous part? If not, would you give me some suggestions on how to do it?

Have a nice weekend !
 Bengt O. Muthen posted on Saturday, February 10, 2007 - 9:11 am
Mplus does not offer the gamma distribution. But using mixtures you can obtain almost any distribution you like. For example, using Mplus to generate data with 2 latent classes, each with a normal distribution, the mixture has a skewed distribution. We used that approach in our 2002 article which you find on our web site:

Muthén, L.K. & Muthén, B.O. (2002). How to use a Monte Carlo study to decide on sample size and determine power. Structural Equation Modeling, 4, 599-620.
 Nicki Bush posted on Monday, March 12, 2007 - 8:06 am
I am examining interactions (between 2 observed variables) within a 2-part model. Are there special considerations on how to probe significant interactions in this context--I mean, given that there are significant interactions predicting both the binary and continuous parts, do I do anything other than calculate log odds for the binary outcome ones and plot 1SD below/above mean on the continuous outcome ones (ala Aiken & West)?

Also, I am unaware of any papers that have examined interactions in a 2-part model. Can you refer me to any writings or examples of this?

Thank you so much!
 Linda K. Muthen posted on Saturday, March 17, 2007 - 9:14 am
There is nothing special about the observed variable interactions in a two-part model. The issues would be the same as for regular regression and the discussion in any regression text would apply.
 Louise Lafortune posted on Thursday, October 04, 2007 - 2:42 pm
I want to run a two-part model with covariates on cross-sectional data (re: analyzing health care costs data with a preponderance on zeros). I already created my binary and continuous outcomes (but I know it is not necessary since V4)

Is e.g 6.16 the best starting point? If yes, how do I specify the model?

Or should I treat my problem as a censored-inflated regression WITH correlation between the two regression (e.g. 3.3)? If so, how can I make sure the censoring limit is zero? For the correlation, do I just add y1 with y1#1 in the model?

 Linda K. Muthen posted on Friday, October 05, 2007 - 9:10 am
The residual covariance between the two parts of a regression model is not identified so it doesn't matter if you run the two parts separately or together. You might was to see work at the Rand Corporation by Naihua Duan on two-part regression modeling.
 Louise Lafortune posted on Friday, October 05, 2007 - 11:15 am
Thank you Linda.
 Steven Haas posted on Tuesday, March 11, 2008 - 4:43 pm
I have a question about the transform option in the data twopart command.

The mplus guide lists a log transformation as the default. It also lists a transform=none as an option.

If I do not include the transform option statement does mplus automatically log the continous part of the model?
 Linda K. Muthen posted on Tuesday, March 11, 2008 - 5:03 pm
 Michael Spaeth posted on Wednesday, February 04, 2009 - 6:43 am
In a two-part growth approach, how would you label the growth curve of the y-part, when "type = missing" was applied? Normally, one labels this curve "frequency of behavior for those engaging in the behavior in question". However, since I'm using type=missing random "nonusers" are not excluded from the estimation of growth means of the y-curve.
If a person was nonuser at t1 and then user at t2, according to FIML (type=missing), a sort of "imputation" is done with the missing y-value at t1 (is that correct), despite the fact that this person is a nonuser!?

Thanks, Michael
 Linda K. Muthen posted on Wednesday, February 04, 2009 - 8:09 am
Please see the following paper for a discussion of these issues:

Olsen, M. K, & Schafer, J., L. (2001). A two-part random effects model for semicontinuous longitudinal data. Journal of the American Statistical Association, 96, 730-745.
 Hanno Petras posted on Sunday, November 14, 2010 - 3:49 pm
Dear Linda and Bengt,

I am trying to setup a twopart growth mixture model with two categorical latent variables, i.e., where c1 depicts developmental heterogeneity in the probability and where c2 depicts this for the conditional mean of y. While the model converges, it seems that my model set up is incorrect because for both c1 and c2 both measurement models (i.e., u and y) are used to define the classes. How can I keep these two part of the model separate and allowing for different number of classes for u and y. Thank you for your feedback. Below is my abreviated syntax:


iu su qu cu | etc.
iy sy qy cy | etc.

c2 on c1;

Model c1:
%c1#1% [ iu su qu cu];
%c1#2% [ iu su qu cu];
%c1#3% [ iu su qu cu];
Model c2:
%c2#1% [ iy sy qy cy];
%c2#2% [ iy sy qy cy];
%c2#3% [ iy sy qy cy];
 Linda K. Muthen posted on Monday, November 15, 2010 - 8:32 am
This looks correct. Try it out and if you have a problem send it to support.
 Lauren Mitchell posted on Monday, February 09, 2015 - 1:50 pm
Hello! Is it possible to run a two-part longitudinal model with individually varying times of observation using Mplus? e.g.,
iu su | b1 b2 b3 b4 AT a1 a2 a3 a4;
iy sy | c1 c2 c3 c4 AT a1 a2 a3 a4;

Thank you!
 Linda K. Muthen posted on Monday, February 09, 2015 - 3:57 pm
Yes, this should be possible.
 Philipp Jugert posted on Wednesday, February 03, 2016 - 12:48 am
Hello! I have done two-part growth modeling for 3 ethnic groups separately. Is it admissible to test for differences in intercepts and slopes between ethnic groups from these separate models using simple t-Tests or ANOVAs? I know it's theoretically possible to use Mixture Modeling with known classes to compare groups. But I don't want to go down that road unless absolutely necessary.
 Bengt O. Muthen posted on Thursday, February 04, 2016 - 7:03 pm
Yes, this is ok because the estimates are uncorrelated.
 Philipp Jugert posted on Friday, February 05, 2016 - 9:57 am
Dear Bengt,
Thanks a lot for your answer. Would you perhaps have a reference to justify this procedure of model comparison?
 Bengt O. Muthen posted on Friday, February 05, 2016 - 5:59 pm
No, but it follows from first principles. You have independent samples in the different groups and you have no across-group parameter equalities which means that the estimates are uncorrelated.

To check that you get the right results you can run the 3 groups in one multi-groups run with no equalities and use MODEL Constraint to express whatever difference you are interested in - it gives you SEs and z-scores.
 Philipp Jugert posted on Tuesday, April 19, 2016 - 3:50 am

We did two-part growth modeling for a continuous outcome (victimization) with a preponderance of zeros (between 30 to 50 % of the sample depending on time point). Now, a reviewer questions the use of this method and suggests to use either bootstrapping or Bayes instead.
However, when we use bootstrapping interesting effects found for the binary part are no longer there. So I would like to keep using two-part. Are there some arguments that could justify the use of two-part over these other suggested methods?
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