Statistics for Nonlinear Surface Plasmon Polaritons: Analytical and Numerical Studies