TSSEM_stage2: issue with CIs. And a multiple effect sizes for the same variable

Hi Khaled,

I suggest using diag.constraints=FALSE, which is usually more robust. Since the metaSEM package uses the OpenMx package to construct the confidence intervals, there is not much I can do. I would suggest checking the recent discussion on the confidence intervals in OpenMx. If it still does not work, you may go back to the Wald confidence intervals.

Regarding the non-independent effect sizes, there is not much research on MASEM. If the amount of non-independent effect sizes is small, it seems reasonable to use the average in the analysis. If there are lots of non-independent effect sizes, either the three-level model or the robust standard error approach may be helpful. The average correlation matrix may still be used with the wls() function in the metaSEM package.

Regards,
Mike

## Use diag.constraints=FALSE rather than TRUE
random2 <- tssem2(random1, Amatrix=A, Smatrix=S, 
                  diag.constraints=FALSE, intervals="LB")
random2 <- rerun(random2)
summary(random2) 
## Inspect why it does not work
summary(random2$mx.fit, verbose = TRUE)
 
## Use Wald CI
random2b <- tssem2(random1, Amatrix=A, Smatrix=S, 
                   diag.constraints=FALSE, intervals="z")
summary(random2b)