R/get_CORCES2020_bulkATACseq_peaks.R
get_CORCES2020_bulkATACseq_peaks.RdEach row represents an individual peak identified in the bulk ATAC-seq data.
get_CORCES2020_bulkATACseq_peaks()Data originally from Corces et al. (bioRxiv) (doi:10.1038/s41588-020-00721-x ), as of May 2020. Specifically: STable2_Features_bulkATAC-seq_Peaks
Other CORCES2020:
CORCES2020_get_ATAC_peak_overlap(),
CORCES2020_get_hichip_fithichip_overlap(),
CORCES2020_prepare_bulkATAC_peak_overlap(),
CORCES2020_prepare_scATAC_peak_overlap(),
get_CORCES2020_cicero_coaccessibility(),
get_CORCES2020_hichip_fithichip_loop_calls(),
get_CORCES2020_scATACseq_celltype_peaks(),
get_CORCES2020_scATACseq_peaks()
if (FALSE) { # \dontrun{
dat <- readxl::read_excel(
file.path(
"~/Desktop/Fine_Mapping/echolocatoR",
"annotations/Coceres_2020",
"STable2_Features_bulkATAC-seq_Peaks.xlsx"
),
skip = 18
)
CORCES2020_bulkATACseq_peaks <- data.table::data.table(dat)
#### piggyback ####
tmp <- file.path(tempdir(), "CORCES2020_bulkATACseq_peaks.tsv.gz")
data.table::fwrite(CORCES2020_bulkATACseq_peaks, tmp, sep = "\t")
piggyback::pb_upload(
file = tmp,
repo = "RajLabMSSM/echoannot"
)
} # }