6

Aim

I am trying to fill the area between two lines in a plot generated with ggplot in R. I would like to fill everything between the lines above of the horizontal line with a different color than below the horizontal line.

I succeeded to fill everything between the two lines with a single color, however, I did not manage to differentiate above and below the vertical line by two different colors.

Code

set.seed(123) 

# Load packages
library(tidyverse)

# Create sample dataframe
df <- data.frame(x=seq(1,50,1),y=runif(50, min = 0, max = 10))

# Generate plot
ggplot(data = df, aes(x = x, y = y)) +
  geom_line() +
  geom_hline(yintercept = 5) +
  theme_classic() +
  geom_ribbon(aes(ymin=5,ymax=y), fill="blue")

Question

How do I fill the space above and below the horizontal line with a different color?

tjebo
  • 21,977
  • 7
  • 58
  • 94
user213544
  • 2,046
  • 3
  • 22
  • 52
  • I do not fully agree. This question concerns values above or below a specific line, while the other concerns y=0. Therefore, the initial steps of data preparation are different (I.e., this would be a more generalizable approach, since it is not dependent on switching of signs). – user213544 Jul 14 '22 at 08:15
  • the specific value in the linked thread is "0". The answers are extendable to any value – tjebo Jul 14 '22 at 08:15
  • Well, the interpolation in the linked thread uses a switch in sign of the number, right? Then I guess it will not work on any other positive value? – user213544 Jul 14 '22 at 08:20
  • see another linked thread - https://stackoverflow.com/questions/72972363/r-ggplot2-fill-only-the-area-between-a-line-and-a-reference-value?noredirect=1&lq=1 - the solution works very well. – tjebo Jul 14 '22 at 08:21
  • I think I am beginning to understand that you're referring to the accepted answer. But the linked thread does not only have this answer, it contains a series of useful answers that are all applicable to the problem - which is effectively the one presented here. Closing this question here as a duplicate does not mean that this here is a bad question - it is not! - but it will ensure that future answers to the actual problem are found concentrated in one thread. – tjebo Jul 14 '22 at 08:26

2 Answers2

8

You can calculate the coordinates of the points where the two lines intersect & add them to your data frame:

m <- 5 # replace with desired y-intercept value for the horizontal line

# identify each run of points completely above (or below) the horizontal
# line as a new section
df.new <- df %>%
  arrange(x) %>%
  mutate(above.m = y >= m) %>%
  mutate(changed = is.na(lag(above.m)) | lag(above.m) != above.m) %>%
  mutate(section.id = cumsum(changed)) %>%
  select(-above.m, -changed)

# calculate the x-coordinate of the midpoint between adjacent sections
# (the y-coordinate would be m), & add this to the data frame
df.new <- rbind(
  df.new,
  df.new %>%
    group_by(section.id) %>%
    filter(x %in% c(min(x), max(x))) %>%
    ungroup() %>%
    mutate(mid.x = ifelse(section.id == 1 | 
                            section.id == lag(section.id), 
                          NA,
                          x - (x - lag(x)) / 
                            (y - lag(y)) * (y - m))) %>% 
    select(mid.x, y, section.id) %>%
    rename(x = mid.x) %>%
    mutate(y = m) %>%
    na.omit())

With this data frame, you can then define two separate geom_ribbon layers with different colours. Comparison of results below (note: I also added a geom_point layer for illustration, & changed the colours because the blue in the original is a little glaring on the eyes...)

p1 <- ggplot(df,
             aes(x = x, y = y)) +
  geom_ribbon(aes(ymin=5, ymax=y), fill="dodgerblue") +
  geom_line() +
  geom_hline(yintercept = m) +
  geom_point() + 
  theme_classic()

p2 <- ggplot(df.new, aes(x = x, y = y)) +      
  geom_ribbon(data = . %>% filter(y >= m),
              aes(ymin = m, ymax = y), 
              fill="dodgerblue") +
  geom_ribbon(data = . %>% filter(y <= m),
              aes(ymin = y, ymax = m), 
              fill = "firebrick1") +
  geom_line() +
  geom_hline(yintercept = 5) +
  geom_point() +      
  theme_classic()

result

Z.Lin
  • 28,055
  • 6
  • 54
  • 94
1

fill three layers in order, for a particularly ugly result:

# Generate plot
ggplot(data = df, aes(x = x, y = y)) +
  geom_line() +
  geom_hline(yintercept = 5) +
  theme_classic() +
  geom_ribbon(aes(ymin=y,ymax=10), fill="green")+
  geom_ribbon(aes(ymin=0,ymax=y), fill="yellow")+
  geom_ribbon(aes(ymin=5,ymax=y), fill="blue")
Jorrit G
  • 509
  • 2
  • 5
  • 6
  • Thanks for the contribution! However, I am trying to color between the lines and not 'outside' of the lines. – user213544 Feb 14 '19 at 10:43