Message/Author 


Hello, I am doing a simulation study. I simulated the data (1000 data files) external to Mplus. I would like to run mixture modeling on each of the files and then save the results of each of the simulated data files in one summary file with individual statistics (no averages). Furthermore, I only want to save the BIC. Is there a wait to specify such. I believe the Runall will care for saving the information in one file. However, is there anyway to specify what to save? Thank you. 


I don't know of a way to do this in Mplus. See the website under HowTo in the left column. You may be able to do it via R. 


Thank you. Is there a Runall compatible with Windows 7? One last question. I cannot get the savedata option to work. Can you please correct coding? Thanks. Title: Model One Data: File=DC9Test.txt; Type=Montecarlo; Variable: Names are X1X4; Classes=c(2:4); Analysis: Type=mixture; Starts 100 15; Model: %Overall% F1 by X1  X4; Savedata: file is "results.txt"; save estimates; 

Stata posted on Thursday, June 09, 2011  12:13 pm



Can I use ordinal variable instead of dichotomous variable for latent class factor analysis? Thank you. 


Renee: We no longer support RUNALL. Please see the user's guide for the correct specification of the ESTIMATES option. 


Stata: Yes. 


None of the save options seem to work with the type=montecarlo. Is there anyway possible to save the results from the individual runs (in one file preferably or several files) instead of just the averages that are given by default. Thank you. 


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


Dear Drs. Muthen, I am running a Monte Carlo Simulation with a LTA model to test pathways between covariates and transitions between classes. The link function is logit, so that means I am testing a log odds, right? I can then exponentiate that value to present an odds ratio, I assume. Thanks, Mary 


Yes. 


Thank you! 


I have another question on this Monte Carlo simulation that I'm running. I have 3 classes and 3 points in time with a binary covariate. I want to compare power estimates between a covariate that has a 50/50 split compared to one with a 90/10 split. There is a line of code: [cg#1*.5]. Does this signify that the binary covariate has a mean of .5, therefore at 50%? Or, do I alter the line of code that has [g$1@15] and change the threshold? Then, the problem is that I don't know what thresholds match with a 90/10 split. Thanks for your help! Mary 


See the description of the Cutpoints option in UG ex 12.1. 


Hi Dr. Muthen, Thanks for your response. I am a little confused, though, about implementing the cutpoint code. My covariate or grouping variable, is defined by the code in Example 8.13 for a binary covariate. I have under the classes command: cg(2). How does this mesh with using the covariate, g, as a continuous variable with one cutpoint at 1 standard deviation above the mean? Thanks, Mary 


Then see the MonteCarlo version of ex 8.13 on our website at http://www.statmodel.com/usersguide/chapter8.shtml 

Back to top 