File:ComplexHeatmap2.png

Summary

# Simulate data
library(ComplexHeatmap)
ng <- 30; ns <- 20
set.seed(1)
mat <- matrix(rnorm(ng*ns), nr=ng, nc=ns)
colnames(mat) <- 1:ns
rownames(mat) <- 1:ng

# color bar on RHS
ind_e <- 1:round(ng/3)
ind_m <- (1+round(ng/3)):ng
epimes <- rep(c("epi", "mes"), c(length(ind_e), length(ind_m)))
row_ha <- rowAnnotation(epimes = epimes,
        col = list(epimes = c("epi" = "orange", "mes" = "darkgreen")))

# color bar on Top
tumortype <- rep(c("carcinoma", "sarcoma"), c(round(ns/2), (ns-round(ns/2))))
bin <- sample(c("resistant", "sensitive"), ns, replace = T)
top_ha <- HeatmapAnnotation(bin = bin,
                        col = list(bin =c("resistant" = "violet", 
                                          "sensitive" = "mediumseagreen")))

# combine
Heatmap(mat, 
        name = "expression",
        column_split = tumortype,
        row_split = epimes,
        top_annotation = top_ha,   # color bat on top
        right_annotation = row_ha, # color bar on RHS
        show_row_names = TRUE,     # FALSE if ng is large
        column_names_rot = 60)