When learning Bayesian statistics, I feel like its common to see some version of the plot that shows a prior density, a density of the likelihood, and a posterior density, somewhere in between (e.g. https://goo.gl/images/Ksg7uR or

)I was trying to build of version of this graph for one of the parameters in my model - using RStan. Getting the prior curve is quite easy (since the prior is just normal with a specific mean/sd, plotting it with `stat_function`

), and the posterior is pretty easy (just extract the sample, and plot it with `stat_density`

) - but is there a way to get the likelihood density? Is there something I’d need to put into `generated quantities`

or something like that?

My first attempt was to just run a new model, with an extremely uninformative prior (like `normal(0, 100)`

, and assume that the posterior of that represents the likelihood of the data, but I’m getting a nonsensical result, where the posterior is higher than either the prior or the likelihood.

Is there a way to do this? Thanks!