I have this heatmap that I generated using DESeq2 and pheatmap(). I had it perform hierarchal clustering on the rows. Now, how do I identify which genes are in the primarily blue band (on top) or the red band (on bottom). I am fairly new to R, so I have been using the code displayed on the vignette. Code below.
library("pheatmap") df <- as.data.frame(colData(dds)[,c("condition","type")]) pheatmap(assay(vsd), cluster_rows=TRUE, show_rownames=FALSE, cluster_cols=FALSE, annotation_col=df)