BigWigSummary.R

# devtools::source_gist('8d889e75c6cc80a1b27c147fedf32982', filename='BigWigSummary.R')

BigWigSummary <- function(files, ranges) {
  library(rtracklayer)
  library(tidyverse)
  require(pbmcapply)
  
  out_list <- pbmclapply(files, function(file) {
    vc <- mean(import.bw(file, selection=BigWigSelection(ranges), as='RleList')[ranges])
    names(vc) <- NULL
    return(vc)
  }, mc.cores = parallel::detectCores()/2 )
  output <- tbl_df(data.frame(coordinates=gsub(':\\*', '', paste0(ranges)), out_list, stringsAsFactors=FALSE))
  return(output)
}