[R-meta] multivariate fixed-effect meta-analysis

Viechtbauer, Wolfgang (SP) wo||g@ng@v|echtb@uer @end|ng |rom m@@@tr|chtun|ver@|ty@n|
Wed Nov 24 14:44:43 CET 2021

>-----Original Message-----
>From: Luke Martinez [mailto:martinezlukerm using gmail.com]
>Sent: Tuesday, 23 November, 2021 20:38
>To: Viechtbauer, Wolfgang (SP)
>Cc: Filippo Gambarota; R meta
>Subject: Re: [R-meta] multivariate fixed-effect meta-analysis
>Dear Wolfgang,
>Strictly, the model is fixed-effects multivariate (i.e., MANOVA type)
>if Filippo has one effect size per outcome, right?

I don't know what you mean by that. If you only specify V and no random effects, one could call it a multivariate fixed-effects model, just like used for example in this chapter:


Whether one has one effect size per outcome or 20 is not relevant as long as V captures the covariance between the sampling errors of the estimates.

>I mean to the extent that this is not the case, then will this model
>diverge from a fixed-effect multivariate model and become more like
>marginal models (i.e., nlme::gls() type)?

Again, I can't follow your reasoning here.

>On Tue, Nov 23, 2021 at 1:22 PM Viechtbauer, Wolfgang (SP)
><wolfgang.viechtbauer using maastrichtuniversity.nl> wrote:
>> With method="FE", 'random' is also ignored. You will see in the output that it
>says "Variance Components: none".
>> If 'cov_mat' captures the sampling error covariances, then this could be argued
>to be a fixed-effects version of a multivariate model.
>> Best,
>> Wolfgang
>> >-----Original Message-----
>> >From: R-sig-meta-analysis [mailto:r-sig-meta-analysis-bounces using r-project.org]
>> >Behalf Of Filippo Gambarota
>> >Sent: Tuesday, 23 November, 2021 20:07
>> >To: R meta
>> >Subject: [R-meta] multivariate fixed-effect meta-analysis
>> >
>> >Hi!
>> >I'm performing a multivariate meta-analysis with metafor, however I'm
>> >not sure how to obtain the fixed-effect version. Given that I have not
>> >enough data I'm not interested in estimating tau for each outcome and
>> >the correlation among outcomes but only taking into account the
>> >sampling error dependence. I'm using this function:
>> >```
>> >rma.mv(
>> >    yi = eff_size,
>> >    V = cov_mat,
>> >    mods = ~ 0 + outcome,
>> >    struct = "UN",
>> >    random = ~ outcome|paper_id,
>> >    method = "FE",
>> >    data = data)
>> >```
>> >Of course, the struct argument is no more relevant (as the warning
>> >message said) but I'm wondering if the result is what I'm looking for
>> >because from the rma.mv documentation the method = "FE" is not
>> >mentioned combined with a multivariate parametrization.
>> >Thank you!
>> >
>> >--
>> >Filippo Gambarota
>> >PhD Student - University of Padova
>> >Department of Developmental and Social Psychology
>> >Website: filippogambarota.netlify.app
>> >Research Group: Colab   Psicostat

More information about the R-sig-meta-analysis mailing list