12

In the following example I create two series of points and plot them using ggplot2. I also highlight several points based on their values

library(ggplot2)
x <- seq(0, 6, .5)
y.a <- .1 * x -.1
y.b <- sin(x)
df <- data.frame(x=x, y=y.a, case='a')
df <- rbind(df, data.frame(x=x, y=y.b, case='b'))
print(ggplot(df) + geom_point(aes(x, y), color=ifelse(df$y<0, 'red', 'black')))

And here is the result

First result

Now I want to separate the two cases into two facets, keeping the highlighting scheme

> print(ggplot(df) + geom_point(aes(x, y), color=ifelse(df$y<0, 'red', 'black')) + facet_grid(case ~. ,))
Error: Incompatible lengths for set aesthetics: colour

How can this be acheived?

tjebo
  • 21,977
  • 7
  • 58
  • 94
Boris Gorelik
  • 29,945
  • 39
  • 128
  • 170
  • related https://stackoverflow.com/questions/15804504/how-to-conditionally-highlight-points-in-ggplot2-facet-plots-mapping-color-to – tjebo Aug 11 '22 at 14:11

2 Answers2

27

You should put color=ifelse(y<0, 'red', 'black') inside the aes(), so color will be set according to y values in each facet independently. If color is set outside the aes() as vector then the same vector (with the same length) is used in both facets and then you get error because length of color vector is larger as number of data points.

Then you should add scale_color_identity() to ensure that color names are interpreted directly.

ggplot(df) + geom_point(aes(x, y, color=ifelse(y<0, 'red', 'black'))) + 
   facet_grid(case ~. ,)+scale_color_identity()

enter image description here

Didzis Elferts
  • 95,661
  • 14
  • 264
  • 201
4

Instead of using scale_..._identity, one can also wrap the color (and fill) aesthetic in I(). It also requires having color defined in aes.

I came across this question, where the OP I guess kind of accidentally made use of I()... ggplot color is not automatically coloring based on group

Not sure I would every make use of that, but I find this kind of fun.

library(ggplot2)
x <- seq(0, 6, .5)
y.a <- .1 * x -.1
y.b <- sin(x)
df <- data.frame(x=x, y=y.a, case='a')
df <- rbind(df, data.frame(x=x, y=y.b, case='b'))

ggplot(df) + 
  geom_point(aes(x, y, color= I(ifelse(y < 0, 'red', 'black')))) +
  facet_grid(case ~. )

Created on 2020-07-01 by the reprex package (v0.3.0)

tjebo
  • 21,977
  • 7
  • 58
  • 94