Vector Autoregressions?

I’ve posted a new related topic here:

vector autoregressions of independent series are doable by passing a 2d set of priors for the rho argument to pm.AR. I’m not sure how we let these series load on each other’s lagged values though.