R/meta_group_func.R

In setgree/ResultsStandardizeR: Some functions that Seth Green & John-Henry Pezzuto wrote to assist with a forthcoming meta-analysis

# meta_group_func

#' # subgroup function
meta_group_func <- function(dataset) {
  print(paste0('number of studies: ', 
               length(unique(dataset$unique_study_id))))
  print(paste0('number of papers: ',
               length(unique(dataset$unique_paper_id))))
  
  # relationship between d and se_d? 
  print('relationship between D and standard error of D:')
  print(summary(lm(d ~ st_err_d, data = dataset)))
  
  print('robustness model: naive meta-model with clustering at paper level') 
  print(metafor::robust(x = rma(yi = d, vi = var_d, data = dataset), 
                        cluster = dataset$unique_paper_id))
  # funnel plot
  metafor::funnel(metafor::rma.mv(yi = d, V = var_d, 
                                  random = ~ 1 | unique_paper_id, 
                                  data = dataset))
}