/
ranges-pairs.R
118 lines (108 loc) · 3.99 KB
/
ranges-pairs.R
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
#' Pair together two ranges objects
#'
#' @param x,y Ranges objects to pair together.
#' @param maxgap,minoverlap The maximimum gap between intervals as an integer
#' greater than or equal to negative one. The minimum amount of overlap between intervals
#' as an integer greater than zero, accounting for the maximum gap.
#' @param suffix A character vector of length two used to identify metadata columns
#' coming from x and y.
#'
#' @details These functions return a DataFrame object, and is one way of
#' representing paired alignments with plyranges.
#'
#' @return a DataFrame with two ranges columns and the corresponding metadata
#' columns.
#'
#' @seealso [join_nearest()][join_overlap_inner()][join_precede()][join_follow()]
#' @examples
#' query <- data.frame(start = c(5,10, 15,20), width = 5, gc = runif(4)) %>%
#' as_iranges()
#' subject <- data.frame(start = 2:6, width = 3:7, label = letters[1:5]) %>%
#' as_iranges()
#'
#' pair_overlaps(query, subject)
#' pair_overlaps(query, subject, minoverlap = 5)
#' pair_nearest(query, subject)
#'
#'
#' query <- data.frame(seqnames = "chr1",
#' start = c(11,101),
#' end = c(21, 200),
#' name = c("a1", "a2"),
#' strand = c("+", "-"),
#' score = c(1,2)) %>%
#' as_granges()
#' subject <- data.frame(seqnames = "chr1",
#' strand = c("+", "-", "+", "-"),
#' start = c(21,91,101,201),
#' end = c(30,101,110,210),
#' name = paste0("b", 1:4),
#' score = 1:4) %>%
#' as_granges()
#'
#' # ignores strandedness
#' pair_overlaps(query, subject, suffix = c(".query", ".subject"))
#' pair_follow(query, subject, suffix = c(".query", ".subject"))
#' pair_precede(query, subject, suffix = c(".query", ".subject"))
#' pair_precede(query, subject, suffix = c(".query", ".subject"))
#' @rdname ranges-pairs
#' @export
pair_overlaps <- function(x, y, maxgap, minoverlap, suffix) {
UseMethod("pair_overlaps")
}
#' @export
pair_overlaps.IntegerRanges <- function(x, y, maxgap = -1L, minoverlap = 0L, suffix = c(".x", ".y")) {
hits <- findOverlaps(x,y, maxgap, minoverlap, type = "any", select = "all")
expand_by_hits(x,y, suffix, hits, return_data_frame = TRUE)
}
#' @export
pair_overlaps.GenomicRanges <- function(x, y, maxgap = -1L, minoverlap = 0L, suffix = c(".x", ".y")) {
hits <- findOverlaps(x,y, maxgap, minoverlap,
type = "any", select = "all", ignore.strand = TRUE)
expand_by_hits(x,y, suffix, hits, return_data_frame = TRUE)
}
#' @rdname ranges-pairs
#' @export
pair_nearest <- function(x, y, suffix) {
UseMethod("pair_nearest")
}
#' @export
pair_nearest.IntegerRanges <- function(x,y, suffix = c(".x", ".y")) {
hits <- nearest(x,y, select = "arbitrary")
expand_by_hits(x,y, suffix, hits, return_data_frame = TRUE)
}
#' @export
pair_nearest.GenomicRanges <- function(x,y, suffix = c(".x", ".y")) {
hits <- nearest(x,y, select = "arbitrary", ignore.strand = TRUE)
expand_by_hits(x,y, suffix, hits, return_data_frame = TRUE)
}
#' @rdname ranges-pairs
#' @export
pair_precede <- function(x, y, suffix) {
UseMethod("pair_precede")
}
#' @export
pair_precede.IntegerRanges <- function(x,y, suffix = c(".x", ".y")) {
hits <- precede(x,y)
expand_by_hits(x,y,suffix, hits, return_data_frame = TRUE)
}
#' @export
pair_precede.GenomicRanges <- function(x,y, suffix = c(".x", ".y")) {
hits <- precede(x,y, ignore.strand = TRUE)
expand_by_hits(x,y, suffix, hits, return_data_frame = TRUE)
}
#' @rdname ranges-pairs
#' @export
pair_follow <- function(x, y, suffix) {
UseMethod("pair_follow")
}
#' @export
pair_follow.IntegerRanges <- function(x,y, suffix = c(".x", ".y")) {
hits <- follow(x,y)
expand_by_hits(x,y,suffix, hits, return_data_frame = TRUE)
}
#' @export
pair_follow.GenomicRanges <- function(x,y, suffix = c(".x", ".y")) {
hits <- follow(x,y, ignore.strand = TRUE)
expand_by_hits(x,y, suffix, hits, return_data_frame = TRUE)
}