Plot with individually-varying time PreviousNext
Mplus Discussion > Growth Modeling of Longitudinal Data >
 Shige Song posted on Friday, January 26, 2007 - 2:51 pm
I am trying to get an estimated growth trajectory out of a model with individually-varying time of observations but had some problems. To simplify the explanation, I will borrow an example from the Mplus manual.

Here is a modified version of Example 6.12 in the Mplus manual:
USEVARIABLE = y1-y4 a11-a14;
TSCORES = a11-a14;
MODEL: i s | y1-y4 AT a11-a14;
SERIES = y1-y4 (s);

According to the manual (Pp. 535), this should produce estimated means and estimated individual values, among many other type of plots. But I got only three type of plot: histograms, scatterplots, and observed individual values. If I choose individual values, I did not get lines but several points on the Y-axis.

Any suggestions?


 Bengt O. Muthen posted on Friday, January 26, 2007 - 3:23 pm
Unfortunately, the Series plot is currently not available with Type = Random, so this would have to be plotted in another program.
 Bengt O. Muthen posted on Friday, January 26, 2007 - 3:25 pm
An approximation is to use

Series = y1(x) y2(y)....;

with x, y, ... being average values (ages) at that occasion.
 Shige Song posted on Friday, January 26, 2007 - 4:17 pm
Hi Bengt,

Thanks for clearing things up. Have a good weekend!

 Sylvana Robbers posted on Tuesday, January 29, 2008 - 1:52 am
Dear Dr. Muthen,

Could you please explain how to plot Shige's example in another program?
Which program would you suggest to use?


 Linda K. Muthen posted on Tuesday, January 29, 2008 - 6:44 am
I am not familiar with any other program that might do this.
 Sylvana Robbers posted on Wednesday, February 06, 2008 - 1:50 am
Thanks Linda.
Is maybe Bengt familiar with any program that might plot this? Are you planning to implement the Seriesplot-feature for this type of analysis in Mplus in the near future?

 Linda K. Muthen posted on Wednesday, February 06, 2008 - 10:42 am
Neither of us knows of such a program. We may add this to Mplus in the future but it is not a high priority.
 Daniel E Bontempo posted on Wednesday, August 27, 2008 - 8:35 pm
If I understand the question, you want to plot what SPSS proc MIXED calls BLUP, and Stata postestimation (e.g., after XTMIXED) predict command calls "fitted". (Both give individual model based trajectories by including fixed and random effects in the prediction). I have used savedata and FSCORES to get a file with individual scores for the growth factors and the time points. I then compute the prediction for each occasion using the individual tscore and the individual's factor score on the growth parameters. Then it is simply a matter of choosing software that handles spagetti plots well. I like Stata's twoway. Both SPSS and Excel try to make a legend entry for each person and that is a mess. SAS is also useable.

This appears to give the emperical Bayes prediction. (See discussion by SInger & Wilett, p.132). The values I get from MPlus inthis manner are identical to values the equivalent SPSS model saves under then name BLUP.

It is actually more tedious to get the group-averaged trajectories, what Stata predict calls "xb" and SPSS calls BLUE. FOr these you need to apply the fixed effects for each person, the growth factor means, and the TSCORES.

Hope this helps.
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