pheatmap order group and set color

 

pheno_data$group = factor(pheno_data$group,level=unique(pheno_data$group))

expr_data = expr_data[, as.character(rownames(pheno_data))]

color = colorRampPalette(c("navy", "white", "firebrick3"))(100),
border=TRUE,
scale="row",
border_color ="black"