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; VARIABLE: NAMES ARE id x m y; USEVARIABLES ARE id x m y; CLUSTER IS id; ANALYSIS: TYPE IS TWOLEVEL RANDOM; MODEL: %WITHIN% 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 %BETWEEN% 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
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.
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?
VARIABLE: NAMES ARE myid pol anger attackm; MISSING ARE ALL (-99); USEVARIABLES ARE myid pol anger attackm; CLUSTER IS myid; ANALYSIS: TYPE IS TWOLEVEL RANDOM; !ALGORITHM = INTEGRATION; !INTEGRATION = MONTECARLO; MODEL: %WITHIN% sa | anger ON pol; sb | attackm ON anger; sc | attackm ON pol; %BETWEEN% 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); MODEL CONSTRAINT: 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 OUTPUT: TECH1 TECH8 CINTERVAL;
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)")?
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.
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?
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?
(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.
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: WITHIN 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: %WITHIN% S | m on x1; m on x2 x3; y1 on m x1 x2 x3 y2 on m x1 x2 x3
%BETWEEN% 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)?
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: %WITHIN% social on x; %BETWEEN% social on x;
Then doing the same with random effect: %WITHIN% sa|social ON x; %BETWEEN% 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: %WITHIN% sa|social ON x; %BETWEEN% 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?
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
MODEL CONSTRAINT: 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
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.
WITHIN = u; BETWEEN = (S_ID) w (C_CN) z; CLUSTER = C_CN S_ID; ANALYSIS: TYPE = THREELEVEL RANDOM; MODEL: %WITHIN% m ON u x (aw); y ON m (bw); y ON u x; %BETWEEN S_ID% y m x; m ON w x (ab1); y ON m (bb1); y ON w x; %BETWEEN C_CN% y m x; m ON z x (ab2); y ON m (bb2); y ON z x; MODEL CONSTRAINT: NEW (indb2 indb1 indw); Indw = aw*bw; Indb1 = ab1*bb1; Indb2 = ab2*bb2; OUTPUT: TECH1 TECH8 CINTERVAL;