Spiking activity in cortical networks is nonlinear in nature. The linear-nonlinear cascade model, some versions of which are also known as point-process generalized linear model, can efficiently capture the nonlinear dynamics exhibited by such networks. Of particular interest in such models are theoretical predictions of spike train statistics. However, due to the moment-closure problem, approximations are inevitable. We suggest here a series expansion that explains how higher-order moments couple to lower-order ones. Our approach makes predictions in terms of certain integrals, the so-called loop integrals. In previous studies these integrals have been evaluated numerically, but numerical instabilities are sometimes encountered rendering the results unreliable. Analytic solutions are presented here to overcome this problem, and to arrive at more robust evaluations. We were able to deduce these analytic solutions by switching to Fourier space and making use of complex analysis, specifically Cauchy's residue theorem. We formalized the loop integrals and explicitly solved them for specific response functions. To quantify the importance of these corrections for spike train cumulants, we numerically simulated spiking networks and compared their sample statistics to our theoretical predictions. Our results demonstrate that the magnitude of the nonlinear corrections depends on the working point of the nonlinear network dynamics, and that it is related to the eigenvalues of the mean-field stability matrix. For our example, the corrections for the firing rates are in the range between 4 % and 21 % on average. Precise and robust predictions of spike train statistics accounting for nonlinear effects are, for example, highly relevant for theories involving spike-timing dependent plasticity (STDP).