How do I extract gene names from different portions of my heatmap?
1
0
Entering edit mode
2.3 years ago
mgranada3 ▴ 50

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.

enter image description here

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)
pheatmap heatmap deseq2 • 1.0k views
ADD COMMENT
2
Entering edit mode
ADD COMMENT
0
Entering edit mode

As detailed in the linked reply above, you need to (a) decide at which level to cut the dendrogram and then (b) obtain the names of the members for the nodes. DAve Tang explained the relationship between (p)heatmaps and dendrograms well, if you need a primer on that.

ADD REPLY

Login before adding your answer.

Traffic: 2722 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6