[R-meta] Egger's test with multilevel meta analysis
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 at maastrichtuniversity.nl> wrote:
Argh, of course :) Well, glad I could be of no assistance. Best, Wolfgang
-----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
_______________________________________________ R-sig-meta-analysis mailing list R-sig-meta-analysis at r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-meta-analysis