user3389288 user3389288 - 2 years ago 103
R Question

Changing the alpha values in R{graphics} while the colour argument is used

I'm so used to doing this in ggplot2 that I'm having a very hard time figuring out how to specify alpha values using R base graphics, while the col= argument in plot() is used for assigning a colour type to a categorical variable.

Using the iris data set (although in this context it doesn't really make sense why we would need to change the alpha values)

data(iris)
library(ggplot2)
g <- ggplot(iris, aes(Sepal.Length, Petal.Length)) + geom_point(aes(colour=Species), alpha=0.5) #desired plot

plot(iris$Sepal.Length, iris$Petal.Length, col=iris$Species) #attempt in base graphics


What about mapping another variable to the alpha value using {graphics}?
For example in ggplot2:

g2 <- ggplot(iris, aes(Sepal.Length, Petal.Length)) + geom_point(aes(colour=Species, alpha=Petal.Width))


Any help is appreciated!

Answer Source

Adjusting alpha is pretty easy with adjustcolor function:

COL <- adjustcolor(c("red", "blue", "darkgreen")[iris$Species], alpha.f = 0.5)
plot(iris$Sepal.Length, iris$Petal.Length, col = COL, pch = 19, cex = 1.5) #attempt in base graphics

enter image description here

Mapping alpha to variable requires a bit more hacking:

# Allocate Petal.Length to 7 length categories
seq.pl <- seq(min(iris$Petal.Length)-0.1,max(iris$Petal.Length)+0.1, length.out = 7)

# Define number of alpha groups needed to fill these
cats <- nlevels(cut(iris$Petal.Length, breaks = seq.pl))

# Create alpha mapping
alpha.mapping <- as.numeric(as.character(cut(iris$Petal.Length, breaks = seq.pl, labels = seq(100,255,len = cats))))

# Allocate species by colors
COLS <- as.data.frame(col2rgb(c("red", "blue", "darkgreen")[iris$Species]))

# Combine colors and alpha mapping
COL <- unlist(lapply(1:ncol(COLS), function(i) {
  rgb(red = COLS[1,i], green = COLS[2,i], blue = COLS[3,i], alpha = alpha.mapping[i], maxColorValue = 255)
  }))

# Plot
plot(iris$Sepal.Length, iris$Petal.Length, col = COL, pch = 19, cex = 1.5) 

enter image description here

Recommended from our users: Dynamic Network Monitoring from WhatsUp Gold from IPSwitch. Free Download