Spatial prediction in a multivariate setting