Дисклеймер: кратко основные вещи по Bioconductor GenomicRanges
library("GenomicRanges")
### Сборка объекта GenomicRanges с нуля
GenRange <- GRanges(
seqnames = Rle(c("chrX", "chrY", "chrZ"), c(1, 1, 2)), # задание "хромосом" и числа блоков на каждой как Rle
ranges = IRanges(start = c(1, 3, 19, 30), end = c(7, 5, 35, 40), names = c('One', '
(
Read more... )