The only issue with this is that this is a toy dataset/model and I think adapting the full version of what I’m looking for to rstanarm/brms is beyond my ability. See this earlier question for a sense of where I am trying to get to.
Will give this a shot and report back!