The subsonic limit of the quadratic direct interaction approximation (DIA) applied to the Zakharov equations is compared with the cubic DIA applied to the nonlinear Schro¨dinger equation, which is the subsonic limit of the Zakharov equations. Comparisons with Monte Carlo simulations of a truncated system show that the first theory more accurately describes the regime of stationary turbulence, while the second theory more accurately describes the subsonic evolution of the modulational instability. The weak turbulence limits of the two theories describe the sonic and subsonic regimes, respectively. The addition of vertex corrections to the DIA leads to a hybrid weak turbulence theory that smoothly interpolates between the sonic and subsonic regimes.