Maximum-likelihood estimation of low-rank signals for multiepoch MEG/EEG analysis

A maximum-likelihood-based algorithm is presented for reducing the effects of spatially colored noise in evoked response magneto- and electro-encephalography data. The repeated component of the data, or signal of interest, is modeled as the mean, while the noise is modeled as the Kronecker product of a spatial and a temporal covariance matrix. The temporal covariance matrix is assumed known or estimated prior to the application of the algorithm. The spatial covariance structure is estimated as part of the maximum-likelihood procedure. The mean matrix representing the signal of interest is assumed to be low-rank due to the temporal and spatial structure of the data. The maximum-likelihood estimates of the components of the lowrank signal structure are derived in order to estimate the signal component. The relationship between this approach and principal component analysis (PCA) is explored. In contrast to prestimulus-based whitening followed by PCA, the maximumlikelihood approach does not require signal-free data for noise whitening. Consequently, the maximum-likelihood approach is much more effective with nonstationary noise and produces better quality whitening for a given data record length. The efficacy of this approach is demonstrated using simulated and real MEG data.