tests_not_run/test_genomic_regions_merge.R

library(GenomicRanges)
library(GetoptLong)

source("/home/guz/project/development/epik/R/genomic_region_merge.R")

gr = GRanges(seqnames = "chr1", ranges = IRanges(start = c(1, 4, 8 ,16), 
    end = c(2, 5, 10, 30)), value = 1:4)
reduce2(gr, gap = bp(2))
reduce2(gr, gap = 0.6)
reduce2(gr, gap = 0.6, max_gap = 4)
jokergoo/epik documentation built on Sept. 28, 2019, 9:20 a.m.