scripts/proc_all_demographs.R

# loads processed files, adds demographics, spits out a bunch of plots

rm(list = ls())

BASE = "~/Desktop"
FIRST_BLOCK = "first_block"
setwd(BASE)

library(aceR)

# define for plots
factors = c("study_name", "age_group", "cohort_group")
variable_patterns = c("rt_", "acc_", "rw_")
variable_skip = c("cohort", "count", "length")

# load & process data by module
all_data_path = "ace_process"
out_path = "ace_process_demographs"
out_path_first_block = paste(out_path, FIRST_BLOCK, sep = "/")
out_path_first_block_plots = paste(out_path_first_block, "plots", sep = "/")

# make ouput directories
make_directory(out_path)
make_directory(out_path_first_block)
# make_directory(out_path_first_block_plots)

all_files = list.files(path = all_data_path, recursive = TRUE)
mods = unique(basename(all_files))
for (module in mods) {
  print(paste("processing", module, sep = " "))
  # load data
  mod_dat = load_files(path = all_data_path, pattern = module, recursive = TRUE)
  mod_name = gsub(".csv", "", module)
  mod_dat$pid = as.character(mod_dat$pid)
  mod_demo = mod_dat[order(mod_dat$pid), ]
  row.names(mod_demo) = NULL
  
  # remove rows with out study_name
  mod_demo = mod_demo[!is.na(mod_demo$study_name),]
  
  # output by tasks - nothing subsetted here
  mod_file_name = paste(out_path, module, sep ="/")
  write.csv(mod_demo, mod_file_name, na = "")
  
  # only look at first block

  mod_demo_first_block = aceR:::subset_first_block(mod_demo)
  mod_demo_first_block_file_name = paste(out_path, "first_block", module, sep = "/")
  write.csv(mod_demo_first_block, mod_demo_first_block_file_name, na = "")
  
  # # make box plots
  # mod_demo_first_block_columns = names(mod_demo_first_block)
  # mod_demo_first_block_variables = aceR:::multi_filter_vec(mod_demo_first_block_columns, variable_patterns)
  # mod_demo_first_block_variables = aceR:::multi_filter_out_vec(mod_demo_first_block_variables, variable_skip)
  # for (group in factors) {
  #   title_prefix = mod_name
  #   file_prefix = paste(mod_name, group, sep = "_")
  #   box_plots = make_box_plot_bulk(mod_demo_first_block, group, mod_demo_first_block_variables, title_prefix, file_prefix, out_path_first_block_plots)
  #   for (boxplot in box_plots) {
  #     plot(boxplot)
  #   }
  # }

}
gazzlab/AceScripts documentation built on May 16, 2019, 10:11 p.m.