У меня есть этот фреймворк:
set.seed(1)
x <- c(rnorm(50, mean = 1), rnorm(50, mean = 3))
y <- c(rep("site1", 50), rep("site2", 50))
xy <- data.frame(x, y)
И я сделал этот график плотности:
library(ggplot2)
ggplot(xy, aes(x, color = y)) + geom_density()
Для site1
Мне нужно затенять область под кривой, чтобы > 1% данных. Для site2
мне нужно затенять область под кривой, 75% данных.
Я ожидаю, что сюжет будет выглядеть примерно так (фотошоп). Пройдя через переполнение стека, я знаю, что другие спрашивали, как затенять часть области под кривой, но я не могу понять, как затенять область под кривой группой.