I am doing missing data study with Bayesian based analysis (Multiple imputation). For this simulation study, MplusAutomation of R program is also imposed. I successfully created and run models. Problem is that the model fit statistics (e.g., Posterior fit statistics and DIC) cannot read within R console. When reading Model output (readModels in R), there were only parameter estimate information. Also, when running summaryStats in R, I failed to obtain any model fit information. But, I found all fit results in the separated Mplus outputs.
Do you have any idea about it? Am I wrong? in what respect?
Second question is "PPP can be interpreted same as p-value of chi-square in ML analysis?" Can they be comparable?
I checked with the author of MplusAutomation and he had not added the fit statistics. He is doing this now and it will be in the next version of MplusAutomation which will be available next week. You can contact the author directly if you have further questions.
Yes, the p-value of PPP can be interpreted the same as the chi-square p-value. For further information, see on the website under Papers:
Muthén, B. & Asparouhov, T. (2011). Bayesian SEM: A more flexible representation of substantive theory. Forthcoming in Psychological Methods.
It is hard to diagnose the problem you described without seeing the Mplus output file or the structure of the R data read in by MplusAutomation. But if you're extracting summary statistics from a multiple imputation analysis, then most of them are reported as means and SDs. So, I would do this first:
to see what's in the data.frame. I would guess you'll see something like:
Classes 'mplus.summaries' and 'data.frame': 1 obs. of 38 variables: $ Mplus.version : chr "6" $ AnalysisType : chr "complex" $ DataType : chr "IMPUTATION" $ Estimator : chr "MLR" $ Parameters : num 15 $ LL_Mean : num -12270 $ LL_SD : num 26.5 $ LL_NumComputations : num 5 $ AIC_Mean : num 24570 $ AIC_SD : num 53.1 $ AIC_NumComputations : num 5 $ BIC_Mean : num 24665 $ BIC_SD : num 53.1
If that's the case, you should be passing things like "DIC_Mean" and "BIC_Mean" to showSummaryTable.
If this doesn't solve the problem, please email me the output file and all of the relevant R code.