consensus_peaks: Consensus peaks

In anilchalisey/chompR: ChIP and ATAC-seq data processing pipeline

Description

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)

Usage

consensus_peaks(peaks, replicates, plot = FALSE)

Arguments

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.

Value

A GRanges object of consensus peaks

Examples

## Not run: 
conPeaks <- consensus_peaks(peaks = peaks, replicates = 2)

## End(Not run)

anilchalisey/chompR documentation built on May 9, 2019, 3:59 a.m.