While I’ve used great packages like
rstanarm in the past, I’ve begun working on writing my own models using
Rstan. Although I’m starting out with just basic models, one issue I’ve already run into is not having the same functions I can rely on for different aspects of my workflow.
For example, what if I want draws from the posterior predictive distribution? Using
rstanarm I’d just use
posterior_predict() (and in
brms I’d use
tidybayes::add_predicted_draws()), so are there analogous functions for models written in