Как создать диаграмму с баром ошибок, используя ggplot для двух категориальных переменных в R?

Question

Я хочу построить гистограмму для двух категориальных переменных, отображающих CI для пропорции, используя ggplot.

Это то, что я пробовал.Переменные сокращения и ясности набора данных diamond имеют больше категорий, чем мои данные, но это не имеет значения.Код графика, приведенный ниже, довольно близок к тому, что я хочу получить в качестве окончательного графика, я просто хочу добавить 95% -й КИ от пропорции и панели ошибок к графику.

N <- diamonds %>%
  ggplot(aes(x = cut, group = clarity)) + 
  geom_bar(aes(y = ..prop.., fill = clarity), stat = "count") +
  geom_text(aes(label = scales:: percent(..prop..), 
                y = ..prop..), stat = "count", vjust = -.5) + 
  labs(y = "Percent", title = "Distribution of diamond clarity by cut",fill = "Clarity") + 
  facet_grid(~clarity) + 
  scale_y_continuous(labels = scales::percent)

N

Заранее спасибо.