We discuss the advantages of a statistical technique, based on the simplex minimization and on the parametric bootstrap, with respect to the standard $\chi^2$-minimization method.
We present the application of this technique for the extraction of the proton scalar dipole dynamical polarizabilities from real Compton scattering data. In particular, we discuss how this technique is able to provide realistic probability distributions of the fitted parameters, with no a-priori assumptions on their shape, and to include the systematic errors in the minimization procedure in a straightforward way.