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 GenomicRanges::reduce() and GenomicRanges::disjoin() for more information on these functions.

Value

Returns a GRanges object

Examples

UnifyPeaks(object.list = list(atac_small, atac_small))
#> Warning: `UnifyPeaks()` was deprecated in Signac 1.17.0.
#> 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