Nonintegrable systems thermalize, leading to the emergence of fluctuating hydrodynamics. Typically, this hydrodynamics is diffusive. We use the effective field theory (EFT) of diffusion to compute higher-point functions of conserved densities. We uncover a simple scaling behavior of correlators at late times, and, focusing on three and four-point functions, derive the asymptotically exact universal scaling functions that characterize nonlinear response in diffusive systems. This allows for precision tests of thermalization beyond linear response in quantum and classical many-body systems. We confirm our predictions in a classical lattice gas.