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.
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