Disease risk mapping: a multivariate geostatistical approach