Modelling geo-referenced categorical data: a multivariate geostatistical approach