[R-meta] MLMA - shared control group
rnorouz|@n @end|ng |rom gm@||@com
Sat Aug 28 17:31:27 CEST 2021
Please see my answers below.
> Hey everyone.
> Regarding MLMA, due to shared control group, I wonder if:
> 1) Is it enough to code "studies/obs" and we are done?
> res_mlma <- rma.mv(yi, vi, random = ~ 1 | studies/obs, data=dat)
Unfortunately, no, using random-effects alone doesn't directly account
for that source of dependency. See:
for a good discussion on this.
> 2) Or after that, do we also need to compute a correlation matrix? I got
> lost in this part.
This type of dependency needs to be specified in the rma.mv() via the
V argument. See the link in the previous answer for details. Also
check out the archives to find several discussions on this.
> 3) When coding for "studies/obs", the best option is to NOT split the
> number of participants in obs?
Not sure, what you mean here, but `obs` usually denotes the id for
each unique row in your data, like:
When you fit a model via rma.mv() and specify the random part as
"studies/obs", then, a unique random effect for each study and a
unique random effect for each row within a study is added to your
model. The former accounts for the effects' variation between studies,
the latter accounts for effects' variation within studies.
> 3.1) Any good literature to support that decision in MLMA? It still seems
> strange to me, as it will inflate the actual real number of participants.
see my previous answer.
> Thanks for your time and best wishes,
> [[alternative HTML version deleted]]
> R-sig-meta-analysis mailing list
> R-sig-meta-analysis using r-project.org
More information about the R-sig-meta-analysis