intersect_tads: Finds intersection between 2 sets of TADs.

In rz6/DIADEM: Hi-C differential analysis via dependency modelling of chromatin interactions with generalized linear models

Description

Intersection is defined for the same chromosomes, for example for 2 TAD sets from chromosome 18 if there is TAD in set1 with coordinates 100, 120 (start, end respectively) and in set2 there is a TAD with coordinates 110, 140 then their intersection is interval with coordinates: 110, 120.

Usage

intersect_tads(tads1, tads2)

Arguments

tads1	data frame of TAD set 1 with columns: start, end
tads2	data frame of TAD set 2 with columns: start, end

Value

data frame with start, end columns containing TAD intersection intervals

See Also

Intervals, interval_intersection for functions used to find intersection between 2 sets of intervals

Examples

# simple artificial data set of TADs
tads1 <- data.frame(start = c(1,10,30), end = c(5,15,35))
tads2 <- data.frame(start = c(3,14,28), end = c(12,18,40))
tads.intersection <- intersect_tads(tads1, tads2)
print(tads.intersection)

rz6/DIADEM documentation built on Dec. 31, 2019, 3:51 a.m.