Robustness of the best fit for the F₀ (top-left), Y_VSG (top-right), Y_PAH (bottom-left), and Y_BG (bottom-right). We use the best fit SED given by the minimisation method (Y_BestFit) and create a new SED choosing a random flux value in each band from a Gaussian distribution having the best fit flux as a mean value and the observational uncertainty as σ. The best fit parameters of the new fit (Y_New) are then compared with those obtained initially from the minimisation method.