getConsensusPeaks: Given a GRangeslist object with peaks for each samples,...

In Mathelab/ALTRE: Workflow for Identifying Regulatory Elements from Chromatin Accessibility Assay Data

Description

Given a GRangeslist object with peaks for each samples, determine the consensus peaks (found in at least N replicates, where N is input by the user) for each sample type

Usage

getConsensusPeaks(samplepeaks, minreps)

Arguments

samplepeaks	A GRangesList object comprising one GRanges object (peaks) for each sample (output of loadBEDFiles() function)
minreps	minimum number of replicate samples that a peak should be contained in to be called as a consensus peak. This cutoff will be applied to both samples.

Value

a list comprising: 1) a GRangeslist with one GRange for each sample type which contains consensus peaks 2) a summary statistic table

Examples

## Not run: 
csvfile <- loadCSVFile("DNaseEncodeExample.csv")
samplePeaks <- loadBedFiles(csvfile)
consensusPeaks <- getConsensusPeaks(samplepeaks = samplePeaks, minreps = 2)

## End(Not run)

Mathelab/ALTRE documentation built on May 7, 2019, 3:41 p.m.