[R-meta] Egger-type test for multi-level meta-analysis

[James Pustejovsky]

Correct. The slope on sqrt(vi) indicates an association between SE and
effect size. The intercept is a "PET"-style estimate of the average effect
size in a population of infinitely large studies (i.e. SE = 0).

James

On Mon, Jun 19, 2023 at 9:21 AM Dr. Guido Schwarzer via R-sig-meta-analysis
<r-sig-meta-analysis using r-project.org> wrote:

> Hi all,
>
> Another question on multi-level models (while I am still waiting for an
> answer on my previous one ;-) ).
>
> I would like to conduct a test for small-study effects for data from a
> three-level model, e.g., for the dataset dat.konstantopoulos2011.
>
> library("metafor")
> library("metadat")
> m.ml <- rma.mv(yi, vi, random = ~ 1 | district / school, data =
> dat.konstantopoulos2011)
>
> If I understand James' comment from February 2018 correctly (
> https://stat.ethz.ch/pipermail/r-sig-meta-analysis/2018-February/000610.html),
> I could conduct an Egger-type test for small study effects by using
> cluster-robust variance estimation following a multi-level meta-regression
> with the standard error as moderator:
>
> sse.ml <- update(m.ml, mods = sqrt(vi))
> library("clubSandwich")
> conf_int(sse.ml, vcov = "CR2")
>
> Did I get this right?
>
> Best wishes,
> Guido
>
>
> _______________________________________________
> R-sig-meta-analysis mailing list @ R-sig-meta-analysis using r-project.org
> To manage your subscription to this mailing list, go to:
> https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis
>

	[[alternative HTML version deleted]]