library("Gviz")
gen <- "hg38"
chr <- "chr1"
start <- 40000
end <- 50000
featureDisplay <- "H3K27me3"
if(interactive()){
regulatoryEvidenceBiomartTrackSingle <- regulatoryEvidenceBiomart_ENSEMBL(gen,chr,start,end,featureDisplay)
plotTracks(regulatoryEvidenceBiomartTrackSingle, from = start, to = end)
} else {
data(regulatoryEvidenceBiomartTrackSingle)
plotTracks(regulatoryEvidenceBiomartTrackSingle, from = start, to = end)
}
#####
library("Gviz")
gen <- "hg38"
chr <- "chr1"
start <- 40000
end <- 100000
featureDisplay <- c("H3K27me3","H3K36me3")
if(interactive()){
regulatoryEvidenceBiomartTrackMultiple<-regulatoryEvidenceBiomart_ENSEMBL (gen,chr,start,end,featureDisplay)
plotTracks(regulatoryEvidenceBiomartTrackMultiple, from = start, to = end)
} else {
data(regulatoryEvidenceBiomartTrackMultiple)
plotTracks(regulatoryEvidenceBiomartTrackMultiple, from = start, to = end)
}
#####
library("Gviz")
gen <- "hg38"
chr <- "chr1"
start <- 50000
end <- 100000
featureDisplay <- "all"
if(interactive()){
regulatoryEvidenceBiomartTrackAll<-regulatoryEvidenceBiomart_ENSEMBL (gen,chr,start,end,featureDisplay)
plotTracks(regulatoryEvidenceBiomartTrackAll, from = start, to = end)
} else {
data(regulatoryEvidenceBiomartTrackAll)
plotTracks(regulatoryEvidenceBiomartTrackAll, from = start, to = end)
}
Run the code above in your browser using DataLab