-----Original Message-----
From: R-sig-meta-analysis [mailto:r-sig-meta-analysis-bounces at r-project.org] On
Behalf Of James Pustejovsky via R-sig-meta-analysis
Sent: Monday, 19 June, 2023 18:23
To: R Special Interest Group for Meta-Analysis
Cc: James Pustejovsky
Subject: Re: [R-meta] Egger-type test for multi-level meta-analysis
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 at 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