Moderation Analysis PreviousNext
Mplus Discussion > Structural Equation Modeling >
 learningmplus posted on Saturday, October 22, 2011 - 3:44 pm
Hi Linda & Bengt!

I'm hoping you'll be able to help me. I am looking at the interaction between type of treatment (COND = family = 1; group = 0) on anxiety after treatment. In my model I have
Estimate P-Value

COND -0.444 0.471
DMARITAL 0.847 0.428
DMALE -0.516 0.460
D_ETHNIC -0.795 0.409
M_RCMAS 0.580 0.000
M_AGE -0.880 0.000
M_AGEXCOND 0.938 0.003

RCMASPPO 7.606 0.000

COND 1=family, 0=group
I'm thinking that the mean of change for group = 7.606; family = 7.606+(-.444) = 7.162.
As age goes up, change in anxiety goes down (i.e., anxiety changes less).
I am not sure how to interpret the product term. The reason I am asking this question here is because I have an additional question about missing data. I got this finding using FIML in MPlus. However, when I imputed the data in SPSS and used other tools to test the interaction I got a non-significant result. I am not sure if it is okay to use FIML when testing the above. Sorry if the questions sound basic, but whoever I asked in my program did not know and I need to figure it out for my dissertation. Thank you very much.
 Linda K. Muthen posted on Sunday, October 23, 2011 - 3:20 pm
FIML and imputation are asymptotically equivalent but may differ in practice with small samples. You may not have the same sample size in Mplus and SPSS. Mplus may delete observations with missing on x.

The interaction between age and cond is signficant in the prediction of anxiety. This means that treatment behaves differently with respect to age. For cond=1, the coefficient is -.880 + .938, a small positive value. For cond=0, it is -.880, a negative value.
 learningmplus posted on Sunday, October 23, 2011 - 6:31 pm
Thank you. I get the full sample in MPlus by bringing in the covariates x-y*; as advised here. I also get the full sample in SPSS when I impute the values, but I don't get the same results. Could clustering be the reason? MPlus accounts for clustering whereas SPSS doesn't. It would be interesting to know in what way, if any, clustering makes a difference in moderation analyses.

Based on the above -- would it be accurate if I say that for cond = 1 -- for every one unit that age goes up, anxiety goes up by .06. This is a small positive value.

for cond = 0 -- for every one unit that age goes up, anxiety goes down by -.880 units.

I'm not sure if the way I wrote this makes sense. I should say that I used

Anxiety Post on Anxiety Pre
Anxiety Post on M_Age
Anxiety Post ON M_AgeXCond as the syntax.

And one more thing -- I'm running some analyses with categorical moderators. How can I get the simple main effects and the means so I put them in a 2X2 table? It would make it easier for me to interpret the findings. I will then be able to report simple main effect contrasts and interaction contrasts.

Thanks very much! This helps me a lot toward my dissertation.
 Bengt O. Muthen posted on Tuesday, October 25, 2011 - 9:03 am
If you do Type = Complex in your Mplus run to deal with clustering, the parameter estimates are not affected, only their SEs. So that wouldn't be the source of the difference visavis SPSS.

Yes, your statements about how age affects anxiety are correct.

You should also have Cond on the right-hand side of your ON statements, which you did in your initial post.

You ask about means in a 2 x 2 table, so you must be asking about a binary moderator with your binary Cond variable. This is regression with two binary covariates. It is straightforward to get the 4 means by plugging in the 0's/1's for the two variables as in regular regression.
 Anke Schmitz posted on Wednesday, February 19, 2014 - 10:17 am
I have a signifikant moderation between two observed variables (one dichotomous and one continous). Where can I see how large the effect is for students with mean values, -1SD and +1 SD?
Regards Anke
 Bengt O. Muthen posted on Wednesday, February 19, 2014 - 11:34 am
See UG ex 3.18.
 Anke Schmitz posted on Thursday, February 20, 2014 - 7:07 am
Dear Bengt,
I have huge problems to apply ex 3.18 on my model.

Here is my command:

DEFINE: int1 = t*g; !Interaction between continuous observed variable an binary observed categorical variable


LV by lv4 lv7 lv9 lv10 lv11 lv13 lv14 lv17 lv21;

VW by vw2 vw3 vw5 vw6 vw7;

int2 | vw XWITH text;

LV on VW t g int1 int2;

There is a significant interaction of int1. I have to know how the effect is when the continuous varibale is -1SD, Mean and +1SD. I don't understand the command in ex. 3.18.
 Bengt O. Muthen posted on Friday, February 21, 2014 - 8:31 am
The LOOP plot is described in the UG on page 696. An application is shown on slides 21-24 in the V7Part1.pdf handout shown at

See Part 1 under

Handouts for New Developments in Mplus Version 7
The handouts for the Mplus Version 7 workshops at Utrecht University on August 27-29, 2012 are posted here in 4-per-page format and in regular format:

Part 1: 4-per-page Regular
Part 2: 4-per-page Regular
Part 3: 4-per-page Regular

You also find a video of that there.
 Anke Schmitz posted on Saturday, February 22, 2014 - 6:35 am
Bengt, thank you very much! It worked quite well. I have one more question: How do I compare a model without interaction (just with main effects) with a model with interaction term?
There are no model criteria in the output. Just those information:

Number of Free Parameters 36


H0 Value -6874.044
H0 Scaling Correction Factor 1.0009
for MLR

Information Criteria

Akaike (AIC) 13820.088
Bayesian (BIC) 13987.030
Sample-Size Adjusted BIC 13872.714
(n* = (n + 2) / 24)

Can I do it with the formula on this website?

Kind regards
 Bengt O. Muthen posted on Saturday, February 22, 2014 - 3:54 pm
Since it is only 1 parameter that differs between the 2 models you can simply use the z test that is given in the output. Using MLR, this takes the non-normality into account.
 Anke Schmitz posted on Thursday, March 06, 2014 - 7:39 am
Thanks for your replys on feb, 22. It worked quite good. Is there any chance to compare moderationsmodels between two observed groups (1 versus 2)? I just want to know if the interaction effects differ between two groups of students or if they are similar. I tried the "grouping is" command, but it did not work. Do I have to split the data file and generate two separate models?
Regards Anke
 Linda K. Muthen posted on Thursday, March 06, 2014 - 2:02 pm
Please send the output showing GROUPING not working and your license number to
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