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

[James Pustejovsky]

There was a typo in the syntax I gave too. So just for the record, here is
the correct pseudo-syntax:

egger_multi <- rma.mv(yi = yi, V = sei^2, random = ~ 1 | studyID /
effectID, mods = ~ sei, data = dat)
coef_test(egger_multi, vcov = "CR2")

On Wed, Dec 9, 2020 at 1:38 PM Viechtbauer, Wolfgang (SP) <
wolfgang.viechtbauer using maastrichtuniversity.nl> wrote:

> Argh, of course :) Well, glad I could be of no assistance.
>
> Best,
> Wolfgang
>
> >-----Original Message-----
> >From: Dylan Johnson [mailto:dylanr.johnson using 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 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
>
> _______________________________________________
> R-sig-meta-analysis mailing list
> R-sig-meta-analysis using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis
>

	[[alternative HTML version deleted]]