Interaction overlaps {InteractionSet} | R Documentation

## Find overlaps between interactions in one or two dimensions

### Description

Find overlaps between interactions and linear intervals, between interactions and pairs of intervals, and between interactions and other interactions in a GInteractions or InteractionSet object.

### Usage

    ## S4 method for signature 'GInteractions,GInteractions'
    findOverlaps(query, subject, maxgap=-1L, minoverlap=0L,
        type=c("any", "start", "end", "within", "equal"),
        select=c("all", "first", "last", "arbitrary"),
        ignore.strand=TRUE, ..., use.region="both")

    ## S4 method for signature 'GInteractions,GInteractions'
    overlapsAny(query, subject, maxgap=-1L, minoverlap=0L,
        type=c("any", "start", "end", "within", "equal"),
        ignore.strand=TRUE, ..., use.region="both")

    ## S4 method for signature 'GInteractions,GInteractions'
    countOverlaps(query, subject, maxgap=-1L, minoverlap=0L,
        type=c("any", "start", "end", "within", "equal"),
        ignore.strand=TRUE, ..., use.region="both")

    ## S4 method for signature 'GInteractions,GInteractions'
    subsetByOverlaps(query, subject, maxgap=-1L, minoverlap=0L,
        type=c("any", "start", "end", "within", "equal"),
        ignore.strand=TRUE, ..., use.region="both")

### Arguments

`query, subject` | A Vector, GInteractions or InteractionSet object, depending on the specified method. At least one of these must be a GInteractions or InteractionSet object. Also, `subject` can be missing if `query` is a GInteractions or InteractionSet object.
`maxgap, minoverlap, type, select` | See `?findOverlaps` in the GenomicRanges package.
`ignore.strand` | A logical scalar indicating whether strand information in `query` or `subject` should be ignored. Note that the default setting here is different to that in `findOverlaps` as genomic interactions are usually unstranded.
`...` | Further arguments to pass to `findOverlaps`. This includes `drop.self` and `drop.redundant` when `subject` is missing.
`use.region` | A string specifying the regions to be used to identify overlaps.

### Value

For `findOverlaps`, a Hits object is returned if `select="all"`, and an integer vector of subject indices otherwise.

For `countOverlaps` and `overlapsAny`, an integer or logical vector is returned, respectively.

For `subsetByOverlaps`, a subsetted object of the same class as `query` is returned.

### Overview of overlaps for GInteractions

All methods can be applied using a GInteractions as either the `query` or `subject`, and a Vector as the other argument. In such cases, the Vector is assumed to represent some region on the linear genome (e.g., GRanges) or set of such regions (GRangesList). An overlap will be defined between the interval and an GInteractions interaction if either anchor region of the latter overlaps the former. This is considered to be a one-dimensional overlap, i.e., on the linear genome.

The same methods can be applied using two GInteractions objects as the `query` and `subject`. In such cases, a two-dimensional overlap will be computed between the anchor regions of the two objects. An overlap is defined if each anchor region of the first object overlaps at least one anchor region of the second object, and each anchor region of the second object overlaps at least one anchor region of the first object, i.e., there are overlapping areas in the two-dimensional interaction space. If `subject` is missing, overlaps will be computed between interactions in `query`.

### Description of overlap methods

When `select="all"`, `findOverlaps` returns a Hits object containing overlapping pairs of queries and subjects (or more specifically, their indices in the supplied objects - see `?findOverlaps` for more details). For other values of `select`, an integer vector is returned with one entry for each element of `query`, which specifies the index of the chosen (first, last or arbitrary) overlapping feature in `subject` for that query. Queries with no overlaps at all are assigned `NA` values.

For the other methods, `countOverlaps` returns an integer vector indicating the number of elements in `subject` that were overlapped by each element in `query`. `overlapsAny` returns a logical vector indicating which elements in `query` were overlapped by at least one element in `subject`. `subsetByOverlaps` returns a subsetted `query` containing only those elements overlapped by at least one element in `subject`.

### Choice of regions to define overlaps

For one-dimensional overlaps, `use.region="both"` by default such that overlaps with either anchor region are considered. If `use.region="first"`, overlaps are only considered between the interval and the first anchor region. Similarly, if `use.region="second"`, only the second anchor region is used.

Equivalent choices are available for two-dimensional overlaps:

- By default, `use.region="both"` such that the order of first/second anchor regions in the query and subject is ignored. This means that the first anchor region in the query can overlap both the first or second anchor regions in the subject. Similarly, the second anchor region in the query can overlap both the first or ssecond anchor regions in the subject.

- If `use.region="same"`, overlaps are only considered between the first anchor regions for the query and the subject, or between the second anchor regions for the query and subject. Overlaps between the first query region and the second subject region, or the second query region and the first subject region, are ignored.

- If `use.region="reverse"`, overlaps are only considered between the first anchor regions for the query and the second anchor regions for the subject, and vice versa. Overlaps between the first query/subject regions or between the second query/subject regions are ignored.

The latter two options tend only to be useful if the order of first/second regions is informative.

### Details for InteractionSet

Each method can also be applied with InteractionSet objects, and the behaviour is largely the same as that described for GInteractions objects. For a given InteractionSet object `x`, the corresponding method is called on the GInteractions object in the `interactions` slot of `x`. The return value is identical to that from calling the method on `interactions(x)`, except for `subsetByOverlaps` for InteractionSet queries (which returns a subsetted InteractionSet object, containing only those rows/interactions overlapping the `subject`).

### Author(s)

Aaron Lun

### See Also

`InteractionSet-class`, `findOverlaps`, `linkOverlaps`

### Examples

    example(GInteractions, echo=FALSE)

    # Making a larger object, for more overlaps.
    Np <- 100
    N <- length(regions(gi))
    all.anchor1 <- sample(N, Np, replace=TRUE)
    all.anchor2 <- sample(N, Np, replace=TRUE)
    gi <- GInteractions(all.anchor1, all.anchor2, regions(gi))

    # GRanges overlaps:
    of.interest <- resize(sample(regions(gi), 2), width=1, fix="center")
    findOverlaps(of.interest, gi)
    findOverlaps(gi, of.interest)
    findOverlaps(gi, of.interest, select="first")
    overlapsAny(gi, of.interest)
    overlapsAny(of.interest, gi)
    countOverlaps(gi, of.interest)
    countOverlaps(of.interest, gi)
    subsetByOverlaps(gi, of.interest)
    subsetByOverlaps(of.interest, gi)

    # GRangesList overlaps:
    pairing <- GRangesList(first=regions(gi)[1:3], second=regions(gi)[4:6],
        third=regions(gi)[7:10], fourth=regions(gi)[15:17])
    findOverlaps(pairing, gi)
    findOverlaps(gi, pairing)
    findOverlaps(gi, pairing, select="last")
    overlapsAny(gi, pairing)
    overlapsAny(pairing, gi)
    countOverlaps(gi, pairing)
    countOverlaps(pairing, gi)
    subsetByOverlaps(gi, pairing)
    subsetByOverlaps(pairing, gi)

    # GInteractions overlaps (split into two):
    first.half <- gi[1:(Np/2)]
    second.half <- gi[Np/2+1:(Np/2)]
    findOverlaps(first.half, second.half)
    findOverlaps(first.half, second.half, select="arbitrary")
    overlapsAny(first.half, second.half)
    countOverlaps(first.half, second.half)
    subsetByOverlaps(first.half, second.half)

    findOverlaps(gi)
    countOverlaps(gi)
    overlapsAny(gi) # trivial result

    #################
    # Same can be done for an InteractionSet object:

    Nlibs <- 4
    counts <- matrix(rpois(Nlibs*Np, lambda=10), ncol=Nlibs)
    colnames(counts) <- seq_len(Nlibs)
    iset <- InteractionSet(counts, gi)

    findOverlaps(of.interest, iset)
    findOverlaps(iset, pairing)
    findOverlaps(iset[1:(Np/2),], iset[Np/2+1:(Np/2),])

    # Obviously returns InteractionSet objects instead
    subsetByOverlaps(of.interest, iset)
    subsetByOverlaps(iset, pairing)
    subsetByOverlaps(iset[1:(Np/2),], iset[Np/2+1:(Np/2),])

    # Self-overlaps
    findOverlaps(iset)
    countOverlaps(iset)
    overlapsAny(iset) # trivial result

[Package InteractionSet version 1.18.1 Index]