Unify genomic ranges

Create a unified set of non-overlapping genomic ranges from multiple Seurat objects containing single-cell chromatin data.

UnifyPeaks(object.list, mode = "reduce")

Arguments

object.list

A list of Seurat objects or ChromatinAssay objects

mode

Function to use when combining genomic ranges. Can be "reduce" (default) or "disjoin". See reduce and disjoin for more information on these functions.

Value

Returns a GRanges object

Examples

UnifyPeaks(object.list = list(atac_small, atac_small))
#> GRanges object with 323 ranges and 0 metadata columns:
#>         seqnames          ranges strand
#>            <Rle>       <IRanges>  <Rle>
#>     [1]     chr1   713460-714823      *
#>     [2]     chr1   752422-753038      *
#>     [3]     chr1   762106-763359      *
#>     [4]     chr1   779589-780271      *
#>     [5]     chr1   804872-805761      *
#>     ...      ...             ...    ...
#>   [319]     chr1 9299648-9300348      *
#>   [320]     chr1 9327071-9327557      *
#>   [321]     chr1 9335457-9336176      *
#>   [322]     chr1 9349019-9350779      *
#>   [323]     chr1 9352328-9354391      *
#>   -------
#>   seqinfo: 1 sequence from an unspecified genome; no seqlengths