We derive an equilibrium statistical theory for the macroscopic description of a ferromagnetic material at positive finite temperatures. Our formulation describes the most-probable equilibrium macrostates that yield a coherent deterministic large-scale picture varying at the size of the domain, as well as it captures the effect of random spin fluctuations caused by the thermal noise. We discuss connections of the proposed formulation to the Landau-Lifschitz theory and to the studies of domain formation based on Monte Carlo lattice simulations.