A Bayesian direct method implementation to fit emission line spectra: application to the primordial He abundance determination