Heat capacity data of many crystalline solids can be described in a physically sound manner by Debye–Einstein integrals in the temperature range from (Formula presented.) to (Formula presented.). The parameters of the Debye–Einstein approach are either obtained by a Markov chain Monte Carlo (MCMC) global optimization method or by a Levenberg–Marquardt (LM) local optimization routine. In the case of the MCMC approach the model parameters and the coefficients of a function describing the residuals of the measurement points are simultaneously optimized. Thereby, the Bayesian credible interval for the heat capacity function is obtained. Although both regression tools (LM and MCMC) are completely different approaches, not only the values of the Debye–Einstein parameters, but also their standard errors appear to be similar. The calculated model parameters and their associated standard errors are then used to derive the enthalpy, entropy and Gibbs energy as functions of temperature. By direct insertion of the MCMC parameters of all (Formula presented.) computer runs the distributions of the integral quantities enthalpy, entropy and Gibbs energy are determined.