Question: Is there a standard to the range/breaks in z-score coloring for microarray expression data?
0
gravatar for biohack92
5.8 years ago by
biohack92150
United States
biohack92150 wrote:

Just wondering how you determine the range of colors at which you represent your expression data.  Is there a standard or is it case-by-case basis (i.e. expecting minute changes between groups, etc)? 

library(gplots) # contains heatmap.2
plotgenes <- matrix(rnorm(1000),100,10) # generate random data
colbreak = c(seq(-2,-.5,length=100),seq(-.5,.5,length=100),seq(.5,2,length=100))
mycolors <- colorRampPalette(c("red", "black", "green"))(n = 299)
par(cex.main=.8)
ht <- heatmap.2(plotgenes, main=title, col=mycolors, breaks=colbreak,
    density.info="none", trace="none", srtCol = 45, margins=c(7,10), cexCol= .7,
    offsetCol = -.6, dendrogram=c("none"), symm=F,symkey=F,symbreaks=T,
    scale="row", Colv = NA)
R • 2.1k views
ADD COMMENTlink modified 5.8 years ago • written 5.8 years ago by biohack92150
1
gravatar for seidel
5.8 years ago by
seidel6.9k
United States
seidel6.9k wrote:

In my experience it's case by case because data sets have different levels of spread and significance. But I usually expect, and am satisfied to see -3 to 3 (encompassing 8-fold expression differences in log2 scale) because a lot of data sets fit in that range on microarray platforms, and it's comfortable biologically. If it's wider than that - good for you; if it's narrower than that, good luck convincing me. By the way, you might add the following two lines to your code above:

library(gplots) # contains heatmap.2
plotgenes &lt;- matrix(rnorm(1000),100,10) # generate random data

(By the way, in your example, the color breaks are not linear across the range of -2 to 2, which should only be done as long as a legend is also drawn to clearly illustrate the color spectrum).

ADD COMMENTlink modified 4 weeks ago by RamRS25k • written 5.8 years ago by seidel6.9k
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: 1056 users visited in the last hour