Stepwise model order estimation in blind source separation applied to ictal EEG

Most algorithms for blind source separation (BSS) or independent component analysis (ICA) assume an equal number of sources as sensors. For multichannel electrophysiological recordings, such as the electroencephalogram (EEG), however, there are often far fewer sources of neurophysiologically relevant activity than the number of sensors. This adds a model order estimation problem to the source separation problem. Conventional estimates of the number of sources are based on the dominant eigenvalues of the data covariance matrix, obtained from principal component analysis (PCA), whose corresponding eigenvectors are also used for prewhitening. It is well known that PCA is susceptible to noise, leading to incorrect model order estimates and data distortion, which in turn limit the accuracy of the source estimates. It is therefore highly desirable to determine the correct number of sources and their spatial topographies directly, without PCA-based data truncation or prewhitening. In this work, we present a stepwise BSS method for extracting only the sources necessary for a sufficiently good least-square fit to the data. This simultaneously yields model order and source estimates, which we examine at different noise levels. We also show how only a few neurophysiologically meaningful components can be extracted from 25-channel ictal EEG.