[BioC] non-symmetric heatmap colors
Arne.Muller at sanofi-aventis.com
Arne.Muller at sanofi-aventis.com
Thu Feb 3 10:30:08 CET 2005
Wolfgang and Sean,
thanks for your hints. The "breaks" argument (described in ?image), did what I wanted. I just had to set all fold change values < -3 to -3 and > 3 to 3 to get a saturation at these limits, otherwise the color beyond these "extreme" values would just be white.
kind regards,
Arne
> -----Original Message-----
> From: bioconductor-bounces at stat.math.ethz.ch
> [mailto:bioconductor-bounces at stat.math.ethz.ch]On Behalf Of Wolfgang
> Huber
> Sent: 03 February 2005 04:49
> Cc: bioconductor at stat.math.ethz.ch
> Subject: Re: [BioC] non-symmetric heatmap colors
>
>
> Hi Arne,
>
> try the following:
>
> library(RColorBrewer)
> library(geneplotter)
> f = colorRamp(brewer.pal(11, "RdBu"))
>
> then f is a functions that takes numbers between 0 and 1 and maps them
> to colors between red and blue.
>
> And
>
> g = function(x) f( (x-x0)/(x1-x0) )
>
> maps numbers between x0 and x1 to these colors.
>
> Best wishes
> Wolfgang
>
>
> Arne.Muller at sanofi-aventis.com wrote:
> > Hello,
> >
> > well, it's about heatmap colors, again ... .
> >
> > Does anybody know a package to create a non-symmetric color
> map with user-choosen mir point and satuaration? I've a
> fold-change matrix, with values from -3 fold (ratio 1/3) to
> +7 fold. I'd like to go from 0.5 (green) to 1 (black) to 3 (red).
> >
> > I think maPallete just assigns the lowest FC to the first
> color and the highest FC to the last color, i.e. an evenly
> distributed color list, right?
> >
> > Any hints are welcome,
> >
> > Arne
>
> -------------------------------------
> Wolfgang Huber
> European Bioinformatics Institute
> European Molecular Biology Laboratory
> Cambridge CB10 1SD
> England
> Phone: +44 1223 494642
> Fax: +44 1223 494486
> Http: www.ebi.ac.uk/huber
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
>
More information about the Bioconductor
mailing list