library(GenomicRanges)
data(cage)
data(cpgi)
data(promoters)
cage$tpm = NULL
targets = GRangesList(cage=cage, cpgi=cpgi)
sml = ScoreMatrixList(targets, promoters, bin.num=10)
kmeans.clust = kmeans(sml$cage,3)
sml.ordered = orderBy(sml, kmeans.clust$cluster)
multiHeatMatrix(sml.ordered)
Run the code above in your browser using DataLab