Given a GRangesList with peaks for each sample, this function determines consensus peaks (found in at least N replicates, where N is defined by the user)
consensus_peaks(peaks, replicates, plot = FALSE)
| peaks | A GRangesList object consisting of one GRanges object for each set of peaks. |
|---|---|
| replicates | The minimum number of replicates in which a peak must be present for it to count as a consensus peak. |
| plot | Boolean to indicate whether or not to generate diagnostic plots. Default = FALSE. |
A GRanges object of consensus peaks
# NOT RUN { conPeaks <- consensus_peaks(peaks = peaks, replicates = 2) # }