: The creep deformation process is constructed by complex interactions of multiple factors, and the measurement of creep deformation requires enormous economic costs and a long experimental time, so there is a small amount of measurement data. In such a situation, multiple models are often proposed to explain the same experimental data. The coexistence of multiple models based on different physical assumptions makes it difficult to understand the creep deformation process. The purpose of this study is to construct a framework to compare and evaluate coexistence models based on measurement data using the Bayesian model selection framework. Basically, in the creep deformation model, basis functions, usually the exponential function eat or power function tb, corresponding to the primary creep stage and the tertiary creep stage are set, and the time series of the whole creep deformation process is regressed as the linear sum of these bases. The parameters a and b of the basis functions are parameters that must also be optimized to regress the time series data. Izuno et al. proposed the Bayesian model selection method for the creep deformation process, and it was confirmed that the model with the linear term corresponding to the secondary creep stage, called the modified theta method, is a good model compared to the model without