data("HOT.spots")
temp <- tempfile() url <- "http://metatracks.encodenets.gersteinlab.org" download.file(file.path(url, "HOT_All_merged.tar.gz"), temp) temp2 <- tempfile() download.file(file.path(url, "HOT_intergenic_All_merged.tar.gz"), temp2) untar(temp, exdir=dirname(temp)) untar(temp2, exdir=dirname(temp)) f <- dir(dirname(temp), "bed$") HOT.spots <- sapply(file.path(dirname(temp), f), toGRanges, format="BED") names(HOT.spots) <- gsub("_merged.bed", "", f) HOT.spots <- sapply(HOT.spots, unname) HOT.spots <- GRangesList(HOT.spots) save(list="HOT.spots", file="data/HOT.spots.rda", compress="xz", compression_level=9)
data(HOT.spots)
elementNROWS(HOT.spots)
Run the code above in your browser using DataLab