The accuracy and precision of high-energy spallation models are key issues for the design and development of new applications and experiments. We present a method to estimate model parameters and associated uncertainties by leveraging the Bayesian version of the Generalised Least Squares method, which enables us to incorporate prior knowledge on the parameter values. This approach is designed to adjust parameters based on experimental data, accounting for experimental uncertainty information, and providing uncertainties for all adjusted parameters. This approach is designed in order both to improve the accuracy of models through the modification of free parameters of these models, which results in a better reproduction of experimental data, and to estimate the uncertainties of these parameters and, by extension, their impacts on the model output. We aim at demonstrating the Generalised Least Square method can be applied in the case of Monte Carlo models. We present a proof-of-concept for Monte Carlo models in the specific case of nuclear physics with the model combination INCL/ABLA. We discuss the challenges in the application of this method to high-energy spallation models, notably the large runtime and the stochasticity of the models. Our results indicate this framework can also be applied to analogous situations where parameters of a computationally expensive Monte Carlo code should be inferred/improved.