Joint estimation and regularized aggregation of brain network in FMRI data