library("Gviz")
gen <- "hg38"
chr <- "chr1"
start <- 10000
end <- 50000
featureDisplay <- "CTCF enriched"
if(interactive()){
regulatorySegmentsBiomartTrackSingle<-regulatorySegmentsBiomart_ENSEMBL(gen,chr,start,end,featureDisplay)
plotTracks(regulatorySegmentsBiomartTrackSingle, from = start, to = end)
} else {
data(regulatorySegmentsBiomartTrackSingle)
plotTracks(regulatorySegmentsBiomartTrackSingle, from = start, to = end)
}
####
library("Gviz")
gen <- "hg38"
chr <- "chr1"
start <- 10000
end <- 50000
featureDisplay <- c("CTCF enriched","Predicted Promoter Flank")
if(interactive()){
regulatorySegmentsBiomartTrackMultiple<-regulatorySegmentsBiomart_ENSEMBL(gen,chr,start,end,featureDisplay)
plotTracks(regulatorySegmentsBiomartTrackMultiple, from = start, to = end)
} else {
data(regulatorySegmentsBiomartTrackMultiple)
plotTracks(regulatorySegmentsBiomartTrackMultiple, from = start, to = end)
}
####
library("Gviz")
gen <- "hg38"
chr <- "chr1"
start <- 10000
end <- 50000
featureDisplay <- "all"
if(interactive()){
regulatorySegmentsBiomartTrackAll<-regulatorySegmentsBiomart_ENSEMBL(gen,chr,start,end,featureDisplay)
plotTracks(regulatorySegmentsBiomartTrackAll, from = start, to = end)
} else {
data(regulatorySegmentsBiomartTrackAll)
plotTracks(regulatorySegmentsBiomartTrackAll, from = start, to = end)
}
Run the code above in your browser using DataLab