In the context of neutron science, Bayesian inference methods have been recently implemented within the MANTID framework [Monserrat D et al. 2015 J. Phys. Conf. Ser. 663 012009 (2015)]. In this contribution, we highlight the advantages of this software package for robust data analysis and subsequent model selection. To this end, we use the celebrated Rosenbrock function to illustrate its merits and strengths relative to classical fitting algorithms. We also introduce the latest additions implemented in MANTID, with a view to increasing its user friendliness as well as stimulating wider use. These include simulated-annealing schemes to reduce the need for initial guesses, as well as new options for multidimensional fitting. © Published under licence by IOP Publishing Ltd.