Question: What package to use in bioconductor to draw heatmap like this
3
gravatar for Xinsen Xu
3.2 years ago by
Xinsen Xu30
US/Boston/Harvard
Xinsen Xu30 wrote:

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!

The sample heatmap link:

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

rna-seq • 2.6k views
ADD COMMENTlink modified 3.2 years ago by TriS3.4k • written 3.2 years ago by Xinsen Xu30
4

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.

ADD REPLYlink written 3.2 years ago by h.mon18k

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.

ADD REPLYlink written 3.2 years ago by Xinsen Xu30
1
gravatar for Sean Davis
3.2 years ago by
Sean Davis24k
National Institutes of Health, Bethesda, MD
Sean Davis24k wrote:

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.  

ADD COMMENTlink modified 3.2 years ago • written 3.2 years ago by Sean Davis24k

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 hightlight 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.

ADD REPLYlink written 3.2 years ago by Xinsen Xu30
1

You can find a lot of tutorials on mkaing heatmaps via a search engine, e.g.: http://sebastianraschka.com/Articles/heatmaps_in_r.html

ADD REPLYlink written 3.2 years ago by Alex Reynolds25k
1
gravatar for TriS
3.2 years ago by
TriS3.4k
United States, Buffalo
TriS3.4k wrote:

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 tho...

ADD COMMENTlink modified 3.2 years ago • written 3.2 years ago by TriS3.4k

Thanks TriS. I found this ggplots tutorial might be useful. I'll try to run and have a look.

http://stackoverflow.com/questions/5687891/r-how-do-i-display-clustered-matrix-heatmap-similar-color-patterns-are-grouped

ADD REPLYlink written 3.2 years ago by Xinsen Xu30
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1018 users visited in the last hour