Non-Gaussian probabilistic MEG source localisation based on kernel density estimation.

There is strong evidence to suggest that data recorded from magnetoencephalography (MEG) follows a non-Gaussian distribution. However, existing standard methods for source localisation model the data using only second order statistics, and therefore use the inherent assumption of a Gaussian distribution. In this paper, we present a new general method for non-Gaussian source estimation of stationary signals for localising brain activity from MEG data. By providing a Bayesian formulation for MEG source localisation, we show that the source probability density function (pdf), which is not necessarily Gaussian, can be estimated using multivariate kernel density estimators. In the case of Gaussian data, the solution of the method is equivalent to that of widely used linearly constrained minimum variance (LCMV) beamformer. The method is also extended to handle data with highly correlated sources using the marginal distribution of the estimated joint distribution, which, in the case of Gaussian measurements, corresponds to the null-beamformer. The proposed non-Gaussian source localisation approach is shown to give better spatial estimates than the LCMV beamformer, both in simulations incorporating non-Gaussian signals, and in real MEG measurements of auditory and visual evoked responses, where the highly correlated sources are known to be difficult to estimate.

MEG can map short and long-term changes in brain activity following deep brain stimulation for chronic pain.

Deep brain stimulation (DBS) has been shown to be clinically effective for some forms of treatment-resistant chronic pain, but the precise mechanisms of action are not well understood. Here, we present an analysis of magnetoencephalography (MEG) data from a patient with whole-body chronic pain, in order to investigate changes in neural activity induced by DBS for pain relief over both short- and long-term. This patient is one of the few cases treated using DBS of the anterior cingulate cortex (ACC). We demonstrate that a novel method, null-beamforming, can be used to localise accurately brain activity despite the artefacts caused by the presence of DBS electrodes and stimulus pulses. The accuracy of our source localisation was verified by correlating the predicted DBS electrode positions with their actual positions. Using this beamforming method, we examined changes in whole-brain activity comparing pain relief achieved with deep brain stimulation (DBS ON) and compared with pain experienced with no stimulation (DBS OFF). We found significant changes in activity in pain-related regions including the pre-supplementary motor area, brainstem (periaqueductal gray) and dissociable parts of caudal and rostral ACC. In particular, when the patient reported experiencing pain, there was increased activity in different regions of ACC compared to when he experienced pain relief. We were also able to demonstrate long-term functional brain changes as a result of continuous DBS over one year, leading to specific changes in the activity in dissociable regions of caudal and rostral ACC. These results broaden our understanding of the underlying mechanisms of DBS in the human brain.

Fusion of magnetometer and gradiometer sensors of MEG in the presence of multiplicative error.

Novel neuroimaging techniques have provided unprecedented information on the structure and function of the living human brain. Multimodal fusion of data from different sensors promises to radically improve this understanding, yet optimal methods have not been developed. Here, we demonstrate a novel method for combining multichannel signals. We show how this method can be used to fuse signals from the magnetometer and gradiometer sensors used in magnetoencephalography (MEG), and through extensive experiments using simulation, head phantom and real MEG data, show that it is both robust and accurate. This new approach works by assuming that the lead fields have multiplicative error. The criterion to estimate the error is given within a spatial filter framework such that the estimated power is minimized in the worst case scenario. The method is compared to, and found better than, existing approaches. The closed-form solution and the conditions under which the multiplicative error can be optimally estimated are provided. This novel approach can also be employed for multimodal fusion of other multichannel signals such as MEG and EEG. Although the multiplicative error is estimated based on beamforming, other methods for source analysis can equally be used after the lead-field modification.

Localizing heart sounds in respiratory signals using singular spectrum analysis.

Respiratory sounds are always contaminated by heart sound interference. An essential preprocessing step in some of the heart sound cancellation methods is localizing primary heart sound components. Singular spectrum analysis (SSA), a powerful time series analysis technique, is used in this paper. Despite the frequency overlap of the heart and lung sound components, two different trends in the eigenvalue spectra are recognizable, which leads to find a subspace that contains more information about the underlying heart sound. Artificially mixed and real respiratory signals are used for evaluating the performance of the method. Selecting the appropriate length for the SSA window results in good decomposition quality and low computational cost for the algorithm. The results of the proposed method are compared with those of well-established methods, which use the wavelet transform and entropy of the signal to detect the heart sound components. The proposed method outperforms the wavelet-based method in terms of false detection and also correlation with the underlying heart sounds. Performance of the proposed method is slightly better than that of the entropy-based method. Moreover, the execution time of the former is significantly lower than that of the latter.

Application of a null-beamformer to source localisation in MEG data of deep brain stimulation.

In this paper, we present an analysis of magnetoencephalography (MEG) signals from a patient with whole-body chronic pain in order to investigate changes in neural activity induced by DBS. The patient is one of the few cases treated using DBS of the anterior cingulate cortex (ACC). Using MEG to reconstruct the neural activity of interest is challenging because of interference to the signal from the DBS device. We demonstrate that a null-beamformer can be used to localise neural activity despite artefacts caused by the presence of DBS electrodes and stimulus pulses. We subsequently verified the accuracy of our source localisation by correlating the predicted DBS electrode positions with their actual positions, previously identified using anatomical imaging. We also demonstrated increased activity in pain-related regions including the pre-supplementary motor area, brainstem periaqueductal gray and medial prefrontal areas when the patient was in pain compared to when the patient experienced pain relief.

A new method for spatiotemporal identification of event-related potential subcomponents.

In this study a novel method for tracking and separation of event-related potential (ERP) subcomponents from trial to trial is considered. The sources of ERP subcomponents are assumed to be electric current dipoles (ECD). The shape of each ERP subcomponent is also supposed to be monophasic wave and modeled using a Gaussian waveform. We are interested in the estimation and tracking of ERP subcomponent locations and parameters (amplitude, latency and width of each Gaussian waveform). Estimation of ECD locations, which have nonlinear relation to the measurement, is performed by particle filtering, and estimation of the amplitude is optimally estimated by a maximum likelihood approach, and finally estimation of latency and width of the Gaussian functions are given by Newton-Raphson technique. New recursive methods are introduced for both maximum likelihood and Newton-Raphson approaches to prevent the divergence of the filtering in the presence of very low signal to noise ratio (SNR). The proposed method was assessed using both simulated and real data and the results verified a successful deployment of the method in ERP analysis.

The application of particle filters in single trial event-related potential estimation.

In this paper, an approach for the estimation of single trial event-related potentials (ST-ERPs) using particle filters (PFs) is presented. The method is based on recursive Bayesian mean square estimation of ERP wavelet coefficients using their previous estimates as prior information. To enable a performance evaluation of the approach in the Gaussian and non-Gaussian distributed noise conditions, we added Gaussian white noise (GWN) and real electroencephalogram (EEG) signals recorded during rest to the simulated ERPs. The results were compared to that of the Kalman filtering (KF) approach demonstrating the robustness of the PF over the KF to the added GWN noise. The proposed method also outperforms the KF when the assumption about the Gaussianity of the noise is violated. We also applied this technique to real EEG potentials recorded in an odd-ball paradigm and investigated the correlation between the amplitude and the latency of the estimated ERP components. Unlike the KF method, for the PF there was a statistically significant negative correlation between amplitude and latency of the estimated ERPs, matching previous neurophysiological findings.

Detection of rhythmic discharges in newborn EEG signals.

This paper presents a scalp eletroencephalogram (EEG) rhythmic pattern detection scheme based on neural networks. rhythmic discharges detection is applicable to the majority of seizures seen in newborns, and is listed as detecting 90% of all the seizures. In this approach some features based on various methods are extracted and compared by a modified multilayer neural network in order to find rhythmic discharges. Statistical performance comparison with seizure detection schemes of Gotman et al. and Liu et al. is performed.

Seizure detection in EEG signals: a comparison of different approaches.

In this paper, the performance of traditional variance-based method for detection of epileptic seizures in EEG signals are compared with various methods based on nonlinear time series analysis, entropies, logistic regression,discrete wavelet transform and time frequency distributions.We noted that variance-based method in compare to the mentioned methods had the best result (100%) applied on the same database.