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", sep = c(":", "-"))

Arguments

object.list	A list of Seurat 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.
sep	Separators to use to extract genomic ranges from object row names. To specify different separators for different objects, pass a list of length equal to the length of object.list.

Value

Returns a GRanges object

Examples

UnifyPeaks(object.list = list(atac_small, atac_small))
#> GRanges object with 100 ranges and 0 metadata columns:
#>         seqnames          ranges strand
#>            <Rle>       <IRanges>  <Rle>
#>     [1]     chr1   565107-565550      *
#>     [2]     chr1   569174-569639      *
#>     [3]     chr1   713460-714823      *
#>     [4]     chr1   752422-753038      *
#>     [5]     chr1   762106-763359      *
#>     ...      ...             ...    ...
#>    [96]     chr1 1553343-1553743      *
#>    [97]     chr1 1555436-1556258      *
#>    [98]     chr1 1557272-1558191      *
#>    [99]     chr1 1559561-1559642      *
#>   [100]     chr1 1562519-1567986      *
#>   -------
#>   seqinfo: 1 sequence from an unspecified genome; no seqlengths