Residual values for B_{ℓ1 + ℓ2 + ℓ3 = N} for N = 30 and N = 60, for the GNILC reconstructed maps as a function of the channel number. We plot different values of m_AIC (top and bottom figures), and the residuals are a few orders of magnitude smaller than those presented in Fig. 14, independently of which m_AIC is chosen for our bispectrum test for sub-optimal foreground extractions parameters.