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"