Statistics for A higher-order perturbation analysis of the nonlinear Schrödinger equation