“…In this case, including that single discrete variable as a moderator in the syntax (i.e., after the mods~element) would suffice. However, it has become rather common to report on the mean effect (as well as significance and confidence interval) of all categories of a discrete potential moderating variable (see, for instance, Assink et al, 2015;Houben et al, 2015;Rapp, Van den Noortgate, Broekaert, & Vanderplasschen, 2014;Van der Hallen, Evers, Brewaeys, Van den Noortgate, & Wagemans, 2015;Van der Stouwe, Asscher, Stams, Dekovic, & Van der Laan, 2014;Weisz et al, 2013 # Determine the potential moderating effect of publication status; # Unpublished studies are now tested against published # studies, so published studies serve as the reference category; # Print the results stored in the object "published" on screen. published <-rma.mv(y, v, mods =~pstatnotpub, random = list(~1 | effectsizeID,1 | studyID), tdist=TRUE, data=dataset) summary(published, digits=3)…”