data(exampleMethyGenoSet)
if (require(TxDb.Hsapiens.UCSC.hg19.knownGene) && require(Gviz)) {
## define data track
exampleMethyGenoSet <- checkChrName(exampleMethyGenoSet)
dTrack <- DataTrack(range=suppressWarnings(as(rowRanges(exampleMethyGenoSet), 'GRanges')), data=t(exprs(exampleMethyGenoSet)),
chromosome='chr21', type='heatmap')
## build annotation tracks
annotationTracks <- buildAnnotationTracks('1826', includeGeneBody = FALSE, genomicFeature = "TxDb.Hsapiens.UCSC.hg19.knownGene", cytobandInfo=NA, CpGInfo=NA)
trackList <- c(annotationTracks, list(dTrack))
plotTracksWithDataTrackInfo(trackList, labels=colnames(exampleMethyGenoSet), grange2show = attr(annotationTracks, 'grange2show'))
}
Run the code above in your browser using DataLab