Message/Author 


I would like to test for mediation using the INDIRECT, CINTERVAL and BOOTSTRAP commands. But I have a complex sample. Will the Bootstrap function properly with weighted clustered data? Am I correct in assuming that I don't have to use the TYPE=COMPLEX specification because the bootstrap will provide the appropriate variance estimates? Thanks. 


BOOTSTRAP cannot be used with COMPLEX. BOOTSTRAP does not take into account lack of independence of observations. You need COMPLEX for that. 

Sarah Ryan posted on Thursday, September 15, 2011  11:01 am



I am using data from the ELS:2002 data set (an NCES data set). Is it possible to use "Resampling Methods in Mplus for Complex Survey Data" (Tihomir Asparouhov and Bengt Muthen, May 4, 2010) as a guide in obtaining bootstrap standard errors? If so, would bootstrap SE's also then be provided for the MODEL INDIRECT output? Like the example given with ECLS data in this paper, the sampling structure is available in the ELS:2002 data, i.e., the strata and PSU (cluster) variable are available in the sample. 


Yes on both questions. You can obtain the bootstrap SE using REPSE=BOOTSTRAP; BOOTSTRAP=500; and "model indirect:" is available as well. 

Sarah Ryan posted on Friday, September 16, 2011  9:10 am



Excellent. Thank you. SR 

Stat posted on Saturday, April 19, 2014  1:32 pm



I try to do a mediation and compute indirect/direct effects using SEM with latent variables and a complex sample design (Cluster, Stratum and weight). Since I cannot use bootstrap with COMPLEX, I tried to use the "Resampling Methods in Mplus for Complex Survey Data". Since the p values changed slightly, I would like to be able to justify this choice. So: 1Is the “Resampling Methods [bootstrap in my case] in Mplus for Complex Survey Data" give SE that are more “accurate”, or does the idea of the resampling methods was implemented to answer specific problems. In other words, when we have complex sample design, is it generally better to use this method, or only in certain (precise) cases? 2Will it still be possible/adequate to use model constraints to compare strength of my direct and indirect effects? Thank you 


You should use MLR or ML if you are using replicate weights. The indirect effect standard errors are typically okay unless the sample size is small. Bootstrap often does not make a difference. 

Stat posted on Monday, April 21, 2014  11:19 am



Thank you! 


Just for clarification, is it or is it not incorrect to use the "repe=bootstrap" command if you do not have replicate weights? Said another way, if I have weights, strata, and clustering design effects but not replicate weights, it is appropriate to use "repse=bootstrap" to obtain bootstapped confidence intervals for indirect effects? 


This option is for replicate weights only. I would not using it for any other purpose. 


Dear Bengt, Linda & Tihomir, what would you recommend for 2 mediatian analyses on 2 levels with TYPE IS TWOLEVEL COMPLEX that is comparable to bootstrapping? I tried to calculate the indirect effects via Estimator = bayes and MODEL CONSTRAINT. However Estimator BAYES is not allowed with TYPE=COMPLEX. Also the resampling methods does not seem to work. Thank you very much Hannah 


You can use Bayes with type=threelevel and model both levels that are on the cluster command. At this point this is the only way to obtain asymmetric distributions for the indirect effects in your situation. 


Thanks a lot! 


I have data from children nested in classes so I use type = complex. Model is with indirect effects (latent and observed variables) . Models without bootstrap and with bootstrap (type =complex with WEIGHT IS weight; and weight =1) give substantially different results not only in estimation of INDIRECT effects but also DIRECT effects. For example: part of result from model WITHOUT bootstrap DEPR ON LIKE 0.247 0.208 1.189 0.234 DISLIKE 0.359 0.359 1.001 0.317 And the same WITH bootstrap DEPR ON LIKE 0.275 0.140 1.966 0.049 DISLIKE 0.421 0.224 1.877 0.061 Which result should I believe? What kind of analysis do you recommend to obtain trustworthy results? 


Bootstrap does not affect the point estimates. There must be another difference between the files (such as different weight variable which can affect the point estimates). Since this is a mixture model it maybe that one of the runs is at incomplete convergence/local maximum (but that is unlikely if you have all other estimation settings the same). To avoid that problem you can use the run with the higher likelihood to obtain good starting values for the other run (using output:svalaues in the higher likelihood run and then use these starting values for the other run with starts=0). If you can't figure it out send your example and data to support@statmodel.com 


Thank you, I add wieght to model without bootstrap (the same as in model with) but it doesn't change it's point estimates (the differences are still there). Now I'll try with svalues. 

Back to top 