Plotting an Interaction PreviousNext
Mplus Discussion > Growth Modeling of Longitudinal Data >
Message/Author
 Andrew Triplett posted on Wednesday, October 19, 2016 - 1:55 pm
I am writing for advice regarding graphing significant interactions terms in equations using latent growth curve modeling (LGCM)

Here are some details:

We are using 4 waves of data from a national survey; these data were collected over a 14 year period at unequal intervals.

The key measures for these analyses include chronic health conditions, self-rated health, perceived financial status, poverty status, and life satisfaction. These variables are assessed consistently across the four data collections.

We are examining whether life satisfaction moderates the negative effects of perceived financial status and actual poverty level on our health related outcomes. We used LGCM within MPlus version 7.4 and assessed individual changes across time. In this analytical approach, repeated measures data were used as indicators of latent variables to model baseline levels (measured as intercepts) and rates of change (measured as slope factors) to describe differences between individual trajectories.

Do you have suggestions as the easiest way to calculate a range of representative values as a way of representing the interactions?
 Bengt O. Muthen posted on Wednesday, October 19, 2016 - 5:47 pm
First, are your variables that interact time-invariant covariates predicting the growth factors? Have you looked at UG ex 3.18 to see how interactions can be plotted? See also Chapter 1 in our new book.
 Andrew Triplett posted on Friday, October 21, 2016 - 6:00 am
Thank you for your quick reply. No they are not time-invariant. Each of the variables was measured at each time point and was able to vary over time. For example, we are looking to see if life satisfaction over time can buffer the effects of changing perceived financial status has on health outcomes.

We have found significant main effects as well as well as a significant interaction term, we are just trying to determine how we can visually show this. Thanks!
 Bengt O. Muthen posted on Friday, October 21, 2016 - 12:53 pm
Ok, so you are interested in the regressions at each time point of the outcome on the tvc's and their interactions. So say that y is the outcome, x is financial status, z is poverty, and w is the life satisfaction moderator,

y = b0+b1*b2*z+b3*w+b4*x*w+b5*z*w+e
= b0+(b1+b4w)x+(b2+b5w)z+b3w+e

You can express (b1+b4w)x and (b2+b5w)z in Model Constraint using LOOP and PLOT as on our Mediation page:

http://www.statmodel.com/download/Loop%20plot%203.18.pdf

This gives a graph of the interaction plotted against the moderator.
 Andrew Triplett posted on Tuesday, November 01, 2016 - 12:01 pm
Dr. Muthen,

Thank you for your response! One last question, is there an option for obtaining/saving predicted values from a model as is the case in other software programs for linear and logistic regression? And then one could access the saved file? And if yes, would it possible to run 4 sub-group models?

Thank you!
 Bengt O. Muthen posted on Wednesday, November 02, 2016 - 5:20 pm
Yes - look for that plot option in the plot menu when using Type=Plot3. You need version 7.4.
 Andrew Triplett posted on Monday, November 14, 2016 - 9:48 am
Hi Dr. Muthen,

I am attempting to do this plot to look at the interaction, but I am getting errors in my gh5 file when trying to read it into R. Could you please check my syntax and let me know if there is something I am doing wrong?

Variable:
names = ID Sex Age Region Educ LS1 LS2 LS3 LS4 FS1 FS2 FS3 FS4 ITN1 ITN2 ITN3
ITN4 Health1 Health2 Health3 Health4 HProb1 HProb2 HProb3 HProb4;

usevariables = LS1 LS2 LS3 LS4 FS1 FS2 FS3 FS4 HProb1 HProb2 HProb3 HProb4
Sex Age Region Educ;

Missing = all (999.00);

Analysis: TYPE = RANDOM;
ALGORITHM = INTEGRATION;
ESTIMATOR = ML;

Model:
f1 BY LS1@0 LS2@.69 LS3@.77 LS4@1;
f2 BY FS1@0 FS2@.69 FS3@.77 FS4@1;
f3 BY HProb1@0 HProb2@.69 HProb3@.77 HProb4@1;
f3 ON f1 (b1)
f2 (b2);
f1xf2 | f1 XWITH f2;
f3 ON f1xf2 (b3);
f3 ON Sex Age Region Educ;

PLOT: TYPE = PLOT3;
SERIES IS LS1-LS4(*) | FS1-FS4(*);

Output:
sampstat residual standardized TECH1 TECH8;
 Bengt O. Muthen posted on Monday, November 14, 2016 - 10:30 am
Please send output and .gh5 file to Support so we can see if there is a problem on the Mplus end of this.
Back to top
Add Your Message Here
Post:
Username: Posting Information:
This is a private posting area. Only registered users and moderators may post messages here.
Password:
Options: Enable HTML code in message
Automatically activate URLs in message
Action: