plot differentially accessible peaks for a given comparison
Arguments
- results
cinaR result object
- comparison
these are created by cinaR from `contrasts` user provided. If not selected the first comparison will be shown!
- ...
additional arguments for heatmap function, for more info `?pheatmap`
Examples
# \donttest{
library(cinaR)
data(atac_seq_consensus_bm) # calls 'bed'
# a vector for comparing the examples
contrasts <- sapply(strsplit(colnames(bed), split = "-", fixed = TRUE),
function(x){x[1]})[4:25]
results <- cinaR(bed, contrasts, reference.genome = "mm10")
#> >> Experiment type: ATAC-Seq
#> >> Matrix is filtered!
#> >> preparing features information... 2024-05-22 10:40:49
#> >> identifying nearest features... 2024-05-22 10:40:49
#> >> calculating distance from peak to TSS... 2024-05-22 10:40:49
#> >> assigning genomic annotation... 2024-05-22 10:40:49
#> >> assigning chromosome lengths 2024-05-22 10:40:52
#> >> done... 2024-05-22 10:40:52
#> >> Method: edgeR
#> FDR:0.05& abs(logFC)<0
#> >> Estimating dispersion...
#> >> Fitting GLM...
#> >> DA peaks are found!
#> >> No `geneset` is specified so immune modules (Chaussabel, 2008) will be used!
#> >> enrichment.method` is not selected. Hyper-geometric p-value (HPEA) will be used!
#> >> Mice gene symbols are converted to human symbols!
#> >> Enrichment results are ready...
#> >> Done!
heatmap_differential(results)
#> Warning: 'comparison' is not set so 'B6_NZO' will be used!
# }