Removal of non-white noise from single trial event-related EEG signals using soft-thresholding

Infrequent stimulation of a subject generates event-related potential (ERP) signals masked by background EEG activity. It is generally assumed that this background activity is normal white noise. Another assumption made is that the underlying evoked signals are deterministic and they do not vary from trial to trial. The signal, thus, is modeled as y/sub i//sup (j)/=x/sub i/+/spl sigma//sub i//sup (j)/, where y/sub i//sup (j)/ is the jth trial and z/sub i/ is a white noise. The signal x/sub i/ is recovered by averaging the observations y/sub i//sup (j)/. In reality, the background activity is "colored" and not always Gaussian. The time samples of the background activity are generally correlated. In addition, the signal x/sub i/ varies across observations. The purpose of this study is to extract single trial ERPs from the EEG. The authors are interested in the trial to trial variation of the ERPs and their clinical applicability. They have investigated the wavelet soft-thresholding method to remove the background noise and separate out the single trial response. The non-white background activity, after wavelet transformation, may concentrate in certain resolution levels. The authors determined these levels by testing the Gaussianity of the wavelet coefficients, using both the /spl chi//sup 2/ and Kolmogorov-Smirnov goodness of fit tests. In the resolution levels at which the null hypothesis was not rejected, the noise level was estimated. The de-noising threshold was then calculated using a level dependent rule T/sub j,N/=/spl radic/(2 log(N))/spl middot/MAD(C/sub j,k/)/0.6745, where C/sub j,k/ are the wavelet coefficients and MAD(C/sub j,k/)=Median(|C/sub j,k|/) is an estimator of the noise level. The resolution levels at which the null hypothesis was rejected were not thresholded.