This function generates a heatmap showing percentage of overlapping peaks between peak files.
overlap_heatmap(
peaklist,
interact = TRUE,
draw_cellnote = TRUE,
fill_diag = NA,
verbose = TRUE
)
A list of peak files as GRanges object.
Files must be listed and named using list()
.
e.g. list("name1"=file1, "name2"=file2)
.
If not named, default file names will be assigned.
Default TRUE. By default heatmap is interactive. If FALSE, heatmap is static.
Draw the numeric values within each heatmap cell.
Fill the diagonal of the overlap matrix.
Print messages.
An interactive heatmap
### Load Data ###
data("encode_H3K27ac") # example peakfile GRanges object
data("CnT_H3K27ac") # example peakfile GRanges object
### Create Named List ###
peaklist <- list("encode"=encode_H3K27ac, "CnT"=CnT_H3K27ac)
### Run ###
my_heatmap <- overlap_heatmap(peaklist = peaklist)
#> --- Running overlap_heatmap() ---
#> Warning: fill_diag must =NULL due to a known bug in the package 'heatmaply'.
#> Genreating precision-recall matrix.
#> overlap_heatmap(): Done in 2.1s.