High-order statistics for the random paraxial wave equation