Search code examples
rbayesianmcmcstanrstanarm

How to extract stan code from rstanarm object

Is there a possibility to extract the stan code used for the MCMC sampling in rstanarm?

I would like to compare my own parametrisation of a model and prior choices to the one used in rstanarm.

Solution

  • You can execute

    library(rstanarm)
example(example_model)
rstan::get_stanmodel(example_model$stanfit)

    to see the Stan code or look at it on GitHub. However, it was not intended to be read by humans and is not particularly readable due to the #include statements and the large number of options supported.