Approximate Bayesian inference for multivariate point pattern analysis in disease mapping