Plotting subsets of parameters with Rstan

I’m using Rstan. I have something like this:

simplex[K] theta[Q];

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. theta[,1:3,c(2,4)])?

So it would be great if I could do plot(input, pars = c("theta[,2]"), or something similar

Edit: fixed typo

I think plot(input, pars = "theta[2]") works.