Skip to content

[R-meta] GOSH plots for multi-level meta (rma.mv)

3 messages · Hellen Mirr, Wolfgang Viechtbauer

#
Hello,

Apologies if this has already been answered, as I could not find any
threads on it.
I was wondering whether it is possible to create a GOSH plot for a
multi-level meta analysis that is an rma.mv object, and how I would go
about that.

Thank you in advance for all the help!

Best wishes
Hellen
#
Dear Hellen,

This is currently not implemented in metafor. In principle, the idea of a GOSH plot does generalize to more complex models although one needs to think about whether one would want to create subsets based on the indiviual estimates or based on some higher-level grouping variable. For example, suppose we have a multilevel structure such as:

study  esid  yi vi
------------------
1      1     .  . 
1      2     .  .
2      1     .  .
3      1     .  .
3      2     .  .
3      3     .  .
4      1     .  .

So, what are, for example, then the subsets of size 2? Are they based just on the rows? Then the estimates in row 1 and 2 would be one such subset. Or does one base the subsets on the studies? Then rows 1, 2, 3 (i.e., studies 1 and 2) would be such a subset.

This could all be implemented (just like cooks.distance() and rstudent() allow for the optional specification of a clustering variable), but I haven't done this.

Aside from that: Fitting rma.mv models can take a bit of time. Doing this 1000's of times (based on all possible subsets) could take a LONG time.

Best,
Wolfgang
#
Dear Wolfgang,

Thank you very much for your clear explanation.

Best,
Hellen

On Fri, Oct 30, 2020 at 11:48 AM Viechtbauer, Wolfgang (SP) <
wolfgang.viechtbauer at maastrichtuniversity.nl> wrote: