library("Gviz")
chr <- "chr1"
start <- 236728
end <- 238778
gen="hg19"
extdata <- system.file("extdata", package="coMET",mustWork=TRUE)
bedFilePath <- file.path(extdata, "RoadMap/CD3-DS17198.hg19.bed")
if(interactive()){
dgfootprints_RoadMapSingle <- dgfootprints_RoadMap(gen,chr,start, end, bedFilePath, tissueGroupDisplay='Blood & T-cell' )
plotTracks(dgfootprints_RoadMapSingle, from = start, to = end)
} else {
data(dgfootprints_RoadMapSingle)
plotTracks(dgfootprints_RoadMapSingle, from = start, to = end)
}
Run the code above in your browser using DataLab