I am running a series of CFAs with censored variables and would like to use the MLR estimator due to missing data. I am not getting regular fit indices (chi-square, CFI, TLI, RMSEA) when I use this estimator. Do these need to be requested?
Maximum-likelihood estimation with censored variables does not have means, variances and covariances as sufficient statistics, but instead raw data, and therefore does not do the usual model test of fit. If you want a test of the fit to the covariance matrix you can use WLSMV. You can also use MLR and create your own unrestricted covariance matrix model in Mplus to test against, that is do a second run, and then compute chi-square as 2 times the log likelihood difference.
I spoke too quickly in my last sentence above. To do what I suggested one must allow a factor behind each censored variable which will lead to too many dimensions of integration with MLR.
Instead, you can either use WLSMV for such testing and hope that the missing data handling in WLSMV is sufficient. Or, better still, use MLR and work with likelihood-ratio chi-square testing of nested neighboring models to see if specific restrictions in your model are well fitting or not.
Dr. Muthen. Thank you for your prompt response. I had originally used the censored with WLSMV, but after reading more about estimators it seemed that MLR might be more appropriate. I appreciate your suggestions.
Cassie Suh posted on Monday, January 25, 2016 - 6:46 am
Hello, I am running a SEM model with a categorical variable as the final outcome and with weights.
MODEL: C by c1 c2 c3 c4 c5; D ON C B A; C ON B A; B ON A ;
I used the estimator MLR because of issues of nonnormality and because my model uses weights. I was wondering if 1) The use of MLR is appropriate for my measurement and structural model with a categorical outcome (I think it is based on other posts I've been) 2) And if so, how would I proceed to explain the results for my measurement model as the MLR estimator with categorical outcomes does not give model fit indices such as CFI, TLI, RMSEA ( I am not comparing my model with any other models). Do you think it would be enough to give results of CFA for just the one variable that is latent? 3) Also, for my structural model (I am not comparing this with any other models), is there any reference that I can use that what I am doing is correct and that no model indices of CFI, TLI and RMSEA are natural?
1. Yes. The default for MLR and the CATEGORICAL option is logistirc regression. 2-3. When numerical integration is required, means, variances, and covariances are not sufficient statistics for model estimation and chi-square and related fit statistics are not available.
Which non-normality are you concerned with? Is it your categorical factor indicators?
Cassie Suh posted on Tuesday, January 26, 2016 - 7:23 am
Thank you so much for your prompt reply- I meant that MLR gives robust standard errors.
I had two follow up questions as well.
1) just to clarify your answer to question above regarding the measurement mode, I would say that model fit indices are not available and move on.
2) I want to run a lagged effects model as well by including the same variables from previous time points as controls where BB cc1-cc5 DD corresponds to B c1-c5 D
The model would be USEVARIABLES ARE A B c1-c5 D BB cc1-cc5 DD weights; CATEGORICAL = D; WEIGHT IS weights; ANALYSIS: type=general; INTEGRATION = MONTECARLO(500); ESTIMATOR=MLR; MODEL: C by c1 c2 c3 c4 c5; CC by cc1 cc2 cc3 cc4 cc5; D ON C B A DD; C ON B A CC; B ON A BB ;
My concern is the effects of A and my final dependent variable is D.
2)-1 Is a latent variable(CC in this case) able to be used solely as a control variable like above. Or should I use CC as a observed variable(by adding/averaging cc1-cc5) and then add it to the model?
2)-2 Also, as I have a lagged effects model, I have compared it with the model with no lagged effects and got better model fot. Do you think this is nessasary as sort of a robustness check or if I should just report the results of the lagged effects models as that is the model with stronger controls and thus more appropriate?