Credible regions of modelled galaxy stellar mass functions: The model number densities as calculated from Equation (3) matched to observational data. The credible regions cover (in decreasing order of colour saturation) 68%, 95% and 99.7% of the posterior distributions. The shaded area in the background marks the regime where the two feedback mechanisms contribute about equally to the stellar mass – halo mass relation given by Equation (4), the white area on the low-mass side is dominated by stellar feedback, $({\left(\raisebox{1ex}{$M_{\mathrm{h}}$}\!\left/ \!\raisebox{-1ex}{$M_{\mathrm{c}}^{\star }$}\right.\right)}^{{\delta }_{\star }+{\gamma }_{\star }}<0.1)$, while the area on the high-mass side is AGN feedback-dominated, ${\left(\raisebox{1ex}{$M_{\mathrm{h}}$}\!\left/ \!\raisebox{-1ex}{$M_{\mathrm{c}}^{\star }$}\right.\right)}^{{\delta }_{\star }+{\gamma }_{\star }}>10$. The regime borders are calculated as the median value from a sample of GSMFs drawn from the posterior.