Attachment | Size |
---|---|
cor_mats.csv [6] | 4.6 KB |
indirectEffect metaSEM forum.R [7] | 935 bytes |
Hi metaSEM users, hi Mike,
I've been attempting to use indirectEffect() to estimate the direct and indirect effects from a series of 3 by 3 correlation matrices. I then want to meta-analyze the resulting direct and indirect effects using either meta() from the metaSEM package or mvmeta() from the mvmeta package in R.
However, I've been running into some (potential) problems doing this. First, meta() doesn't seem to converge on a stable meta-analytic solution when I attempt to find the meta-analytic estimates of the direct and indirect effects returned by indirectEffect. I've tried using rerun() to find good start values and re-fitting the model using those start values, but so far that hasn't helped.
I've also tried fitting the same model in mvmeta(), which, to my knowledge, does arrive at a stable solution (please correct me if I'm wrong). What's strange about the results from mvmeta() is that the standard error of the meta-analytic indirect effect seems, to me, anomalously small (.0020), especially in comparison to the standard error of the direct effect (.0277). However, I'm not sure how to trouble-shoot this problem, or even if this result necessarily indicates a problem.
Would anyone be able to give me some advice? I've attached sample data and R code that reproduces the problems for reference.