Models fit using algorithm='sampling', "meanfield", or
"fullrank" are compatible with a variety of plotting functions from
the rstan package. Each function returns at least one
ggplot object that can be customized further using the
ggplot2 package. The plotting functions described here can be called
using the plot method for stanreg objects
without loading the rstan package. For example usage see
plot.stanreg.