Method for selecting an operating mode of devices (111, 112, 113) for thermally regulating a building, from an accumulative first operating mode, configured to maximize external thermal inputs and minimize thermal losses, and a dissipative second operating mode, configured to minimize external thermal inputs and maximize thermal losses, each thermally regulating device being intended to regulate the temperature inside the building, called the interior temperature, the method being executed by a control unit of the building and comprising, for each of the operating modes: - obtaining estimates of interior temperature in light of predefined parameters defined beforehand using a model for predicting interior temperature, - defining a prediction error representative of discrepancies between measured interior-temperature values and the obtained estimates of interior temperature, each discrepancy being defined at a time at which the measured interior temperature is determined when said operating mode is selected, the prediction error being dependent on the predefined parameters, - determining a set of optimum values associated with said operating mode, containing one value associated with each predefined parameter, and allowing the prediction error to be minimized, - determining (301, 311) a prediction sequence associated with said operating mode, the prediction sequence containing estimates of interior temperature calculated, from an initial time t i, using the model for predicting interior temperature and with application of the determined set of optimum values associated with the operating mode, - determining (303, 313) a comfort criterion associated with the operating mode, said comfort criterion being representative of a feeling of comfort and characterizing whether or not a minimum comfort-temperature threshold and a maximum comfort-temperature threshold are crossed by the estimates of interior temperature of the prediction sequence, the method further comprising: - selecting (321, 322) the operating mode for which the associated determined comfort criterion is representative of the greatest feeling of comfort, characterized in that the predicting model is dependent on a first predefined parameter Φ, which is characteristic of a phase shift, and on a second predefined parameter ε corresponding to a damping factor, the estimate of the internal temperature T int estim at a time t n, located between the initial time t i and a final time t f, being obtained on the basis of an interior-temperature value T int mes known at the initial time t i, of an exterior-temperature value T ext known or estimated at the initial time t i and of a predicted exterior-temperature value T ext at time t n, according to: T int estim t n = T int mes t i + ε. T ext t n + Φ T ext t i. - obtaining estimates of interior temperature in light of predefined parameters defined beforehand using a model for predicting interior temperature, - defining a prediction error representative of discrepancies between measured interior-temperature values and the obtained estimates of interior temperature, each discrepancy being defined at a time at which the measured interior temperature is determined when said operating mode is selected, the prediction error being dependent on the predefined parameters, - determining a set of optimum values associated with said operating mode, containing one value associated with each predefined parameter, and allowing the prediction error to be minimized, - determining (301, 311) a prediction sequence associated with said operating mode, the prediction sequence containing estimates of interior temperature calculated, from an initial time t i, using the model for predicting interior temperature and with application of the determined set of optimum values associated with the operating mode, - determining (303, 313) a comfort criterion associated with the operating mode, said comfort criterion being representative of a feeling of comfort and characterizing whether or not a minimum comfort-temperature threshold and a maximum comfort-temperature threshold are crossed by the estimates of interior temperature of the prediction sequence, the method further comprising: - selecting (321, 322) the operating mode for which the associated determined comfort criterion is representative of the greatest feeling of comfort, - selecting (321, 322) the operating mode for which the associated determined comfort criterion is representative of the greatest feeling of comfort, characterized in that the predicting model is dependent on a first predefined parameter Φ, which is characteristic of a phase shift, and on a second predefined parameter ε corresponding to a damping factor, the estimate of the internal temperature T int estim at a time t n, located between the initial time t i and a final time t f, being obtained on the basis of an interior-temperature value T int mes known at the initial time t i, of an exterior-temperature value T ext known or estimated at the initial time t i and of a predicted exterior-temperature value T ext at time t n, according to: T int estim t n = T int mes t i + ε. T ext t n + Φ T ext t i.