I’m looking at posterior predictive distributions for three models. As part of the comparison, I produced boxplots using ppc_boxplot for the observed data and 10 randomly posterior predictive distributions for each model. My problem is that for two of the models, the boxplots all have the same limits on the vertical axis, while for the third model I get different limits. Is there a way to specify the y-axis limits (like the ylim= parameter for general graphing in R)?
The objects produced by the bayesplot package are ggplot2 objects, and so you can customize them the same way you would other plots produced by ggplot2. In particular, in this case I think you can do something like
bayesplot::ppc_boxplot(model1) + ylim(0, 5) bayesplot::ppc_boxplot(model2) + ylim(0, 5) bayesplot::ppc_boxplot(model3) + ylim(0, 5)
to have the vertical axis run from zero to five for all three plots. For more information, see the ggplot2 documentation on
That solved it… thanks!!!