Ferroelectric relaxors are complex materials with distinct properties. The understanding of their dielectric susceptibility, which strongly depends on both temperature and probing frequency, has been a challenge to researchers for many years. Here we report a macroscopic and phenomenological approach based on statistical modeling to investigate how the dielectric response of a relaxor depends on temperature. Employing the MaxwellBoltzmann distribution and considering temperature dependent dipolar orientational polarizability, we propose a minimum statistical model and specific equations to understand and fit numerical and experimental dielectric responses versus temperature. We show that the proposed formula can successfully fit the dielectric response of typical relaxors, including Ba(Zr,Ti)O 3 , 0.87Pb(Zn 1/3 Nb 2/3 ) 0.87 O 3 -0.13PbTiO 3 , 0.95Pb(Mg 1/3 Nb 2/3 )O 3 -0.05Pb(Zr 0.53 Ti 0.47 )O 3 , and Bi-based compounds, which demonstrate the general applicability of this approach.