1-1-1 multi-level mediation with a mo... PreviousNext
Mplus Discussion > Multilevel Data/Complex Sample >
 Tsachi Ein-Dor posted on Thursday, February 10, 2011 - 1:44 am
I have a 1-1-1 mediation (see below an MSEM code). My problem is that I have a level 2 moderator as well. As far as I could understand, level 2 measures can influence level 2 constructs (between), but not level 1 (within). Am I right so far? If so can someone help with the implementation of that into the code below?

TITLE: 1-1-1 mediation (MSEM)
DATA: FILE IS mydata.dat;
id x m y;
id x m y;
m ON x(aw); ! regress m on x, call the slope "aw"
y ON m(bw); ! regress y on m, call the slope "bw"
y ON x; ! regress y on x
x m y; ! estimate Level-2 (residual) variances for x, m, and y
m ON x(ab); ! regress m on x, call the slope "ab"
y ON m(bb); ! regress y on m, call the slope "bb"
y ON x; ! regress y on x
MODEL CONSTRAINT: ! section for computing indirect effects
NEW(indb indw); ! name the indirect effects
indw=aw*bw; ! compute the Within indirect effect
indb=ab*bb; ! compute the Between indirect effect
OUTPUT: TECH1 TECH8 CINTERVAL; ! request parameter specifications, starting values,
! optimation history, and confidence intervals for all effects

 Bengt O. Muthen posted on Thursday, February 10, 2011 - 11:20 am
I am not into the numbering approach of characterizing two-level mediation modeling, but looking at Bauer et al (2006) in Psych Methods Figure 2, bottom, it looks to me that 1-1-1 mediation has a mediation model on the Within level and unrestricted random effects on the Between level. You pose a model for the random effects on Between. It happens to be a saturated model so it becomes the same as unrestricted.

Then you say "level 2 measures can influence level 2 constructs (between), but not level 1 (within)." Let's scrutinize the second part of that statement. It sounds like you want "2-2-1" as in top of Figure 2. Using an education example, you have a teacher background variable influencing a teacher dependent variable (2-2) which you want to influence a student dependent variable (2-2-1). This is done by letting the teacher dependent variable influence the between-level part of the student variable on Between, that is the random intercept (varying over teachers/classrooms) of the student variable. The model you have written in your Mplus input does this with both M and Y having between-level random intercepts.

I don't see any moderator in your Mplus input. A moderator is an interaction.
 Bengt O. Muthen posted on Thursday, February 10, 2011 - 11:51 am
You may also want to take a look at the Topic 7 handout on our web site. Slides 44-45 talk about random intercepts and random slopes (cross-level interaction).
 Sean Lane posted on Wednesday, August 24, 2011 - 9:47 am
I'm trying to fit a model similar to that listed above (using the Bauer et al syntax), but I get an error message saying that Monte Carlo integration is needed. When I specify Monte Carlo integration I get another error message saying that within-person variables cannot be on the right hand side of an ON statement. However, since I didn't specify any variables as BETWEEN I should be able to use them as WITHIN or BETWEEN, right?

myid pol anger attackm;
USEVARIABLES ARE myid pol anger attackm;
sa | anger ON pol;
sb | attackm ON anger;
sc | attackm ON pol;
sa sb sc pol anger attackm;
sa WITH sc pol anger attackm;
sa WITH sb(cab);
sb WITH sc pol anger attackm;
sc WITH pol anger attackm;
anger ON pol(ab);
attackm ON anger(bb);
attackm ON pol;
[sa](aw); [sb](bw);
NEW(a b indb indw);
a=aw+ab; ! compute Between a path
b=bw+bb; ! compute Between b path
indw=aw*bw+cab; ! compute the Within indirect effect
indb=a*b; ! compute the Between indirect effect

Any help is greatly appreciated. Thanks!
 Linda K. Muthen posted on Wednesday, August 24, 2011 - 1:05 pm
Please send the full output and your license number to support@statmodel.com.
 Anonymous posted on Friday, February 20, 2015 - 12:32 pm
Dear Drs. Muthén,
I want to run a twolevel mediation model with x, m, and y on level 1. On level 2 I have w, a predictor for the random slopes ("x on m" and "y on x").
I have two questions:


Here is my peferred solution: Would it be correct...

(step 1) to run the model with fixed slopes (l-1: "y on m x;", "m on x;") and test the indirect effect ("model indirect: y ind x")

(step 2) to include random slopes (for "m on x" and "y on x")

(step 3) to predict the random slopes with w?

Would this still be called a "moderated mediation model"?
 Anonymous posted on Friday, February 20, 2015 - 12:33 pm

Alternatively, Preacher, Zyphur, & Zhang (2010) argue that the measurement of the indirect effect is conflated/biased and recommend for the 1-1-1 model to test the indirect effect via the means of the random slopes on level 2 (see syntax H in http://www.google.de/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&cad=rja&uact=8&ved=0CCcQFjAA&url=http%3A%2F%2Fsupp.apa.org%2Fpsycarticles%2Fsupplemental%2Fa0020141%2FMET-Preacher20090012-RR-S1.doc&ei=55XnVOWoDs7iO7LtgSg&usg=AFQjCNGkr_M3CHSyYKlL5NVXXtcB9hwZ5Q&bvm=bv.86475890,d.ZWU ).

Because of the characteristcs of my data, I can only model the indirect effect on level 1, as described above.
Would there be a problem with conflation or bias? Would it be ok to use the Sobel-Test under "model indirect" and the bias-corrected confidence intervals ("output: cinterval(bcbootstrap)")?

Than you very much in advance!
 Bengt O. Muthen posted on Friday, February 20, 2015 - 2:54 pm
I would do step 1 - step 3, but in step 3 I would allow all random effects, including the 2 random intercepts to correlate freely. And, you might not find random slopes - they aren't that common - particularly not for both slopes (if they are both random special indirect formulas are needed).

You don't have to model Level-2 random effects, but can leave them free to correlate. But if you are interested in level-2 mediation processes, you can structure that as a mediation model too. The Preacher et al (2010) article is more related to that latter case.
 Bengt O. Muthen posted on Friday, February 20, 2015 - 2:55 pm
And, we ask that postings be limited to one window.
 Anonymous posted on Monday, February 23, 2015 - 8:54 am
Thank you very much for your quick response, Bengt. I have one follow up question:

In the Preacher et al. (2010) syntax document "H. 1-1-1 model (unconflated MLM)", they use groupmean centering for x and m.
I want to run a twolevel moderated mediation model. As far as I understood from Hox (2010), you should use grandmean centering for variables involved in cross-level interaction effects.

My question: Which would be the appropriate way to center my x-variable (cross-level interaction effect refers to the random slopes of x [m on x, y on x])? Would grandmean centering be appropriate or would I then run into problems regarding the indirect effect?

Again, thank you very much for your help!
 Bengt O. Muthen posted on Monday, February 23, 2015 - 9:50 am
But Preacher et al don't group-mean center a variable that has a random slope, right?
 Anonymous posted on Tuesday, February 24, 2015 - 4:41 pm
Dear Bengt, three questions:

(1) That is correct, Preacher et al. don't have random slopes. I found slides from Andrew Hayes where he applied the Preacher syntax for random slopes. He also used groupmean centering there, so I had the impression it had something to do with the indirect effect.
I now believe to understand that the groupmean centered x and m are used to correct the conflation in 1-1-1 models (xmean and mmean in Preacher's syntax). I think I understand this idea, but I am still confused how to analyse a cross-level interaction effect for my random slopes when I use groupmean centering - not grandmean-centering as recommended by Hox (2010).

(2) In the UG example 9.5 you analyse cross-level interaction but you don't use centering at all. Don't the variables involved in the interaction NEED to be centered? Maybe I am mistaken in that point.

(3) In your post from Friday, February 20, 2015 - 2:54 pm you talk about special indirect formulas for modeling several random slopes. In the UG example 9.5 you model two random slopes using a very basic syntax. What would be the alternative syntax and do I really need one?

Thank you so much for your help.
 Bengt O. Muthen posted on Wednesday, February 25, 2015 - 1:18 pm
(1) Do you mean that you are confused by how to do it in Mplus, or about the meaning of it? For the latter, I recommend turning to Multilevelnet.

(2) They don't NEED to be centered. but as with any interactions it can help computationally (reducing collinearity) and in interpretations (see the example on slide 33 and on of our Topic 7 handout and video from our Johns Hopkins Mplus short courses on our website).

(3)In UG ex 9.5 we don't estimate an indirect effect. To do that, see slide 80 of the Topic 7 handout and video.
 Annamaria Kubovcikova posted on Monday, May 09, 2016 - 1:58 am

I have a very similar problem to the person above. I also want to run a twolevel mediation model. I however have x1, x2, x3 and they are all mediated by the same m to two dependent variables y1 and y2.
So my model is:
m on x1 x2 x3
y1 on m x1 x2 x3
y2 on m x1 x2 x3

I have managed to extend the Preacher, Zyphur, & Zhang (2010) unconflated approach for my model (their model "I" with fixed effects). That means basically estimating the same set of mediations both WITHIN and BETWEEN.

I suspect one of the mediations to be random.
Under the classical "conflated" approach I was able to make the slope between x1 and m random and then explain its variance between in the following way:
S | m on x1;
m on x2 x3;
y1 on m x1 x2 x3
y2 on m x1 x2 x3

S on w;
S WITH y1 y2;
But how do I do this with the unconflated approach?
I can fit one random mediation as Preacher, Zyphur, & Zhang (2010) suggest.
They propose that all the paths are random between "y on m" "y on x" "m on x".
Can I make just one slope "m on x1" random?
Bengt, you have mentioned that random slopes are very seldom, how do I test for their randomness (so I have justification for freeing just certain slopes and keeping the rest fixed)?

Thank you!
 Bengt O. Muthen posted on Tuesday, May 10, 2016 - 6:44 pm
First see if you have a random slope by just doing the 2-level regression for the mediator as DV.
 Annamaria Kubovcikova posted on Thursday, May 12, 2016 - 6:04 am
Hello Bengt,
Thank you very much for your quick response. Keeping my unconflated ambitions, I have tried a very simple analysis as you advised. First fixed effect model, using the mediator (social) as DV:
social on x;
social on x;

Then doing the same with random effect:
sa|social ON x;
social on Dane;
sa WITH social Dane;

In the random model, in the between part of output I get the variance for both sa and Dane, the variance for Dane is significant the variance for sa is insignificant. Looking at the likelihood ratio test, 2x(-1614.136+1612.121)=4.03 and 6 degrees of freedom, the test is insigificant. SO there should be no random coefficient. Then why can I find significant effect when I try to explain the random slope in this model:
sa|social ON x;
social on x;
sa on moderator (varies only between);
sa WITH social x;

It seems like I am explaining non-existing random slope.
What am I doing wrong?
 Annamaria Kubovcikova posted on Thursday, May 12, 2016 - 6:07 am
Sorry I am posting in one more window, but one more thing might be relevant, all my variables are uncentered, the variable social is continuous, but x and moderator are dummy variables (0/1).

Thank you again!
 Bengt O. Muthen posted on Friday, May 13, 2016 - 1:35 pm
It is a common finding that when you add a covariate to predict a random slope you get significance. I would think it is due to increased power.
 Ok-young, JI posted on Monday, January 16, 2017 - 11:35 pm
hi Dr. Muthens
i have a multilevel data and my model consits of 1-1-1-1 mediation (all within level).
i was wondering that my estimates are all fixed effect then sould i call these estimates are beta or gamma?
And could you pleas let me know about my model's mediation equation????
here is my syntax

title: this is twolevel analysis with day and level variables
data: file is small model9_final.txt ;
variable: names = x1 x2 m1 m2 y1 y2 y3 w id ;
usevariables = x1 x2 m1 m2 y1 id ;
within = x1 x2 m1 m2 y1 ;
between = ;
cluster = id ;
analysis: type=twolevel ;
estimator = ML ;
model: %WITHIN% ! model for within effects / fixed effect
m1 ON x1 (aw1) ; ! x1 -> m1
m1 ON x2 (aw2) ; ! x2 -> m1
m2 ON m1 (bw) ; ! m1 -> m2
y1 ON m2 (cw) ; ! m2 -> y1
y1 ON x1 ; ! x1 -> y1 direct effect
y1 ON x2 ; ! x2 -> y1 direct effect
%BETWEEN% ! model for between effects / fixed effect

NEW(indw1 indw2) ;
indw1=aw1*bw*cw ; ! compute the within level indirect effect1 for m_WRT ->m_PD->m_SED->e_JE
indw2=aw2*bw*cw ; ! compute the within level indirect effect2 for m_WSA ->m_PD->m_SED->e_JE

pleaese, help me!!
 Bengt O. Muthen posted on Tuesday, January 17, 2017 - 9:07 am
Please send your output to Support along with your license number.
 Derek Boy posted on Monday, January 15, 2018 - 9:22 pm
Hi, Dr. Muthen

I need to analyze a multilevel mediation scenario involving six variables (z, w, u, x, m, y), having outcome y modeled at three levels of L1, L2, and L3, mediator m at all three levels, predictor x at all three levels, predictor u at only L1, predictor w at only L2, and predictor z at only L3.

I actually wrote my codes basing on Preacher, Zyphur, & Zhang’s (2010, 2011) accompanying material, “Mplus syntax files for single- and multilevel mediation models: 1-1-1 model with fixed slopes (MSEM)”, which had originally intended for a TWOLEVEL RANDOM model, but I adapted and expanded it to suit my need of the above-mentioned THREELEVEL RANDOM situation.

The codes I wrote are as below. Please kindly let me know if I am doing it right or not. And, my main interest is the indirect effects between y, m, and x.

BETWEEN = (S_ID) w (C_CN) z;
m ON u
x (aw);
y ON m (bw);
y ON u x;
y m x;
m ON w
x (ab1);
y ON m (bb1);
y ON w x;
y m x;
m ON z
x (ab2);
y ON m (bb2);
y ON z x;
NEW (indb2 indb1 indw);
Indw = aw*bw;
Indb1 = ab1*bb1;
Indb2 = ab2*bb2;
 Bengt O. Muthen posted on Tuesday, January 16, 2018 - 2:29 pm
This looks fine.
 Derek Boy posted on Tuesday, January 16, 2018 - 4:02 pm
Thanks for your kind comment.
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