Message/Author 


Hello, I am a new Mplus userthanks for any help or advice! I'm trying to run a measurement model + multiple paths (SEM). I've pasted my model syntax below. There are 8,037 observations (no missing data). The outcome (o1) is binary; I'm using ESTIMATOR = ML. I have three latent variables (f1, f2, f3), and all of my indicators are binary (a1a16) or ordinal/categorical (b1b9, c1c5). I am testing 4 different measured mediators (m1m4), and I have 6 measured covariates (x1x6). I have two questions: 1. MPlus has been running for almost 2 hours now. I'm wondering if this is usualperhaps given the complexity, number of dimensions in my modelor if I need to adjust any settings to help the model converge? 2. Are there any issues with combining binary indicators (latent variable f1) and ordinal indicators (latent variables f2 and f3)? E.g., would the interpretation of the coefficients change for these different factors? Thanks, Travis MODEL: f1 BY a1a16; f2 BY b1b9; f3 BY c1c5; o1 ON f1 f2 f3 m1m4 x1x6; m1 ON f1 f2 f3 x1x6; m2 ON f1 f2 f3 x1x6; m3 ON f1 f2 f3 x1x6; m4 ON f1 f2 f3 x1x6; f1 ON x1x6; f2 ON x1x6; f3 ON x1x6; 


1. How many dimensions of integration did your screen printing say that you have? 2. No problem doing this. The coefficients are in the logit scale so your interpretations of them are in line with logistic regression. 


Hi, Thank you for your reply. The screen printing says I have 3 dimensions of integration. 


That's 15*15*15 = 3375 integration points which shouldn't be too bad if you have several processors (like 8)  but it does get slower with many observations like you have. If you have a slow computer you can try integration = montecarlo(500); but keep an eye on TECH8 to make sure you don't get negative ABS changes due poor precision. 1000 points should still take only a third of the time. 


I'm a student & am running the analysis on my own laptop, which has 1 processor (dual core), so perhaps this explains the long running time. I tried your suggestion with the integration=montecarlo(500) command. This reduced the running time from 4.5 hours to 1 hour. I did check the TECH8 output. The ABS change was negative after the sixteenth step. Can you tell me how to interpret this? I also noticed that while most of the coefficients are approximately the same, some were fairly different. Thank you for your help. 


Negative ABS changes are a sign of numerical imprecision, which is in turn due to a small number of integration points (500 here). Increase to 1000 and see how that comes out. But you could use a faster computer for these sorts of analyses. Also, Bayesian estimation may be faster. 


I'll give it a try. Thank you again for your help. 

mboer posted on Saturday, April 25, 2020  4:58 am



Dear Prof. Muthen, I would like to run a multilevel logistic regression model using 5 datasets with imputed data. I also have 3 dimensions of integration, which takes a lot of computation time (especially since I the estimation runs 5 times). I have used 'integration = montecarlo(1000)' to speed up the computation time. In 2 out of 5 replications, the ABS change was negative. Are my results unreliable? Do you recommend increasing the integration points? Thank you in advance. 


See the FAQ on our web site: TECH8 – negative ABS changes Often 5000 points are needed. 

Back to top 