You may want to use ggplot's facet_wrap() function, with a symptom variable for the subplots. See section 3.5 in R for Data Science. https://r4ds.had.co.nz/index.html
1- ggplot2 library for make the bar plot and add errors bars
2- ggstatsplot library for add top bars and * simbol
3- like said @EconProf facet_wrap() function of ggplot2 or gridExtra packages