Bayesian neural networks with correlating residuals

In a multivariate regression problem it is often assumed that residuals of outputs are independent of each other. In many applications a more realistic model would allow dependencies between the outputs. In this paper we show how a Bayesian treatment using the Markov chain Monte Carlo method can allow for a full covariance matrix with multilayer perceptron neural network.