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

[Dylan Johnson]

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)<mailto:wolfgang.viechtbauer using maastrichtuniversity.nl>
Sent: December 9, 2020 2:23 PM
To: Dylan Johnson<mailto:dylanr.johnson using mail.utoronto.ca>
Cc: R meta<mailto:r-sig-meta-analysis using r-project.org>
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 using 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 using 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

	[[alternative HTML version deleted]]