Inferring the position of functionally active regions from a multi-subject fMRI dataset involves the comparison of the individual data and the inference of a common activity model. While voxel-based analyzes, e.g. Random Effect statistics, are widely used, they do not model each individual activation pattern. Here, we develop a new procedure that extracts structures individually and compares them at the group level. For inference about spatial locations of interest, a Dirichlet Process Mixture Model is used. Finally, inter-subject correspondences are computed with Bayesian Network models. We show the power of the technique on both simulated and real datasets and compare it with standard inference techniques.