For what is regarding the fitting procedure, in order to compare q-Gaussians and pseudo-Voigtians, the procedure must be the same. For comparison with a Voigt function, we must consider a further intrinsic difficulty. When we use a Voigt function, we necessarily need a numerical approximation of the convolution integral and consequently we can obtain different results by changing this approximation.