I am not sure how to make the custom plot you ask for, but the histogram seems to do a good job already. The density plot isn’t good for this data (as you noted).
I would be interested to see how the zero-one-mid-inflated beta model works on this data from the vignette, if that spike near the midpoint is right at 0.5 and represents, say, neutral responses from a slider scale.
Note this comment from Michael Betancourt here regarding mixture models:
"In other words because the the inflated and non-inflated points are essentially modeled by separate data generating processes the entire model decomposes into a binomial model for the total inflated observations and a baseline model for the non-zero observations.
Because these two models are independent we can fit them independently or jointly, whichever is more convenient. In particular if the inflated counts are just a nuisance then we can ignore them entirely and just fit the non-inflated observations directly without any consideration of the inflation!"
That’s worth thinking about when you are fitting some slider scale data where people tend to plop the slider at units of 10 or whatever, with few responses in between.