This is a test version of Biostars. For the public version, visit https://www.biostars.org.
What package to use in bioconductor to draw heatmap like this

Hi all, I'm new in this field, and I tried to cluster some TCGA cancer patients for further survival analysis. Just as the TCGA nature paper, I used the k-means method to cluster the normalized RNA-seq data. After that, how can I draw a heatmap like this, to cluster the samples, and do the further survival analysis? Is this done by some specific package? I think the heatmap.2 and pheatmap might not be able to draw the heatmap like this.

Thank you so much!

image: heatmap

The sample heatmap link:

http://www.nature.com/nature/journal/v497/n7447/fig_tab/nature12113_F3.html

image: heatmap from paper

rna-seq
• 6,400 views • link
updated by Ram Ram • written by Xinsen Xu • Xinsen Xu

It would be useful to see how "this" is like without having to search for the paper. Edit your post to provide a link, please.

• 0 views • link
written by h.mon h.mon

So sorry that it's my first time to post here. I tried to upload an image of the heatmap, however failed.

Now I've added the link. Please check. Thanks.

• 0 views • link
updated by Ram Ram • written by Xinsen Xu • Xinsen Xu

2 answers

Well, the first step is going to be determining what genes to include. This will likely involve a hypothesis test (differential expression). After that, heatmap.2, pheatmap, or even ComplexHeatmap packages will allow you to portray the data. Note that it is quite possible that your samples/genes will not cluster as nicely as the figure you give in your post.

• 0 views • link
updated by Ram Ram • written by Sean Davis Sean Davis

Thank you so much Sean!

I wanna to cluster the cancer patients into 3 or 4 different groups, with obvious different overall survival of each group. I just try to imitate the nature paper as the link:

http://www.nature.com/nature/journal/v497/n7447/fig_tab/nature12113_F3.html

Ideally, I wanna the heatmap to show the obvious different patient groups, with the highlight of the characteristic genes in each group, just as the four red rectangle in the nature heatmap.

Could you give me some ideas about that? Thanks.

• 0 views • link
updated by Ram Ram • written by Xinsen Xu • Xinsen Xu

heatmap.2 would be my guess too. they probably scaled/center the data before plotting the image (you can use the scale() function in R on both rows and cols). for the "bars" on the side and on top they could have done by using the hclust() function with the RColorBrewer package to assign a color to each cluster and passing then the data to the heatmap.2() function. I'm not sure about the white lines though.

• 0 views • link
updated by Ram Ram • written by TriS TriS

Log in to answer this question.

More posts like this