Interactions among neural signals in different frequency components have become a focus of strong interest in biomedical signal processing. The bispectrum is a method to detect the presence of quadratic phase coupling (QPC) between different frequency bands in a signal. The traditional way to quantify phase coupling is by means of the bicoherence index (BCI), which is essentially a normalized bispectrum. The main disadvantage of the BCI is that the determination of significant QPC becomes compromised with noise. To mitigate this problem, a statistical approach that combines the bispectrum with an improved surrogate data method to determine the statistical significance of the phase coupling is introduced. The method was first tested on two simulation examples. It was then applied to the human EEG signal that has been recorded from the scalp using international 10–20 electrodes system. The frequency domain method, based on normalized spectrum and bispectrum, describes frequency interactions associated with nonlinearities occurring in the observed EEG.