RPM 2 BW

Script

This is an R function that can turn an RPM coverage object (rle-list) into a BigWig file that can be loaded into igv for track viewing or other purposes.

library(rtracklayer)
library(base)
RDatafiles <- dir(pattern="rpm.RData$")

for (i in 1:length(RDatafiles)) {
    cat("Loading ", RDatafiles[i], "\n")
    load(RDatafiles[i])
}

rpmfiles <- ls(pattern = "rpm$")
label <- paste0(gsub("_150_rpm$", "", rpmfiles), ".bw")
rpmfile.list <- mget(rpmfiles)

for (i in 1:length(rpmfiles)) {
    cat("converting ", rpmfiles[i], "to bigwig", "\n")
    export.bw(object=rpmfile.list[i][[1]], con=label[i])
}