I wanted to overlay the priors of my custom model onto the posterior distributions after running the
rstan function. However, I have not found a way to plot the priors of a Stan model (class object stanfit). There are functions for S3 objects generated using
rstanarm but I could not apply them to a stanfit object (e.g.,
posterior_vs_prior()). I have also tried to access prior information in the same manner as for the extraction of the posterior draws (e.g.,
extract()) but I have not been successful.
Is there a R function that can be used to extract prior distributions and/or plot prior distributions on top of the posteriors of a stanfit object?
If there is not such a function, can I just use the following R code to recreate the beta distribution that I have defined in Stan?
grid = seq(0, 1, .01) plot(grid,dbeta(grid, shape1 = 6, shape2 = 10), type="l", xlab="x", ylab="f(x)")