Prepare ASTARR Input MACS Peaks (Replicates) 02

Create metadata

Set environment

Code

suppressMessages(suppressWarnings(source("../run_config_project_sing.R")))
show_env()

You are working on        Singularity: singularity_proj_encode_fcc 
BASE DIRECTORY (FD_BASE): /data/reddylab/Kuei 
REPO DIRECTORY (FD_REPO): /data/reddylab/Kuei/repo 
WORK DIRECTORY (FD_WORK): /data/reddylab/Kuei/work 
DATA DIRECTORY (FD_DATA): /data/reddylab/Kuei/data 

You are working with      ENCODE FCC 
PATH OF PROJECT (FD_PRJ): /data/reddylab/Kuei/repo/Proj_ENCODE_FCC 
PROJECT RESULTS (FD_RES): /data/reddylab/Kuei/repo/Proj_ENCODE_FCC/results 
PROJECT SCRIPTS (FD_EXE): /data/reddylab/Kuei/repo/Proj_ENCODE_FCC/scripts 
PROJECT DATA    (FD_DAT): /data/reddylab/Kuei/repo/Proj_ENCODE_FCC/data 
PROJECT NOTE    (FD_NBK): /data/reddylab/Kuei/repo/Proj_ENCODE_FCC/notebooks 
PROJECT DOCS    (FD_DOC): /data/reddylab/Kuei/repo/Proj_ENCODE_FCC/docs 
PROJECT LOG     (FD_LOG): /data/reddylab/Kuei/repo/Proj_ENCODE_FCC/log 
PROJECT REF     (FD_REF): /data/reddylab/Kuei/repo/Proj_ENCODE_FCC/references

Set global variables

Code

TXT_FOLDER_REGION = "fcc_astarr_macs_narrowpeak"

Define file labeling

Code

### set directory
txt_folder = TXT_FOLDER_REGION
txt_fdiry  = file.path(FD_RES, "region", txt_folder)
txt_fglob  = file.path(txt_fdiry, "*narrowpeak*")

### get file names
vec_txt_fpath = sort(Sys.glob(txt_fglob))
vec_txt_fname = basename(vec_txt_fpath)
vec_txt_label = paste0("fcc_astarr_macs_input_rep", 1:6)

print(vec_txt_fname)

[1] "K562.hg38.ASTARR.macs.KS91.Input.rep1.narrowpeak.bed.gz"
[2] "K562.hg38.ASTARR.macs.KS91.Input.rep2.narrowpeak.bed.gz"
[3] "K562.hg38.ASTARR.macs.KS91.Input.rep3.narrowpeak.bed.gz"
[4] "K562.hg38.ASTARR.macs.KS91.Input.rep4.narrowpeak.bed.gz"
[5] "K562.hg38.ASTARR.macs.KS91.Input.rep5.narrowpeak.bed.gz"
[6] "K562.hg38.ASTARR.macs.KS91.Input.rep6.narrowpeak.bed.gz"

Code

### init info table
dat = data.frame(
    "Folder" = txt_folder,
    "FName"  = vec_txt_fname,
    "Label"  = vec_txt_label,
    "FPath"  = vec_txt_fpath
)

### assign and show
dat_region_label = dat
fun_display_table(dat)

Folder	FName	Label	FPath
fcc_astarr_macs_narrowpeak	K562.hg38.ASTARR.macs.KS91.Input.rep1.narrowpeak.bed.gz	fcc_astarr_macs_input_rep1	/data/reddylab/Kuei/repo/Proj_ENCODE_FCC/results/region/fcc_astarr_macs_narrowpeak/K562.hg38.ASTARR.macs.KS91.Input.rep1.narrowpeak.bed.gz
fcc_astarr_macs_narrowpeak	K562.hg38.ASTARR.macs.KS91.Input.rep2.narrowpeak.bed.gz	fcc_astarr_macs_input_rep2	/data/reddylab/Kuei/repo/Proj_ENCODE_FCC/results/region/fcc_astarr_macs_narrowpeak/K562.hg38.ASTARR.macs.KS91.Input.rep2.narrowpeak.bed.gz
fcc_astarr_macs_narrowpeak	K562.hg38.ASTARR.macs.KS91.Input.rep3.narrowpeak.bed.gz	fcc_astarr_macs_input_rep3	/data/reddylab/Kuei/repo/Proj_ENCODE_FCC/results/region/fcc_astarr_macs_narrowpeak/K562.hg38.ASTARR.macs.KS91.Input.rep3.narrowpeak.bed.gz
fcc_astarr_macs_narrowpeak	K562.hg38.ASTARR.macs.KS91.Input.rep4.narrowpeak.bed.gz	fcc_astarr_macs_input_rep4	/data/reddylab/Kuei/repo/Proj_ENCODE_FCC/results/region/fcc_astarr_macs_narrowpeak/K562.hg38.ASTARR.macs.KS91.Input.rep4.narrowpeak.bed.gz
fcc_astarr_macs_narrowpeak	K562.hg38.ASTARR.macs.KS91.Input.rep5.narrowpeak.bed.gz	fcc_astarr_macs_input_rep5	/data/reddylab/Kuei/repo/Proj_ENCODE_FCC/results/region/fcc_astarr_macs_narrowpeak/K562.hg38.ASTARR.macs.KS91.Input.rep5.narrowpeak.bed.gz
fcc_astarr_macs_narrowpeak	K562.hg38.ASTARR.macs.KS91.Input.rep6.narrowpeak.bed.gz	fcc_astarr_macs_input_rep6	/data/reddylab/Kuei/repo/Proj_ENCODE_FCC/results/region/fcc_astarr_macs_narrowpeak/K562.hg38.ASTARR.macs.KS91.Input.rep6.narrowpeak.bed.gz

Save results

Code

txt_folder = TXT_FOLDER_REGION
txt_fdiry  = file.path(FD_RES, "region", txt_folder, "summary")
txt_fname  = "metadata.label.tsv"
txt_fpath  = file.path(txt_fdiry, txt_fname)

dir.create(txt_fdiry, showWarnings = FALSE)
dat = dat_region_label
write_tsv(dat, txt_fpath)