library("Gviz")
gen <- "hg19"
chr <- "chr7"
start <- 38290160
end <- 38303219
if(interactive()){
BROWSER.SESSION="UCSC"
mySession <- browserSession(BROWSER.SESSION)
genome(mySession) <- gen
track.name="Broad ChromHMM"
tablestrack<-tableNames(ucscTableQuery(mySession, track=track.name))
table.name<-tablestrack[1]
dnasetrack<-DNAseUCSC(gen,chr,start,end,mySession)
plotTracks(dnasetrack, from = start, to =end)
}else {
data(dnasetrack)
plotTracks(dnasetrack, from = start, to =end)
}
Run the code above in your browser using DataLab