[R] Heatmap construction problems

Struchtemeyer, Chris strucht at okstate.edu
Thu Oct 28 16:42:05 CEST 2010


I am very new to R and don't have any computer program experience
whatsoever.  I am trying to generate a heatmap of the following data:

   
   
 Phylum,AI,AJT,BY,GA,Grt,Sm
 Acidobacteria,0.5,0.7,2.7,0.1,2.6,1.0
 Actinobacteria,33.7,65.1,9.7,2.0,3.9,2.1
 Bacteroidetes,9.7,5.6,0.7,13.2,41.1,21.6
 CCM11b,0.0,0.0,0.0,0.0,0.0,0.1
 Chlamydiae,0.1,0.1,0.0,0.0,1.0,0.2
 Chlorobi,0.0,0.0,0.0,0.0,0.7,1.0
 Chloroflexi,0.1,0.2,0.6,0.2,0.8,0.6
 Cyanobacteria,18.7,0.0,1.0,1.5,0.9,0.3
 Ellusimicrobiales,0.0,0.0,0.0,0.0,0.0,0.0
 Firmicutes,1.0,7.6,8.3,31.9,2.1,6.9
 Gemmatimonadetes,5.0,0.3,0.3,0.0,0.1,0.0
 GN02,0.0,0.0,0.0,0.0,0.0,0.5
 Nitrospirae,0.0,0.2,1.1,0.0,0.0,0.0
 NKB19,0.0,0.0,0.9,0.0,0.0,0.0
 OP8,0.0,0.1,0.0,0.0,0.0,0.0
 OP10,0.6,0.2,0.5,0.0,0.6,0.6
 Planctomycetes,0.9,0.5,6.5,2.2,2.0,2.3
 Alphaproteobacteria,7.8,10.7,21.8,12.2,5.3,26.8
 Betaproteobacteria,9.9,2.8,8.9,21.7,8.3,21.9
 Deltaproteobacteria,0.5,0.2,1.8,2.0,1.2,7.1
 Epsilonproteobacteria,0.0,0.0,0.0,0.1,0.0,0.2
 Gammaproteobacteria,4.0,2.5,8.0,9.4,24.7,5.4
 SC4,0.0,0.0,0.0,0.0,0.7,0.0
 SM2F11,0.0,0.0,0.0,0.0,0.2,0.0
 SPAM,0.0,0.1,0.0,0.0,0.1,0.1
 Synergistes,0.0,0.0,0.0,0.1,0.0,0.0
 Deinococcus-Thermus,0.0,0.0,0.0,0.0,0.0,0.0
 TM6,0.1,0.0,0.0,0.0,0.1,0.0
 TM7,0.0,0.1,0.4,0.0,0.4,0.1
 Verrucomicrobia,3.8,2.1,23.2,2.9,1.3,0.5
 WPS-2,0.0,0.0,0.1,0.0,0.0,0.0
 WS3,0.0,0.0,0.0,0.0,0.0,0.1
 Uncl Bacteria,3.7,1.2,3.7,0.4,1.9,0.8

I am a microbiologist.  What I want to do is construct a heatmap showing the
relative abundance of each phylum.  The far left column of my table contains
all of the phylum names I observed in a set of 6 water samples and each of
the columns to the right contains the relative abundance (%) of each phylum
in each water sample.  I have tried constructing a heatmap using the ggplot
guidelines listed at the following site:
http://learnr.wordpress.com/2010/01/26/ggplot2-quick-heatmap-plotting/

I can generate a heatmap using this method, but would like to alter the
scale.  I would like it so that I can have a little more complex gradient
ranging from 0% to the highest relative abundance that I observe in the
above table (65.1%).  The default scale I get using the link above is just a
relative intensity scale ranging from 1 to 5 (where white represent low
percentages and steelblue represented high percentages).  This is alright
but for phyla that are present at relative abundance of less than 5% all
appear to be white (or non-existant).  Is there anyway to fix this?  Any
help would be greatly appreciated.

Thanks,

Chris



More information about the R-help mailing list