In jianhong/colorBlindness: Safe Color Set for Color Blindness

suppressPackageStartupMessages({
  library(colorBlindness)
  library(ggplot2)
  library(reshape2)
})
knitr::opts_chunk$set(warning=FALSE, message=FALSE, eval=TRUE)

Introduction

Color blindness is color vision deficiency, which means people with color blindness may be not aware of differences among colors that are obvious to the people without it. This package is designed to provide a collection of safe colors for plots such as heatmaps, pie chart, and try to find a solution for acceptable figures for all the end users.

Safe colorsets

There are lots of packages already collected lots of safe colors. I list some of them, but not limited here.

• munsell

• viridis and RColorBrewer

• dichromat

• colorblindr

• shades

• ggsci

For this package, I will try to collect as much as possible safe color sets ready to use.

Color vision deficiency (CVD) simulator

The tool used to simulate the color vision deficiency is called CVD simulator. I list some of packages here with CVD simulator.

• Hans Brettel 1997[@brettel1997computerized]

• Francoise Vienot 1999[@vienot1999digital] (implemented in this package and package dichromat )

• Gustavo M. Machado 2009[@machado2009physiologically] (implemented in package colorspace)

• Claire D. McWhite and Claus O. Wilke. colorblindr

How to use this package

Use safe color to plot

To plot figures by colors of this package, there are three steps:

1. select correct safe color set.

2. plot with the color set.

3. check the plot by CVD simulator.

For example if we want to plot a heatmap.

library(colorBlindness)
mat <- matrix(1:81, nrow = 9, ncol = 9)

library(ggplot2)
library(reshape2)
mat1 <- melt(t(mat[9:1, ]))
len <- length(Blue2DarkRed12Steps)-1
mat1$v2 <- cut(mat1$value,
              breaks = seq(0,ceiling(81/len)*len, 
                           length.out = len+1))
ht <- ggplot(mat1) + 
  geom_tile(aes(x=Var1, y=Var2, fill=v2)) + 
  scale_fill_manual(values=Blue2DarkRed12Steps) + 
  theme_bw()
# check the plot by CVD simulator
cvdPlot(ht)

replace the color of a plot to safe color

If you have figures which is not safe to color blindness people, you can try to replace the colors by replacePlotColor function.

cvdPlot(replacePlotColor(displayColors(c("Red", "Green", "blue"))))

If you have a fluorescence image and want to change the colors, you can read the figure and view the effect of color changed figure. If it is acceptable, you can draw it into a file.

library(grid)
library(png)
path <- system.file("extdata", package = "colorBlindness")
img <- readPNG(file.path(path, "fluorescence.microscopy.images.png"))
g <- rasterGrob(img, interpolate=TRUE)
grid.draw(cvdPlot(g, 
                  layout = c("origin", "deuteranope", "protanope",
                             "enhanced", "enhanced.deuteranope", 
                             "enhanced.protanope")))
newImg <- replacePlotColor(g)
outfile <- tempfile(fileext = ".png")
png(outfile)
grid.draw(newImg)
dev.off()

Best approach to make figures

To plot figures for literature, there serveral things we need to consider:

1. Figure format and size;

2. Font family and size;

3. Proper figure colors that are distinguishable by people with colorblindness and for printed in black and white.

setPDFopt function can provide you preset figure size for PDF outputs.

pdf.options()[c("width", "height", "pointsize")]
setPDFopt()
pdf.options()[c("width", "height", "pointsize")]

Collection of safe colors

There are more than 15 Color palettes for color blindness from Martin Krzywinski and et al.

See also GETTING INTO VISUALIZATION OF LARGE BIOLOGICAL DATA SETS and "Controversial Color Use on Maps"[@cynthia1997spectral].

displayAvailablePalette(color="white")

The displayAllColors funciton can be used to show the effect for color vision deficiency. Here we dispaly color collection from Bang Wong[@wong2011points].

displayAllColors(safeColors, color="white")

replacePlotColor effect for Ishihara Color Test.

## replace the color of a file

ishihara29 <- system.file("extdata", "Plate3.png", package = "colorBlindness")
library(png)
img <- readPNG(ishihara29)
g <- rasterGrob(img, interpolate=TRUE)
grid.draw(cvdPlot(g, 
                  layout = c("origin", "deuteranope", "protanope",
                             "enhanced", "enhanced.deuteranope", 
                             "enhanced.protanope")))

Session Info

sessionInfo()

References

jianhong/colorBlindness documentation built on April 17, 2021, 1:31 p.m.