rm(list=ls())

# meta-analysis code example for metafor
library(metafor)
library(lme4)

# citing metafor package
citation(package='metafor')
packageVersion('metafor') # to know current version number

df <- read.table('test.csv', header = TRUE, sep = ',')
res_IG_Y <- rma.glmm(measure="PLO", xi=num, ni=denom, data=df, slab=UniqueID)

### fit random-effects model in the two subgroups
res_IG_Y_frameless <- rma.glmm(measure="PLO", xi=num, ni=denom, data=df, slab=UniqueID, subset=(igns=="Y" & frame_type=='frameless')) # zero
res_IG_Y_frame <- rma.glmm(measure="PLO", xi=num, ni=denom, data=df, slab=UniqueID, subset=(igns=="Y" & frame_type=="frame"))