Message/Author 


Dear Linda and Bengt, I estimate a twolevel model with mediation relationship on the macrolevel. I use Bayesian estimator. It looks like the "Model indirect" statement is not available in this case. Therefore, I wonder: 1. As far as Bayesian analysis is concerned, would it be appropriate to calculate the indirect effects in a usual way by multiplying the path coefficients. 2. If yes, how can I estimate whether the coefficients for the indirect effects are significant? Can I use Delta parameterization to derive the coefficients, as Bollen describes in the book or would it be something else here? Thank you, Dmitriy 


A quick followup on my own post: can I manually calculate standard errors of the estimates when using Bayesian approach? If this were the case, then I could (supposedly) calculate Sobel's test. 


What are the scales of your mediator and final outcome? 


Dear Linda, Thank you for your reply. To clarify first, I test mediation effects only on the betweenlevel. My exogenous variables are continuous. My mediators are the variables that function on both within and betweenlevels. They are modeled as categorical on the withinlevel (whether respondent considers a problem as nationally important; yes/no), but on the macrolevel they become continuous variables (proportion of the respondents in a given year who named a certain problem nationally important). Therefore, the mediators on the betweenlevel are also continuous. My final outcome is a latent variable on the between level, which is modeled by four categorical indicators. 


Because all of your variables are continuous (on Between), you can use the product formula for the indirect effect. You handle this via Model Constraint where you express the effect as a New parameter using labels for the Model parameters involved. The SEs are automatically computed by Mplus using the Delta method that Sobel used. 


I just remembered that you said you wanted to use Bayes. In this case the SEs are not computed via the Delta method but the New parameter with the indirect effect gets a posterior distribution where you can read off the 95% interval  that also takes into account any nonnormality of the indirect effect estimate distribution. 


Dear Bengt, Thank you, this helped. I did use "model constraint" statement and got the necessary estimates. I was wondering though  there is no STDYX output for the indirect effects. Of course, I can compute the standardized values manually, but I am curious why Mplus won't give them by default? 


I think you are asking wy we don't have Model Indirect yet for Twolevel  no particular reason except many other things to implement. For a New parameter in Model Constraint that you create as an indirect effect product you get STDYX there. 


Not really. In my case I only get the unstandardized effects for the new parameters. This is what surprised me. This is not much relevant however, since I can calculate the standardized coefficients manually. 


So you were asking why we don't standardize New parameters in Model Constraint? That would require another level of sophistication where we would need to go back to the Model statements to keep track of what model parameter ingredients the new parameter uses and which observed and latent variables these parameters connect with, in order to standardize the New parameter correctly. 

Johan Ng posted on Wednesday, June 13, 2012  8:21 pm



Dear Bengt & Linda To continue with the above discussion, although one can manually calculate the standardized point estimates for the indirect effects, but is there a way to do so for the credibility intervals? Would standardizing the inputs do the trick? Many thanks! 


You can get CIs for the standardized estimates that you define in Model Constraint because Model Constraint also gives you SEs. 

Johan Ng posted on Thursday, June 14, 2012  6:51 am



Dear Bengt Thanks for your swift response (as always!). Apologies for omitting a piece of important info in my previous question  I was using the BAYES estimator. Sorry for not making that clear previously. So the phrase my question again: When using the BAYES estimator, together with Model Constraint, is there a way to get standardized credibility intervals? Thanks again! 


If you use Bayes you always get CIs, also for parameters that you define in Model Constraint, including parameters that you define as standardized. It's a good way to take care of nonnormal distributions for indirect effects, for instance. 

MT posted on Monday, August 27, 2012  7:58 am



Dear Linda and Bengt, I also estimate a twolevel model with mediators using Bayes estimation. In this model, I would like to test the indirect effect of X on Y via M1 and M2. Can this be computed by simply using model constraint: indirect = a*b*c? Thank you very much! 


If your model is x > m1 > m2 > y this is correct in some cases. See the following paper on the website for details: Muthén, B. (2011). Applications of causally defined direct and indirect effects in mediation analysis using SEM in Mplus. 


Dear Mplus team, I'm running a twolevel path model using Bayes estimation. I have multiple parallel mediators on the within level and the mediations are also only modeled on the within level. However, when I specify the indirect effect via the usual MODEL CONSTRAINT command, the new parameters for the indirect effects appear in the between level and not in the within level section of the output. What does that mean? On which level are they modeled then? I'll be glad for any advice! Best regards, Anika. 


Please send the output and your license number to support@statmodel.com. 

Back to top 