Monte Carlo model for unbalanced 3-le... PreviousNext
Mplus Discussion > Multilevel Data/Complex Sample >
 Fredrik Falkenström posted on Tuesday, March 24, 2015 - 4:21 am
Dear professors Muthén,

I am trying to run a Monte Carlo simulation for a three-level model, with repeated measurements at Level-1, patients at Level-2, and therapists as Level-3. I would like to be able to specify different cluster sizes both for Level-2 (i.e. patients attending different numbers of sessions) and for Level-3 (therapists treating different numbers of patients). Is this possible in Mplus? I haven't been able to find any examples of such analyses?

Best wishes,

Fredrik Falkenström
 Linda K. Muthen posted on Tuesday, March 24, 2015 - 10:30 am
See pages 783-784 of the user's guide. The NCSIZES and CSIZES options for TYPE=THREELEVEL are described there. See also the Monte Carlo counterparts of the THREELEVEL examples in the user's guide.
 Fredrik Falkenström posted on Thursday, March 26, 2015 - 10:40 am
Dear Linda, thank you very much. A quick follow-up, is there a limit for how many cluster sizes that can be specified? It seemed to me that I got problems when I got over 60, but it may be that I had done something wrong.


 Linda K. Muthen posted on Thursday, March 26, 2015 - 11:23 am
Please send the output and your license number to
 Fredrik Falkenström posted on Wednesday, April 01, 2015 - 2:46 am
Thank you so much for all your help! Do you know if there is any way of using the Monte Carlo simulation procedure to test for bias and power in comparing two models, e.g. using the Likelihood Ratio test or just the difference in AIC/BIC?


 Tihomir Asparouhov posted on Wednesday, April 01, 2015 - 10:20 am
You can use the RESULTS= option of the MONTECARLO command to save AIC/BIC/log-likelihood across the different replications. You can also use the same Model Population to generate the data in two different runs with two different Model commands. Once you get the results for the two models for all the replications you can manipulate them, ex. compare BIC etc, using excel spreadsheet or another program.
 Fredrik Falkenström posted on Wednesday, April 01, 2015 - 11:45 pm

 Fredrik Falkenström posted on Friday, April 10, 2015 - 5:50 am
Hi again, I am trying to estimate a misspecified model in which one covariate (treatment length) is used in the Model Population part but not in the estimation model. However, my model doesn’t converge when I do this, and I am thinking I may have done something wrong. My input is too long to paste here, but what I did was just to omitting this covariate from the estimation model. Is there some other way of doing this?


Fredrik Falkenström
 Bengt O. Muthen posted on Friday, April 10, 2015 - 6:54 am
Are you fixing its coefficient at zero?
 Fredrik Falkenström posted on Friday, April 10, 2015 - 1:36 pm
I had not done that since I thought leaving it out would have the same effect, but now I tested to explicitly fix its coefficients at zero and that didn't seem to help.

 Bengt O. Muthen posted on Friday, April 10, 2015 - 3:37 pm
Please send to support.
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