# Bartl et al. (submitted) metaregression script # install and load metafor package (Viechtbauer, 2010) install.packages("metafor") library(metafor) # load data (source: researchgate.net due to problems loading from OSF link) moddata<-read.csv("https://www.researchgate.net/profile/Gergely_Bartl/publication/335464885_metaregdata/data/5d678d65299bf11adf29ddec/metaregdata.csv") # run model with no moderators nomods<-rma(yi=TE, sei=seTE, data=moddata) nomods forest(nomods) # run model with moderators: current density and blinding # double-blind: 1, single-blind: 0 mods<-rma(yi=TE, sei=seTE, data=moddata, mods=cbind(dblind, density)) forest(mods) mods # Egger's regression test for funnel plot assymetry regtest(mods) # calculating lower and upper bound tau squared for both models # 95% CI, using 1.96SE # model with no moderators nomodstausqlb<-0.1421-1.96*0.0631 nomodstausqub<-0.1421+1.96*0.0631 nomodstausqlb nomodstausqub # model with moderators modstausqlb<-0.1107-1.96*0.0585 modstausqub<-0.1107+1.96*0.0585 modstausqlb modstausqub