Hi @Arvindh_Somanathan,
I used some fake data, but I think this gets you most of the way there:
library(tidyverse)
tibble(
y = letters,
x = runif(length(letters), 0, 4),
gp = rep(paste0('Group ', 1:5), length.out = length(letters))
) %>%
ggplot(aes(x, y)) +
annotate('rect', xmin = 0.5, xmax = 2, ymin = -Inf, ymax = Inf,
fill = 'black', alpha = 0.1) +
geom_vline(xintercept = 1, lty = 2) +
geom_point(aes(color = x > 1), size = 3) +
geom_text(aes(label = round(x, 1)), nudge_x = 0.15) +
scale_x_continuous(
breaks = seq(0, 4, by = 0.5)
) +
scale_color_manual(values = c('darkorange', 'darkblue')) +
labs(y = NULL, x = 'National Benefits Ratio vs. CMAQ') +
facet_grid(rows = vars(gp), switch = 'y', scales = 'free') +
theme_minimal(14) +
theme(
legend.position = 'none',
strip.placement = 'outside',
panel.spacing = unit(0, 'in'),
strip.text.y.left = element_text(angle = 0, vjust = 1),
panel.border = element_rect(fill = 'transparent')
)
