Heatmap of RPKM analysis
1
0
Entering edit mode
Jojocc • 0
@jojocc-19971
Last seen 5.8 years ago

I made an RPKM analysis of my data and I generated this result

head((RPKM_analysis))
                    a          b          c          d          e          f              
gene0         8.5099049 5.5732562   5.0730217  6.4162466  25.075330   19.163763       
gene1         2.6056052 0.4082469   4.4356482  0.9731503   3.508077    41.594664        
gene2         0.4727903 0.3678376   0.3397144  0.1221217   4.475466    2.936647         
gene3         0.8390143 0.3199489   1.2082854  0.1856811  14.063714   4.129234        
gene4         2.2246130 2.3719177   2.3547895  0.7456295  22.917552   7.871053        
gene6       150.0407639 9.5132677 266.7380691 76.9122813 100.517840  61.350209 
                       g           i       j           k                             
gene0             34.667379 13.056078 90.2201082    8.3383381
gene1              4.642868  2.658317  0.5275397     0.6345467
gene2              1.864922  1.487531  8.1337453     24.3506878
gene3              2.594216  2.194451  6.2621722     49.8774268
gene4              4.922619  5.944980 41.6279899    119.1270556
gene6            582.504495 12.147723 74.0451493   11.8167852

It is possible to create a heatmap of this analysis? I used this script

heatmap_RPKM <- heatmap(as.matrix(RPKM_analysis))

but the problem is that it doesn't show me rownames and colnames. Anyone can help me?

normalization • 2.1k views
ADD COMMENT
0
Entering edit mode
Kevin Blighe ★ 4.0k
@kevin
Last seen 21 days ago
Republic of Ireland

You are using the basic heatmap() function. If my memory serves me correctly, by default, row and column names are not activated in this function. You can specify them with labRow and labCol.

You are also using as.matrix(), but I believe you should use data.matrix().

If you want to 'upgrade' to heatmap.2(), then take a look at my answer on Biostars: https://www.biostars.org/p/284721/#284723

heatmap

Kevin

ADD COMMENT
0
Entering edit mode

If the columns of RPKM_analysis are numeric, then as.matrix() or data.matrix give the same result.

If that's not the case, then making a heatmap probably isn't appropriate.

ADD REPLY

Login before adding your answer.

Traffic: 730 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