[R-meta] R-sig-meta-analysis Digest, Vol 37, Issue 23

Send R-sig-meta-analysis mailing list submissions to
        r-sig-meta-analysis at r-project.org

To subscribe or unsubscribe via the World Wide Web, visit
        https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis
or, via email, send a message with subject or body 'help' to
        r-sig-meta-analysis-request at r-project.org

You can reach the person managing the list at
        r-sig-meta-analysis-owner at r-project.org

When replying, please edit your Subject line so it is more specific
than "Re: Contents of R-sig-meta-analysis digest..."


Today's Topics:

   1. Re: weight in rmv metafor (Viechtbauer, Wolfgang (SP))
   2. Re: weight in rmv metafor (Emerson Del Ponte)
   3. Re: weight in rmv metafor (Norman DAURELLE)

----------------------------------------------------------------------

Message: 1
Date: Thu, 11 Jun 2020 13:33:02 +0000
From: "Viechtbauer, Wolfgang (SP)"
        <wolfgang.viechtbauer at maastrichtuniversity.nl>
To: Norman DAURELLE <norman.daurelle at agroparistech.fr>
Cc: r-sig-meta-analysis <r-sig-meta-analysis at r-project.org>
Subject: Re: [R-meta] weight in rmv metafor
Message-ID: <1b8c1463bcdf43baaa39788aa2a859c1 at UM-MAIL3214.unimaas.nl>
Content-Type: text/plain; charset="iso-8859-1"

Dear Norman,

To give a simple example: When (some of the) studies contribute multiple estimates, the dataset has a multilevel structure (with estimates nested within studies). A common way to deal with this is to fit a multilevel model with random effects for studies and estimates within studies. Like this:

http://www.metafor-project.org/doku.php/analyses:konstantopoulos2011

However, multiple estimates from the same study are actually often computed based on the same sample of subjects. In that case, the sampling errors are also correlated. The multilevel model does not capture this. For this, one would ideally want to fit a model that also allows for correlated sampling errors. Like this:

http://www.metafor-project.org/doku.php/analyses:berkey1998

However, computing the covariances between the sampling errors within a study is difficult and requires information that is often not available.

We can ignore those correlations and use the multilevel model as a working model that is an approximation to the model that also accounts for correlated sampling errors. After fitting the multilevel model with rma.mv(), one can then use cluster robust inference methods to 'fix things up'.

Quite a bit of this has been discussed at length in previous posts on this mailing list. You might want to search the archives for some of these posts.

Best,
Wolfgang

>-----Original Message-----
>From: Norman DAURELLE [mailto:norman.daurelle at agroparistech.fr]
>Sent: Thursday, 11 June, 2020 15:05
>To: Viechtbauer, Wolfgang (SP)
>Cc: r-sig-meta-analysis
>Subject: Re: [R-meta] weight in rmv metafor
>
>Thank you.
>I am not sure I understand exactly what you mean by " if the working model
>is only an approximation and doesn't cover all dependencies ".
>Could you please explain it ?
>
>For now I used the rma() function to synthesize the available literature
>existing on the blackleg - oil seed rape disease-yield relationship, using
>slopes as effect-sizes.
>the models that gave me the slopes I used in the meta-analysis are all Y = a
>+ bX, simple linear regressions with Y being the yield and X being the
>diseqse severity.
>So my slopes, b, are all negative, and I have obtained a "summary" effect
>size through the rma() function.
>
>But I indeed have two studies that for now contribute to most of the effect-
>sizes that are included in my meta-analysis.
>
>So why exactly is it necessary to use the rma.mv() function ?
>What exactly does the "multivariate" qualificative refer to ?
>
>Thank you,
>Norman.
>
>________________________________________
>De: "Wolfgang Viechtbauer" <wolfgang.viechtbauer at maastrichtuniversity.nl>
>?: "Norman DAURELLE" <norman.daurelle at agroparistech.fr>, "r-sig-meta-
>analysis" <r-sig-meta-analysis at r-project.org>
>Envoy?: Jeudi 11 Juin 2020 22:34:55
>Objet: RE: [R-meta] weight in rmv metafor
>
>Dear Norman,
>
>If you only used rma(), then this is not correct. rma.mv() with an
>appropriately specified model (plus clubSandwich::coef_test() if the working
>model is only an approximation and doesn't cover all dependencies) would be
>more appropriate.
>
>Best,
>Wolfgang
>
>>-----Original Message-----
>>From: Norman DAURELLE [mailto:norman.daurelle at agroparistech.fr]
>>Sent: Thursday, 11 June, 2020 14:13
>>To: r-sig-meta-analysis
>>Cc: Viechtbauer, Wolfgang (SP)
>>Subject: Re: [R-meta] weight in rmv metafor
>>
>>Hi all,
>>
>>I read this discussion and one question came to my mind : I also had some
>>studies that contributed multiple effect sizes in the meta-analysis that I
>>recently ran thanks to Dr Viechtbauer's advice.
>>For now I only used the rma function, but should I have used rma.mv because
>>of these stuides that had multiple effect sizes ?
>>
>>Thank you !
>>
>>Norman
>>
>>________________________________________
>>De: "James Pustejovsky" <jepusto at gmail.com>
>>?: "Wolfgang Viechtbauer" <wolfgang.viechtbauer at maastrichtuniversity.nl>
>>Cc: "r-sig-meta-analysis" <r-sig-meta-analysis at r-project.org>, "Huang Wu"
>><huang.wu at wmich.edu>
>>Envoy?: Mercredi 10 Juin 2020 05:08:09
>>Objet: Re: [R-meta] weight in rmv metafor
>>
>>Hi Huang,
>>
>>I've written up some notes that add a bit of further intuition to the
>>discussion that Wolfgang provided. The main case that I focus on is a model
>>that is just a meta-analysis (i.e., no predictors) and that includes random
>>effects to capture both between-study and within-study heterogeneity. I
>>also say a little bit about meta-regression models with only study-level
>>predictors.
>>
>>https://www.jepusto.com/weighting-in-multivariate-meta-analysis/
>>
>>Best,
>>James
>>
>>On Sun, Jun 7, 2020 at 4:11 PM Viechtbauer, Wolfgang (SP) <
>>wolfgang.viechtbauer at maastrichtuniversity.nl> wrote:
>>
>>> Of course the weights "impact the estimated fixed effects". But whether
>>> studies with multiple effect sizes tend to receive more weight depends on
>>> various factors, including the variances of the random effects and the
>>> sampling error (co)variances.
>>>
>>> A more detailed discussion around the way weighting works in rma.mv
>>> models can be found here:
>>>
>>> http://www.metafor-project.org/doku.php/tips:weights_in_rma.mv_models
>>>
>>> Note that weights(res, type="rowsum") currently only works in the 'devel'
>>> version of metafor, so follow
>>> https://wviechtb.github.io/metafor/#installation if you want to reproduce
>>> this part as well.
>>>
>>> I hope this clarifies things.
>>>
>>> Best,
>>> Wolfgang




------------------------------

Message: 2
Date: Thu, 11 Jun 2020 10:43:05 -0300
From: Emerson Del Ponte <delponte at ufv.br>
To: Norman DAURELLE <norman.daurelle at agroparistech.fr>
Cc: Wolfgang Viechtbauer
        <wolfgang.viechtbauer at maastrichtuniversity.nl>,  r-sig-meta-analysis
        <r-sig-meta-analysis at r-project.org>
Subject: Re: [R-meta] weight in rmv metafor
Message-ID:
        <CAFgcnpJEEXJdSw0ySSTz+WhLtm+D3RgWm9Y0e8RqFLj3UezjUw at mail.gmail.com>
Content-Type: text/plain; charset="utf-8"

Dear Norman

You may want to check reproducible examples of my previous work on this
exact application context as a starting point.

https://emdelponte.github.io/paper-white-mold-meta-analysis/

https://emdelponte.github.io/paper-FHB-yield-loss/code_meta_analysis.html

Emerson




On Thu, 11 Jun 2020 at 10:06 Norman DAURELLE <
norman.daurelle at agroparistech.fr> wrote:

>
> Thank you.
> I am not sure I understand exactly what you mean by " i f the working
> model is only an approximation and doesn't cover all dependencies ".
> Could you please explain it ?
>
> For now I used the rma() function to synthesize the available literature
> existing on the blackleg - oil seed rape disease-yield relationship, using
> slopes as effect-sizes.
> the models that gave me the slopes I used in the meta-analysis are all Y =
> a + bX, simple linear regressions with Y being the yield and X being the
> diseqse severity.
> So my slopes, b, are all negative, and I have obtained a "summary" effect
> size through the rma() function.
>
> But I indeed have two studies that for now contribute to most of the
> effect-sizes that are included in my meta-analysis.
>
> So why exactly is it necessary to use the rma.mv() function ?
> What exactly does the "multivariate" qualificative refer to ?
>
> Thank you,
> Norman.
>
>
> De: "Wolfgang Viechtbauer" <wolfgang.viechtbauer at maastrichtuniversity.nl>
> ?: "Norman DAURELLE" <norman.daurelle at agroparistech.fr>,
> "r-sig-meta-analysis" <r-sig-meta-analysis at r-project.org>
> Envoy?: Jeudi 11 Juin 2020 22:34:55
> Objet: RE: [R-meta] weight in rmv metafor
>
> Dear Norman,
>
> If you only used rma(), then this is not correct. rma.mv() with an
> appropriately specified model (plus clubSandwich::coef_test() if the
> working model is only an approximation and doesn't cover all dependencies)
> would be more appropriate.
>
> Best,
> Wolfgang
>
> >-----Original Message-----
> >From: Norman DAURELLE [mailto:norman.daurelle at agroparistech.fr]
> >Sent: Thursday, 11 June, 2020 14:13
> >To: r-sig-meta-analysis
> >Cc: Viechtbauer, Wolfgang (SP)
> >Subject: Re: [R-meta] weight in rmv metafor
> >
> >Hi all,
> >
> >I read this discussion and one question came to my mind : I also had some
> >studies that contributed multiple effect sizes in the meta-analysis that
> I
> >recently ran thanks to Dr Viechtbauer's advice.
> >For now I only used the rma function, but should I have used rma.mv
> because
> >of these stuides that had multiple effect sizes ?
> >
> >Thank you !
> >
> >Norman
> >
> >________________________________________
> >De: "James Pustejovsky" <jepusto at gmail.com>
> >?: "Wolfgang Viechtbauer" <wolfgang.viechtbauer at maastrichtuniversity.nl>
> >Cc: "r-sig-meta-analysis" <r-sig-meta-analysis at r-project.org>, "Huang
> Wu"
> ><huang.wu at wmich.edu>
> >Envoy?: Mercredi 10 Juin 2020 05:08:09
> >Objet: Re: [R-meta] weight in rmv metafor
> >
> >Hi Huang,
> >
> >I've written up some notes that add a bit of further intuition to the
> >discussion that Wolfgang provided. The main case that I focus on is a
> model
> >that is just a meta-analysis (i.e., no predictors) and that includes
> random
> >effects to capture both between-study and within-study heterogeneity. I
> >also say a little bit about meta-regression models with only study-level
> >predictors.
> >
> >https://www.jepusto.com/weighting-in-multivariate-meta-analysis/
> >
> >Best,
> >James
> >
> >On Sun, Jun 7, 2020 at 4:11 PM Viechtbauer, Wolfgang (SP) <
> >wolfgang.viechtbauer at maastrichtuniversity.nl> wrote:
> >
> >> Of course the weights "impact the estimated fixed effects". But whether
> >> studies with multiple effect sizes tend to receive more weight depends
> on
> >> various factors, including the variances of the random effects and the
> >> sampling error (co)variances.
> >>
> >> A more detailed discussion around the way weighting works in rma.mv
> >> models can be found here:
> >>
> >> http://www.metafor-project.org/doku.php/tips:weights_in_rma.mv_models
> >>
> >> Note that weights(res, type="rowsum") currently only works in the
> 'devel'
> >> version of metafor, so follow
> >> https://wviechtb.github.io/metafor/#installation if you want to
> reproduce
> >> this part as well.
> >>
> >> I hope this clarifies things.
> >>
> >> Best,
> >> Wolfgang
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-meta-analysis mailing list
> R-sig-meta-analysis at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis
>
--
*Emerson M. Del Ponte*
Universidade Federal de Vi?osa, Brazil
Chair of the Graduate Studies <http://www.dfp.ufv.br/graduate/> in Plant
Pathology
EIC for Tropical Plant Pathology <http://sbfitopatologia.org.br/tpp/>
Co-Founder of Open Plant Pathology <https://www.openplantpathology.org/>
My websites: Twitter <https://twitter.com/edelponte> | GitHub
<https://github.com/emdelponte> | Google Scholar
<https://scholar.google.com.br/citations?user=a1rPnI0AAAAJ> | ResearchGate
<https://www.researchgate.net/profile/Emerson_Del_Ponte>
Tel +55 31 36124830

        [[alternative HTML version deleted]]




------------------------------

Message: 3
Date: Fri, 12 Jun 2020 01:08:31 +0200 (CEST)
From: Norman DAURELLE <norman.daurelle at agroparistech.fr>
To: Emerson Del Ponte <delponte at ufv.br>
Cc: Wolfgang Viechtbauer
        <wolfgang.viechtbauer at maastrichtuniversity.nl>,  r-sig-meta-analysis
        <r-sig-meta-analysis at r-project.org>
Subject: Re: [R-meta] weight in rmv metafor
Message-ID:
        <1316672195.29298546.1591916911202.JavaMail.zimbra at agroparistech.fr>
Content-Type: text/plain; charset="utf-8"


Dear all, Dr Viechtbauer, Dr Del Ponte,thank you for your answers ! I will look into what you advised me to read, and also go and read more of the archive.If I have more questions I will come back and ask them.Once again, thank you for developing the package metafor in R Dr Viechtbauer, and for creating this mailing list.Have a nice week-end,Norman
----- Mail d'origine -----
De: Emerson Del Ponte <delponte at ufv.br>
?: Norman DAURELLE <norman.daurelle at agroparistech.fr>
Cc: Wolfgang Viechtbauer <wolfgang.viechtbauer at maastrichtuniversity.nl>, r-sig-meta-analysis <r-sig-meta-analysis at r-project.org>
Envoy?: Thu, 11 Jun 2020 15:43:05 +0200 (CEST)
Objet: Re: [R-meta] weight in rmv metafor

Dear Norman
You may want to check reproducible examples of my previous work on this exact application context as a starting point.
https://emdelponte.github.io/paper-white-mold-meta-analysis/
https://emdelponte.github.io/paper-FHB-yield-loss/code_meta_analysis.html
Emerson



On Thu, 11 Jun 2020 at 10:06 Norman DAURELLE <norman.daurelle at agroparistech.fr> wrote:


Thank you.

I am not sure I understand exactly what you mean by " i f the working model is only an approximation and doesn't cover all dependencies ".

Could you please explain it ?


For now I used the rma() function to synthesize the available literature existing on the blackleg - oil seed rape disease-yield relationship, using slopes as effect-sizes.

the models that gave me the slopes I used in the meta-analysis are all Y = a + bX, simple linear regressions with Y being the yield and X being the diseqse severity.

So my slopes, b, are all negative, and I have obtained a "summary" effect size through the rma() function.


But I indeed have two studies that for now contribute to most of the effect-sizes that are included in my meta-analysis.


So why exactly is it necessary to use the rma.mv() function ?

What exactly does the "multivariate" qualificative refer to ?


Thank you,

Norman.



De: "Wolfgang Viechtbauer" <wolfgang.viechtbauer at maastrichtuniversity.nl>

?: "Norman DAURELLE" <norman.daurelle at agroparistech.fr>, "r-sig-meta-analysis" <r-sig-meta-analysis at r-project.org>

Envoy?: Jeudi 11 Juin 2020 22:34:55

Objet: RE: [R-meta] weight in rmv metafor


Dear Norman,


If you only used rma(), then this is not correct. rma.mv() with an appropriately specified model (plus clubSandwich::coef_test() if the working model is only an approximation and doesn't cover all dependencies) would be more appropriate.


Best,

Wolfgang


>-----Original Message-----

>From: Norman DAURELLE [mailto:norman.daurelle at agroparistech.fr]

>Sent: Thursday, 11 June, 2020 14:13

>To: r-sig-meta-analysis

>Cc: Viechtbauer, Wolfgang (SP)

>Subject: Re: [R-meta] weight in rmv metafor

>

>Hi all,

>

>I read this discussion and one question came to my mind : I also had some

>studies that contributed multiple effect sizes in the meta-analysis that I

>recently ran thanks to Dr Viechtbauer's advice.

>For now I only used the rma function, but should I have used rma.mv because

>of these stuides that had multiple effect sizes ?

>

>Thank you !

>

>Norman

>

>________________________________________

>De: "James Pustejovsky" <jepusto at gmail.com>

>?: "Wolfgang Viechtbauer" <wolfgang.viechtbauer at maastrichtuniversity.nl>

>Cc: "r-sig-meta-analysis" <r-sig-meta-analysis at r-project.org>, "Huang Wu"

><huang.wu at wmich.edu>

>Envoy?: Mercredi 10 Juin 2020 05:08:09

>Objet: Re: [R-meta] weight in rmv metafor

>

>Hi Huang,

>

>I've written up some notes that add a bit of further intuition to the

>discussion that Wolfgang provided. The main case that I focus on is a model

>that is just a meta-analysis (i.e., no predictors) and that includes random

>effects to capture both between-study and within-study heterogeneity. I

>also say a little bit about meta-regression models with only study-level

>predictors.

>

>https://www.jepusto.com/weighting-in-multivariate-meta-analysis/

>

>Best,

>James

>

>On Sun, Jun 7, 2020 at 4:11 PM Viechtbauer, Wolfgang (SP) <

>wolfgang.viechtbauer at maastrichtuniversity.nl> wrote:

>

>> Of course the weights "impact the estimated fixed effects". But whether

>> studies with multiple effect sizes tend to receive more weight depends on

>> various factors, including the variances of the random effects and the

>> sampling error (co)variances.

>>

>> A more detailed discussion around the way weighting works in rma.mv

>> models can be found here:

>>

>> http://www.metafor-project.org/doku.php/tips:weights_in_rma.mv_models

>>

>> Note that weights(res, type="rowsum") currently only works in the 'devel'

>> version of metafor, so follow

>> https://wviechtb.github.io/metafor/#installation if you want to reproduce

>> this part as well.

>>

>> I hope this clarifies things.

>>

>> Best,

>> Wolfgang


        [[alternative HTML version deleted]]


_______________________________________________

R-sig-meta-analysis mailing list
R-sig-meta-analysis at r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis
--
Emerson M. Del Ponte
Universidade Federal de Vi?osa, Brazil
Chair of the Graduate Studies in Plant PathologyEIC for Tropical Plant Pathology Co-Founder of Open Plant PathologyMy websites: Twitter | GitHub | Google Scholar | ResearchGateTel +55 31 36124830



        [[alternative HTML version deleted]]




------------------------------

Subject: Digest Footer

_______________________________________________
R-sig-meta-analysis mailing list
R-sig-meta-analysis at r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis


------------------------------

End of R-sig-meta-analysis Digest, Vol 37, Issue 23
***************************************************