The high- and low-energy limits of a chain of coupled rotators are integrable and correspond respectively to a set of free rotators and to a chain of harmonic oscillators. For intermediate values of the energy, numerical calculations show the agreement of finite time averages of physical observables with their Gibbsian estimate. The boundaries between the two integrable limits and the statistical domain are analytically computed using the Gibbsian estimates of dynamical observables. For large energies the geometry of nonlinear resonances enables the definition of relevant 1.5-degrce-of-freedom approximations of the dynamics. They provide resonance overlap parameters whose Gibbsian probability distribution may be computed. Requiring the support of this distribution to be right above the large-scale stochasticity threshold of the 1.5-degree-of-freedom dynamics yields the boundary at the large-energy limit. At the low-energy limit, the boundary is shown to correspond to the energy where the specific heat departs from that of the corresponding harmonic chain.