I am running a 2-level MSEM and want to test indirect effects using Bayesian credible intervals rather than the tests provided by MODEL INDIRECT because my understanding is that the MODEL INDIRECT tests, like Sobel tests, assume that the indirect effects are normally distributed over repeated sampling, and they can produce biased significance tests as a result. However, I also have equality constraints placed on some parameters. When I try to run the model, I receive this error:
*** FATAL ERROR VARIANCE COVARIANCE MATRIX IS NOT SUPPORTED WITH ESTIMATOR=BAYES. WITHIN A VARIANCE COVARIANCE BLOCK TWO OR MORE PARAMETERS ARE EQUAL. ONLY FULL VARIANCE COVARIANCE BLOCKS WITH UNEQUAL PARAMETERS ARE ALLOWED. USE ALGORITHM=MH TO RESOLVE THIS PROBLEM.
Algorithm=MH doesn't work for a two-level model.
Do you have a suggestion for how to obtains Bayes CIs for indirect effects in a 2-level MSEM featuring equality constraints?