The measurand value, the conclusions, and the decisions inferred from measurements may depend on the models used to explain and to analyze the results. In this paper, the problems of identifying the most appropriate model and of assessing the model contribution to the uncertainty are formulated and solved in terms of Bayesian model selection and model averaging. As computational cost of this approach increases with the dimensionality of the problem, a numerical strategy, based on multimodal ellipsoidal nested sampling, to integrate over the nuisance parameters and to compute the measurand post-data distribution is outlined. In order to illustrate the numerical strategy, by use of MATHEMATICA an elementary example concerning a bimodal, two-dimensional distribution has also been studied.