Example: loading data from a bigWig file

load-bigwig.r

library(IRanges)
library(GenomicRanges)
library(rtracklayer)

# Select a BigWig file
bw.dir <- "/home/chromatin/roadmap/DNase_hypersensitivity/brain_fetal"
bw.file <- dir(bw.dir, full.names = TRUE, pattern = "*.bigWig")[1]

# Specify a genomic range
selection <- GRanges(seqnames = "chr4",
  ranges = IRanges(start =  46251581, end = 46392421))

# Extract BigWig data
bw.data <- import.bw(con = BigWigFile(bw.file),
  selection = BigWigSelection(selection))