was doing this with bayesplot::mcmc_trace using the `variables` argument but it first extracts all samples and then subsets them

plot.brmsfit does not have this limitation. thanks for the answer, Stephen!