This function outputs a table summarizing information on the peak files. Provides the total number of peaks and the percentage of peaks in blacklisted regions.
peak_info(peaklist, blacklist)
A summary table of peak information
### Load Data ###
data("encode_H3K27ac") # example peakfile GRanges object
data("CnT_H3K27ac") # example peakfile GRanges object
data("hg19_blacklist") # example blacklist GRanges object
### Named Peaklist ###
peaklist <- list("encode"=encode_H3K27ac, "CnT"=CnT_H3K27ac)
### Run ###
df <- peak_info(peaklist = peaklist,
blacklist = hg19_blacklist)
#> --- Running peak_info() ---
#> Done.