A Subspace Method for Dynamical Estimation of Evoked Potentials

It is a challenge in evoked potential (EP) analysis to incorporate prior physiological knowledge for estimation. In this paper, we address the problem of single-channel trial-to-trial EP characteristics estimation. Prior information about phase-locked properties of the EPs is assesed by means of estimated signal subspace and eigenvalue decomposition. Then for those situations that dynamic fluctuations from stimulus-to-stimulus could be expected, prior information can be exploited by means of state-space modeling and recursive Bayesian mean square estimation methods (Kalman filtering and smoothing). We demonstrate that a few dominant eigenvectors of the data correlation matrix are able to model trend-like changes of some component of the EPs, and that Kalman smoother algorithm is to be preferred in terms of better tracking capabilities and mean square error reduction. We also demonstrate the effect of strong artifacts, particularly eye blinks, on the quality of the signal subspace and EP estimates by means of independent component analysis applied as a prepossessing step on the multichannel measurements

Identification of audio evoked response potentials in ambulatory EEG data

Electroencephalography (EEG) is commonly used for observing brain function over a period of time. It employs a set of invasive electrodes on the scalp to measure the electrical activity of the brain. EEG is mainly used by researchers and clinicians to study the brain’s responses to a specific stimulus - the event-related potentials (ERPs). Different types of undesirable signals, which are known as artefacts, contaminate the EEG signal. EEG and ERP signals are very small (in the order of microvolts); they are often obscured by artefacts with much larger amplitudes in the order of millivolts. This greatly increases the difficulty of interpreting EEG and ERP signals.Typically, ERPs are observed by averaging EEG measurements made with many repetitions of the stimulus. The average may require many tens of repetitions before the ERP signal can be observed with any confidence. This greatly limits the study and useof ERPs. This project explores more sophisticated methods of ERP estimation from measured EEGs. An Optimal Weighted Mean (OWM) method is developed that forms a weighted average to maximise the signal to noise ratio in the mean. This is developedfurther into a Bayesian Optimal Combining (BOC) method where the information in repetitions of ERP measures is combined to provide a sequence of ERP estimations with monotonically decreasing uncertainty. A Principal Component Analysis (PCA) isperformed to identify the basis of signals that explains the greatest amount of ERP variation. Projecting measured EEG signals onto this basis greatly reduces the noise in measured ERPs. The PCA filtering can be followed by OWM or BOC. Finally, crosschannel information can be used. The ERP signal is measured on many electrodes simultaneously and an improved estimate can be formed by combining electrode measurements. A MAP estimate, phrased in terms of Kalman Filtering, is developed using all electrode measurements.The methods developed in this project have been evaluated using both synthetic and measured EEG data. A synthetic, multi-channel ERP simulator has been developed specifically for this project.Numerical experiments on synthetic ERP data showed that Bayesian Optimal Combining of trial data filtered using a combination of PCA projection and Kalman Filtering, yielded the best estimates of the underlying ERP signal. This method has been applied to subsets of real Ambulatory Electroencephalography (AEEG) data, recorded while participants performed a range of activities in different environments. From this analysis, the number of trials that need to be collected to observe the P300 amplitude and delay has been calculated for a range of scenarios

Detecting single-trial EEG evoked potential using a wavelet domain linear mixed model: application to error potentials classification

Objective. The main goal of this work is to develop a model for multi-sensor signals such as MEG or EEG signals, that accounts for the inter-trial variability, suitable for corresponding binary classification problems. An important constraint is that the model be simple enough to handle small size and unbalanced datasets, as often encountered in BCI type experiments. Approach. The method involves linear mixed effects statistical model, wavelet transform and spatial filtering, and aims at the characterization of localized discriminant features in multi-sensor signals. After discrete wavelet transform and spatial filtering, a projection onto the relevant wavelet and spatial channels subspaces is used for dimension reduction. The projected signals are then decomposed as the sum of a signal of interest (i.e. discriminant) and background noise, using a very simple Gaussian linear mixed model. Main results. Thanks to the simplicity of the model, the corresponding parameter estimation problem is simplified. Robust estimates of class-covariance matrices are obtained from small sample sizes and an effective Bayes plug-in classifier is derived. The approach is applied to the detection of error potentials in multichannel EEG data, in a very unbalanced situation (detection of rare events). Classification results prove the relevance of the proposed approach in such a context. Significance. The combination of linear mixed model, wavelet transform and spatial filtering for EEG classification is, to the best of our knowledge, an original approach, which is proven to be effective. This paper improves on earlier results on similar problems, and the three main ingredients all play an important role

Multichannel EEG : towards applications in clinical neurology.

Electroencephalogram (EEG) measures the electric activity produced by the brain with electrodes placed on the scalp. It is used for monitoring or as diagnostic tool for neurological disorders. In practice a maximum of 21 electrodes are generally used for a clinical EEG recording. However, EEG systems with 128 and 256 electrodes are also available and used for fundamental research. In this thesis we investigate whether the extra information obtained with 128-channel recordings is clinically relevant. We have focused on evoked potentials (EPs). EP is the electric activity of the brain caused by a stimulus (e.g. a flashlight). We showed that a measure often used for evoked potentials, the peak amplitude, can be estimated more accurately by using 128 channels recordings than by conventional recordings. Therefore this technique might be more sensitive to pathological changes. In addition, we developed a new technique to estimate EP symmetry (similarity of EPs generated in left and right hemisphere). This technique might be useful for diagnosis of neurological disorders with brain damage in one hemisphere. Both methods have been applied to a group of patients with parkinsonism; neurological symptoms typical for Parkinson’s disease. No differences could be observed in amplitude or symmetry between patients with different parkinsonian disorders. Therefore, (so far) these methods cannot be used as diagnostic tool for neurological disorders. Future research will show whether small adaptations to the stimulation method or analysis technique will result in an improvement of the diagnostic value and whether these methods are useful for other neurological disorders.

Single-Trial Extraction of Pure Somatosensory Evoked Potential Based on Expectation Maximization Approach

published_or_final_versio

Single Epoch Analysis and Bi-hemisphere Study of Magnetoencephalographic (MEG) Signals using Vector Signal Transformation V3 and Magnetic Field Tomography (MFT)

The biomagnetic inverse problem has no unique solution, nevertheless even a cursory look at the features shown in raw signal can often suffice to highlight strong superficial activity. To do a proper single epoch analysis is normally prohibitively expensive in terms of computing demands. Hence the original aim of this thesis was to use simple efficient signal transformations to characterize superficial generators and contrast the single epoch signature with that extracted from the average signal. The results have intrigued us sufficiently to go beyond the original goal and extract very preliminary estimates of activity across the cerebral hemisphere in single trials. The original tool, and one that we have used for much of the work, is a simple vector signal transformation called V3. This signal transformation highlights nearby sources; it is a crude but quick estimator of generators directly from the raw MEG signals. Together with Magnetic Field Tomography (MFT), which relies on distributed source analysis of the MEG signals, we have tackled the following specific problems relating to aspects of normal brain function: efficient estimation of generators of magnetic fields; relationship between the average signal and single trials; and interhemispheric differences and relationship between the activity in the left and right hemispheres of the brain. During the project, we have used as examples auditory evoked MEG measurements obtained from two multichannel systems and applied the V3 and MFT analysis to both the average and single trial signals. In particular, we chose the 40-Hz (or gamma band) auditory response as the study subject. We found that in single epochs similar patterns of high frequency activity are observed in the area around the auditory cortex well before, close to and well after stimulus onset; the sequence of events observed in the average can only represent the evolution of events in single trials in a statistical way; and deep and central areas of the brain may be the seeds for the main deflections observed in the auditory responses

Analysis of Dynamic Brain Imaging Data

Modern imaging techniques for probing brain function, including functional Magnetic Resonance Imaging, intrinsic and extrinsic contrast optical imaging, and magnetoencephalography, generate large data sets with complex content. In this paper we develop appropriate techniques of analysis and visualization of such imaging data, in order to separate the signal from the noise, as well as to characterize the signal. The techniques developed fall into the general category of multivariate time series analysis, and in particular we extensively use the multitaper framework of spectral analysis. We develop specific protocols for the analysis of fMRI, optical imaging and MEG data, and illustrate the techniques by applications to real data sets generated by these imaging modalities. In general, the analysis protocols involve two distinct stages: `noise' characterization and suppression, and `signal' characterization and visualization. An important general conclusion of our study is the utility of a frequency-based representation, with short, moving analysis windows to account for non-stationarity in the data. Of particular note are (a) the development of a decomposition technique (`space-frequency singular value decomposition') that is shown to be a useful means of characterizing the image data, and (b) the development of an algorithm, based on multitaper methods, for the removal of approximately periodic physiological artifacts arising from cardiac and respiratory sources.Comment: 40 pages; 26 figures with subparts including 3 figures as .gif files. Originally submitted to the neuro-sys archive which was never publicly announced (was 9804003