I’m using Rstan. I have something like this:
And I can use Rstan’s
rstan::plot(input, pars = c("theta")) . But, this is not good when any dimension of theta gets large, and I’d like to do comparisons.
Is there any way I can get a similar “boxplot” for, in R:
theta[,3] where I index later? What about for a range of indices (i.e.
So it would be great if I could do
plot(input, pars = c("theta[,2]"), or something similar
Edit: fixed typo