Bayesian modelling strategies for spatially varying regression coefficients: A multivariate perspective for multiple outcomes