In our novel approach, we use the previous VBScan recursively on the map with case-control point event data. At each recursive step, we compute two functions: (i) the likelihood ratio of the multiple components and (ii) the likelihood ratio increase since the previous step. As the first function always increases monotonically with every added component to the partition, the last function is used as a measure of the cost-benefit of adding a further region to the partition. This is done employing a multicriteria decision process, determining the nondominated partition solutions. Through Monte Carlo replications under null hypothesis, we compute the significance of the nondominated solutions and choose the best partition.