A maximum-likelihood estimator for trial-to-trial variations in noisy MEG/EEG data sets

The standard procedure to determine the brain response from a multitrial evoked magnetoencephalography (MEG) or electroencephalography (EEG) data set is to average the individual trials of these data, time locked to the stimulus onset. When the brain responses vary from trial-to-trial this approach is false. In this paper, a maximum-likelihood estimator is derived for the case that the recorded data contain amplitude variations. The estimator accounts for spatially and temporally correlated background noise that is superimposed on the brain response. The model is applied to a series of 17 MEG data sets of normal subjects, obtained during median nerve stimulation. It appears that the amplitude of late component (30-120 ms) shows a systematic negative trend indicating a weakening response during stimulation time. For the early components (20-35 ms) no such a systematic effect was found. The model is furthermore applied on a MEG data set consisting of epileptic spikes of constant spatial distribution but varying polarity. For these data, the advantage of applying the model is that positive and negative spikes can be processed with a single model, thereby reducing the number of degrees of freedom and increasing the signal-to-noise ratio.

[1]  J C de Munck,et al.  The use of an MEG device as 3D digitizer and motion monitoring system. , 2001, Physics in medicine and biology.

[2]  Lourens J. Waldorp,et al.  Estimating stationary dipoles from MEG/EEG data contaminated with spatially and temporally correlated background noise , 2001, NeuroImage.

[3]  T. Sejnowski,et al.  Dynamic Brain Sources of Visual Evoked Responses , 2002, Science.

[4]  D Van 't Ent,et al.  Spike cluster analysis in neocortical localization related epilepsy yields clinically significant equivalent source localization results in magnetoencephalogram (MEG) , 2003, Clinical Neurophysiology.

[5]  S. Bressler,et al.  Trial-to-trial variability of cortical evoked responses: implications for the analysis of functional connectivity , 2002, Clinical Neurophysiology.

[6]  F. L. D. Silva,et al.  Event-related EEG/MEG synchronization and desynchronization: basic principles , 1999, Clinical Neurophysiology.

[7]  W. Klimesch.  EEG alpha and theta oscillations reflect cognitive and memory performance: a review and analysis , 1999, Brain Research Reviews.

[8]  P. Jaskowski,et al.  Amplitudes and latencies of single-trial ERP's estimated by a maximum-likelihood method , 1999, IEEE Transactions on Biomedical Engineering.

[9]  T. Gasser,et al.  Variable latencies of noisy signals: Estimation and testing in brain potential data , 1987 .

[10]  J.C. de Munck,et al.  A random dipole model for spontaneous brain activity , 1992, IEEE Transactions on Biomedical Engineering.

[11]  J. Magnus,et al.  Matrix Differential Calculus with Applications in Statistics and Econometrics , 1991 .

[12]  D. G. Wastell,et al.  Potentiation of the habituation of human brain potentials , 1980, Biological Psychology.

[13]  Fetsje Bijma,et al.  A mathematical approach to the temporal stationarity of background noise in MEG/EEG measurements , 2003, NeuroImage.

[14]  Hellmuth Obrig,et al.  Habituation of the Visually Evoked Potential and Its Vascular Response: Implications for Neurovascular Coupling in the Healthy Adult , 2002, NeuroImage.

[15]  B. Pfleiderer,et al.  Visualization of Auditory Habituation by fMRI , 2002, NeuroImage.

[16]  James M Kilner,et al.  Event-related brain dynamics , 2002, Trends in Neurosciences.