Prepare data to plot overlap between data.table of SNPs and cell-type-specific epigenomic peaks and coaccessibility data.

CORCES2020_prepare_scATAC_peak_overlap(
  merged_DT,
  FDR_filter = NULL,
  snp_filter = "Consensus_SNP==TRUE",
  add_cicero = TRUE,
  annotate_genes = TRUE,
  return_counts = TRUE,
  verbose = TRUE
)

Source

https://doi.org/10.1038/s41588-020-00721-x

merged_DT <- echodata::get_Nalls2019_merged() dat_melt <- echoannot:::CORCES2020_prepare_scATAC_peak_overlap( merged_DT = merged_DT)