The Sagdeev-Zaslavski (SZ) equation for wave turbulence is analytically derived, both in terms of generating function and of multi-point pdf, for weakly interacting waves with initial random phases. When also initial amplitudes are random, the one-point pdf equation is derived. Such analytical calculations remarkably agree with results obtained in totally different fashions. Numerical investigations of the two-dimensional nonlinear Schrödinger equation (NLSE) and of a vibrating plate prove that: (i) generic Hamiltonian 4-wave systems rapidly attain a random distribution of phases independently of the slower dynamics of the amplitudes, vindicating the hypothesis of initially random phases; (ii) relaxation of the Fourier amplitudes to the predicted stationary distribution (exponential) happens on a faster timescale than relaxation of the spectrum (Rayleigh-Jeans distribution); (iii) the pdf equation correctly describes dynamics under different forcings: the NLSE has an exponential pdf corresponding to a quasi-gaussian solution, like the vibrating plates, that also show some intermittency at very strong forcings.Introduction Dispersive waves are ubiquitous in nature, and their nonlinear interactions make them intriguing and challenging [1,2]. Wave Turbulence is the theory that describes the statistical properties of large numbers of incoherent interacting waves, with tools such as the wave kinetic equation analytically derived in the late sixties. This equation describes the evolution of the wave spectrum in time, when homogeneity and weak nonlinearity are assumed [3][4][5]. It has been applied to numerous phenomena, including ocean waves [6][7][8], capillary waves [9, 10] Alfvén waves [11], optical waves [12] and solid oscillations [13][14][15][16][17][18]. It is the analogue of the Boltzmann equation for classical particles and it allows the RayleighJeans equilibrium state as well as non-equilibrium solutions, in terms of Kolmogorov-Zakharov (KZ) spectra [19].To characterise the invariant measure of the dynamics, that is to find the complete statistical description concerning all quantities of interest, an important step has been taken by Sagdeev and Zaslavski [20], who obtained the Brout-Prigogine equation for the probability density function (pdf) of wave turbulence [21]. More recently, this statistical framework has been nicely revisited using the diagrammatic technique [4] and performing analytical calculations, in the 3-wave case [22][23][24]. Interestingly, many experimental and theoretical results have shown that deviations from wave-turbulence predictions can be found for rare events, e.g. intermittency [8,[25][26][27][28][29]. This seems to be the case when a more general theoretical framework [30][31][32][33][34] is required, because the nonlinearities are not small [35,36].In this communication, the complete wave-turbulence theory is developed for a fully general 4-wave system, whose hamiltonian is expressed by the following canoni-