require(colorpatch) require(ggplot2) require(grid) require(gridExtra) vplayout <- function(x, y) viewport(layout.pos.row = x, layout.pos.col = y) set.seed(173)
An important step in analyzing high dimensional data is the inspection of visual maps before the application of automatic analysis techniques. Here we present a method for visualizing fold changes and confidence values within a single diagram. Fold changes (or ratios) naturally occur when comparing a measurement value $B$ with a control condition $A > 0$ as it occurs in analyzing gene expression, agricultural, or financial data. Usually fold changes $r$ are defined by: $$ r = \frac{B - A}{A} $$ or (especially in gene expressions): $$ r = \log_2\frac{B}{A} $$ High dimensional data such as gene expression profiles of different conditions are traditionally visualized as a patch grid showing fold changes (in this case the log ratios) of different genes and multiple samples. The inspection of these maps is known to be prone to errors, if no other information than the fold changes is taken into account [@bilban2002]. The absolute (logarithmic) intensities can be seen as a confidence measure for the observed ratios: $$ a = \frac{1}{2} log_2{(A\cdot B)} \qquad A > 0, B > 0 $$ Other possibilities for computing confidence values may include statistical models.
The colorpatch
package introduces a new bi-variate patch grid visualization for showing fold changes $r_{ij}$ of different samples $j=1\ldots m$ among multiple conditions $i=1\ldots n$ (e.g. genes) together with confidence values $a_{ij}$ within a single visual map.
A psychophysically optimized palette [colorpatch::OptimGreenRedLAB] is used with this
visualization scheme for an optimal visual performance.
The package also contains the code for the optimization of bi-colored color palettes [see @kestler06optim]. As the generation of these palettes is time consuming in the R some of them are pre-computed in the data directory (use the data()
function for loading these palettes):
GreenRedRGB
- linearly scales the green channel and the red channelOptimGreenRedLAB
- perceptually optimized green/red palette in the LAB color spaceOptimBlueYellowLAB
- perceptually optimized blue/yellow palette in the LAB color spaceRe-generation of the palettes can be performed with the following call:
GeneratePalettes()
The colorpatch
package provides color grids of different types:
In the following a random data set is generated
dat <- CreateClusteredData(ncol.clusters = 3, nrow.clusters = 3, nrow = 25, ncol = 15, alpha = 50)
ordered, and pre-processed into a data-frame:
dat <- OrderData(dat) df <- ToDataFrame(dat)
All three approaches are used to visualize the same data set. Cutoff values for fold changes (ratios) and confidence values are set to $0.5$:
thresh.ratio <- 0.5 * max(abs(dat$ratio)) thresh.conf <- 0.5 * max(dat$conf)
For rendering the data the colorpatch
package extends the ggplot2
package with two new statistics stat_colorpatch
and stat_bicolor
:
p <- ggplot(df, aes(ratio = ratio, conf = conf, x = x, y = y)) p <- p + theme_colorpatch(plot.background = "white") + coord_fixed(ratio = 1) p + stat_colorpatch(aes(ratio = ratio, conf = 1, x = x, y = y), thresh.ratio = thresh.ratio, color.fun = ColorPatchColorFun("GreenRedRGB")) + ggtitle("(a) standard green/red") p + stat_bicolor(thresh.ratio = thresh.ratio, thresh.conf = thresh.conf) + ggtitle("(b) HSV bivariate") p + stat_colorpatch(thresh.ratio = thresh.ratio, thresh.conf = thresh.conf) + ggtitle("(c) patch grid")
In the following the uniformity within the LAB color space for the standard RGB palette and the OptimGreenRedLAB palettes are displayed.
data("GreenRedRGB") data("OptimGreenRedLAB") grid.newpage() pushViewport(viewport(layout = grid.layout(2, 1), gp = gpar(fill = "black", col = "black", lwd = 0))) p0 <- PlotUniformity(GreenRedRGB) + ggtitle("GreenRedRGB Uniformity") p1 <- PlotUniformity(OptimGreenRedLAB) + ggtitle("OptimGreenRedLAB Uniformity") print(p0, vp = vplayout(1, 1)) print(p1, vp = vplayout(2, 1)) popViewport()
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.