Limitations in the rapid extraction of evoked potentials using parametric modeling

A multi-dimensional functional principal components analysis of EEG data

The electroencephalography (EEG) data created in event-related potential (ERP) experiments have a complex high-dimensional structure. Each stimulus presentation, or trial, generates an ERP waveform which is an instance of functional data. The experiments are made up of sequences of multiple trials, resulting in longitudinal functional data and moreover, responses are recorded at multiple electrodes on the scalp, adding an electrode dimension. Traditional EEG analyses involve multiple simplifications of this structure to increase the signal-to-noise ratio, effectively collapsing the functional and longitudinal components by identifying key features of the ERPs and averaging them across trials. Motivated by an implicit learning paradigm used in autism research in which the functional, longitudinal, and electrode components all have critical interpretations, we propose a multidimensional functional principal components analysis (MD-FPCA) technique which does not collapse any of the dimensions of the ERP data. The proposed decomposition is based on separation of the total variation into subject and subunit level variation which are further decomposed in a two-stage functional principal components analysis. The proposed methodology is shown to be useful for modeling longitudinal trends in the ERP functions, leading to novel insights into the learning patterns of children with Autism Spectrum Disorder (ASD) and their typically developing peers as well as comparisons between the two groups. Finite sample properties of MD-FPCA are further studied via extensive simulations.

A MISO-ARX-Based Method for Single-Trial Evoked Potential Extraction

In this paper, we propose a novel method for solving the single-trial evoked potential (EP) estimation problem. In this method, the single-trial EP is considered as a complex containing many components, which may originate from different functional brain sites; these components can be distinguished according to their respective latencies and amplitudes and are extracted simultaneously by multiple-input single-output autoregressive modeling with exogenous input (MISO-ARX). The extraction process is performed in three stages: first, we use a reference EP as a template and decompose it into a set of components, which serve as subtemplates for the remaining steps. Then, a dictionary is constructed with these subtemplates, and EPs are preliminarily extracted by sparse coding in order to roughly estimate the latency of each component. Finally, the single-trial measurement is parametrically modeled by MISO-ARX while characterizing spontaneous electroencephalographic activity as an autoregression model driven by white noise and with each component of the EP modeled by autoregressive-moving-average filtering of the subtemplates. Once optimized, all components of the EP can be extracted. Compared with ARX, our method has greater tracking capabilities of specific components of the EP complex as each component is modeled individually in MISO-ARX. We provide exhaustive experimental results to show the effectiveness and feasibility of our method.

Single-Trial Sparse Representation-Based Approach for VEP Extraction

Sparse representation is a powerful tool in signal denoising, and visual evoked potentials (VEPs) have been proven to have strong sparsity over an appropriate dictionary. Inspired by this idea, we present in this paper a novel sparse representation-based approach to solving the VEP extraction problem. The extraction process is performed in three stages. First, instead of using the mixed signals containing the electroencephalogram (EEG) and VEPs, we utilise an EEG from a previous trial, which did not contain VEPs, to identify the parameters of the EEG autoregressive (AR) model. Second, instead of the moving average (MA) model, sparse representation is used to model the VEPs in the autoregressive-moving average (ARMA) model. Finally, we calculate the sparse coefficients and derive VEPs by using the AR model. Next, we tested the performance of the proposed algorithm with synthetic and real data, after which we compared the results with that of an AR model with exogenous input modelling and a mixed overcomplete dictionary-based sparse component decomposition method. Utilising the synthetic data, the algorithms are then employed to estimate the latencies of P100 of the VEPs corrupted by added simulated EEG at different signal-to-noise ratio (SNR) values. The validations demonstrate that our method can well preserve the details of the VEPs for latency estimation, even in low SNR environments.

Identifying longitudinal trends within EEG experiments

Differential brain response to sensory stimuli is very small (a few microvolts) compared to the overall magnitude of spontaneous electroencephalogram (EEG), yielding a low signal-to-noise ratio (SNR) in studies of event-related potentials (ERP). To cope with this phenomenon, stimuli are applied repeatedly and the ERP signals arising from the individual trials are averaged at the subject level. This results in loss of information about potentially important changes in the magnitude and form of ERP signals over the course of the experiment. In this article, we develop a meta-preprocessing step utilizing a moving average of ERP across sliding trial windows, to capture such longitudinal trends. We embed this procedure in a weighted linear mixed effects model to describe longitudinal trends in features such as ERP peak amplitude and latency across trials while adjusting for the inherent heteroskedasticity created at the meta-preprocessing step. The proposed unified framework, including the meta-processing and the weighted linear mixed effects modeling steps, is referred to as MAP-ERP (moving-averaged-processed ERP). We perform simulation studies to assess the performance of MAP-ERP in reconstructing existing longitudinal trends and apply MAP-ERP to data from young children with autism spectrum disorder (ASD) and their typically developing counterparts to examine differences in patterns of implicit learning, providing novel insights about the mechanisms underlying social and/or cognitive deficits in this disorder.

Abnormal dynamics of EEG oscillations in schizophrenia patients on multiple time scales

Neuronal oscillations reflect the activity of neuronal ensembles engaged in integrative cognition, and may serve as a functional measure for the cognitive impairment in schizophrenia. This study aims to reveal the abnormal amplitude dynamics of electroencephalogram (EEG) oscillations in schizophrenia patients on multiple time scales. EEGs were recorded from schizophrenia patients ( n = 19) and healthy controls ( n = 16) while they were at resting state with eyes closed, at resting state with eyes open, and at watching video. Detrended fluctuation analysis and measures of life-time and waiting-time were used to characterize the abnormal dynamics of EEG oscillations on both long (1-20 s) and short (≤1 s) time scales. Abnormal dynamics of EEG oscillations in alpha and beta bands were observed. In particular, compared with healthy controls, schizophrenia patients have smaller DFA exponent (implying weaker long-range temporal correlation) in the left fronto-temporal area and smaller DFA exponent, smaller life-time (indicating shorter oscillation burst), and smaller waiting-time in the occipital area in beta band at resting state with eyes open. In addition, schizophrenia patients have larger DFA exponent, larger life-time, and larger waiting-time at some clustered channels in the temporo-parietal area in alpha band at watching video. The present results provide new insights for cognitive deficits and the underlying neuronal dysfunction in schizophrenia.

Noninvasive imaging of the high frequency brain activity in focal epilepsy patients

High-frequency (HF) activity represents a potential biomarker of the epileptogenic zone in epilepsy patients, the removal of which is considered to be crucial for seizure-free surgical outcome. We proposed a high frequency source imaging (HFSI) approach to noninvasively image the brain sources of the scalp-recorded HF EEG activity. Both computer simulation and clinical patient data analysis were performed to investigate the feasibility of using the HFSI approach to image the sources of HF activity from noninvasive scalp EEG recordings. The HF activity was identified from high-density scalp recordings after high-pass filtering the EEG data and the EEG segments with HF activity were concatenated together to form repetitive HF activity. Independent component analysis was utilized to extract the components corresponding to the HF activity. Noninvasive EEG source imaging using realistic geometric boundary element head modeling was then applied to image the sources of the pathological HF brain activity. Five medically intractable focal epilepsy patients were studied and the estimated sources were found to be concordant with the surgical resection or intracranial recordings of the patients. The present study demonstrates, for the first time, that source imaging from the scalp HF activity could help to localize the seizure onset zone and provide a novel noninvasive way of studying the epileptic brain in humans. This study also indicates the potential application of studying HF activity in the presurgical planning of medically intractable epilepsy patients.