Multivariate disease risk mapping: a geostatistical approach