Multivariate spatial models for lattice data with INLA