Добавление сводной информации к графику плотности, созданному с помощью ggplot

У меня есть график плотности, и я хотел бы добавить некоторую сводную информацию, такую ​​как нанесение линии на медиане и закрашивание 90% вероятных интервалов (5-й и 95-й квантили). Есть ли способ сделать это в ggplot?

Вот тип сюжета, который я хотел бы резюмировать:

Я могу понять, как провести линию от y = 0 до y = плотности (медиана (x)), но мне не ясно, смогу ли я заштриховать график с 90% доверительным интервалом. В качестве альтернативы я мог бы добавить горизонтальную коробчатую диаграмму над графиком плотности, но неясно, как повернуть прямоугольную диаграмму отдельно, не вращая график плотности вместе с ним.

x <- as.vector(rnorm(10000))
d <- as.data.frame(x=x)
library(ggplot2)
ggplot(data = d) + theme_bw() + 
  geom_density(aes(x=x, y = ..density..), color = 'black')