[R-meta] Egger's test with multilevel meta analysis

-----Original Message-----
From: Dylan Johnson [mailto:dylanr.johnson at mail.utoronto.ca]
Sent: Wednesday, 09 December, 2020 20:26
To: Viechtbauer, Wolfgang (SP)
Cc: R meta
Subject: RE: [R-meta] Egger's test with multilevel meta analysis

The issue what that coeftest is a different function than the one I should
have been using (I.e. coef_test).

Thanks!

Dylan Johnson, MSc
MA?Student,?School and Clinical Child Psychology
Department of?Applied Psychology?and?Human Development
University of Toronto
252 Bloor Street West
Toronto, ON M5S 1V6

From: Viechtbauer, Wolfgang (SP)
Sent: December 9, 2020 2:23 PM
To: Dylan Johnson
Cc: R meta
Subject: RE: [R-meta] Egger's test with multilevel meta analysis

EXTERNAL EMAIL:

Try:

coeftest(egger_multi, vcov = "CR2", cluster = dataset$COHORT_ID)

If this doesn't help, what version of clunSandwich are you using? (see
sessionInfo()).

Best,
Wolfgang

-----Original Message-----
From: Dylan Johnson [mailto:dylanr.johnson at mail.utoronto.ca]
Sent: Wednesday, 09 December, 2020 20:01
To: Viechtbauer, Wolfgang (SP)
Cc: R meta
Subject: RE: [R-meta] Egger's test with multilevel meta analysis

Unfortunately the error persists even after making your correction

Dylan Johnson, MSc
MA Student, School and Clinical Child Psychology
Department of Applied Psychology and Human Development
University of Toronto
252 Bloor Street West
Toronto, ON M5S 1V6

From: Viechtbauer, Wolfgang (SP)
Sent: December 9, 2020 1:56 PM
To: Dylan Johnson
Cc: R meta
Subject: RE: [R-meta] Egger's test with multilevel meta analysis

That call to rma.mv() doesn't look right. Maybe you meant:

egger_multi <- rma.mv(HEDGE_G, HEDGE_VAR, random = ~ 1 |
COHORT_ID/EFFECT_ID, mods = ~ STD_ERR, data = dataset)

Best,
Wolfgang

-----Original Message-----
From: R-sig-meta-analysis [mailto:r-sig-meta-analysis-bounces at r-

project.org]

On Behalf Of Dylan Johnson
Sent: Wednesday, 09 December, 2020 19:37
To: James Pustejovsky; Tobias Saueressig
Cc: R meta
Subject: Re: [R-meta] Egger's test with multilevel meta analysis

I have tried the following:

egger_multi <- rma.mv(HEDGE_G, HEDGE_VAR, random = ~ 1 | COHORT_ID,
EFFECT_ID, mods = ~ STD_ERR, data = dataset)

coeftest(egger_multi, vcov = "CR2")

When I run the coeftest I receive the error:

Error in diag(se) : invalid 'nrow' value (too large or NA)
In addition: Warning message:
In diag(se) : NAs introduced by coercion

Dylan Johnson, MSc

MA Student, School and Clinical Child Psychology
Department of Applied Psychology and Human Development

University of Toronto
252 Bloor Street West

Toronto, ON M5S 1V6