

Number of integration points in Monte... 

Message/Author 


Hej Bengt & Linda! I'm trying to find the optimal number of exploratory dimensions to explain covariance in a data set with 20 psychopathology items and (highly skewed) 4point ordered categorical responses to each. Though I have over 30 000 response sets, I don't want to lose any information, and missingness seems to be random, so I'm using MLR of all available categorical responses under EFA. I've tried up to a sevendimensional model, which is still interpretable, though memory constraints (16 GB) limit the adaptive integration to 3 integration points per dimension (3^7=2187). I seem to get a better fit with the Monte Carlo integration of the same number of points, so my question is this: If I want to compare e.g. the 6 and 7dimensional models with the BIC, how many integration points should I use for the adaptive MC integration? Always as many as available RAM allows, or scale it by the number of dimensions, for instance to 3^6 and 3^7 points? Cheers, Sebastian 


Sebastian I would recommend using the same number of integration points for the two runs you are comparing. Also I would recommend using the WLSMV estimator just for comparison (that estimator assumes missing completely at random) Tihomir 


Thank you for the response, Tihomir! As a related issue: Is there any way of obtaining the estimated thresholds from the MLR EFA, similarly as in the correlation matrix based estimation methods (WLSMV)? As I understand it, the program estimates (categories1) thresholds for each variable, so they are available "under the hood"? Best wishes, Sebastian 


Run the model as an ESEM model. See Example 5.24. Just leave out the covariates. 


Hi Linda, Thank you for your reply! I've run the equivalent ESEM before, but it doesn't allow numeric integration (which seems to provide slightly better fit in this problem). Is there some way of getting both ML and thresholds? Thank you in advance, Sebastian 


You can ask for TYPE=BASIC and ML and you will get the proportions. You may get those with the regular analysis. Then you can covert them to the thresholds using the following formula: logit = log (category 2/category1) 

Back to top 

