Bayesian multivariate normal modelling