Extract a block from bigWig files

extract.block.R

extract.block <- function(files, chr, start, end, verbose = TRUE){
		rl <- IRanges::RangesList(IRanges::IRanges(start=start, end=end))
		names(rl) <- chr
		cat("[extract.block] Reading in the files \n")
		rles <- lapply(files, function(xx) {
			#if (verbose){
			#	cat("Reading", xx, "\n")
			#}
	    	import(xx, as = "Rle", format = "bw", selection = BigWigSelection(rl))
		})
		cat("[extract.block] Building the block \n")
		megaMatrix <- do.call(cbind, lapply(rles, function(xx) as.numeric(xx[[chr]][start:end])))
		megaMatrix
	}