Shinystan style plots directly in R

Hello,

I love some of the shinystan plots, but would like to customize them for some printed reports. Would be great to figure out a way to re-create something similar in R (maybe with ggplot2.)

The main one I want is the plot shown on the “Parameters plot” tab of the “Estimate” section. It has each parameter’s density represented a s single line, with a thicker portion for the 50% interval. That’s a nice way to communicate the sample variance to non-stats people.

I’m not sure what that plot is called, so don’t know where to look.

Any suggestions on how to create this type of plot?

bayesplot::mcmc_intervals

Thanks Ben!