An exploratory data analysis of electroencephalograms using the functional boxplots approach

Many model-based methods have been developed over the last several decades for analysis of electroencephalograms (EEGs) in order to understand electrical neural data. In this work, we propose to use the functional boxplot (FBP) to analyze log periodograms of EEG time series data in the spectral domain. The functional bloxplot approach produces a median curve—which is not equivalent to connecting medians obtained from frequency-specific boxplots. In addition, this approach identifies a functional median, summarizes variability, and detects potential outliers. By extending FBPs analysis from one-dimensional curves to surfaces, surface boxplots are also used to explore the variation of the spectral power for the alpha (8–12 Hz) and beta (16–32 Hz) frequency bands across the brain cortical surface. By using rank-based nonparametric tests, we also investigate the stationarity of EEG traces across an exam acquired during resting-state by comparing the spectrum during the early vs. late phases of a single resting-state EEG exam.

[1]  Yudi Pawitan,et al.  Nonparametric Spectral Density Estimation Using Penalized Whittle Likelihood , 1994 .

[2]  W. Freeman,et al.  Frequency analysis of olfactory system EEG in cat, rabbit, and rat. , 1980, Electroencephalography and clinical neurophysiology.

[3]  G. Pfurtscheller,et al.  Evaluation of event-related desynchronization (ERD) preceding and following voluntary self-paced movement. , 1979, Electroencephalography and clinical neurophysiology.

[4]  Gustavo Deco,et al.  Functional connectivity dynamics: Modeling the switching behavior of the resting state , 2015, NeuroImage.

[5]  Terrence J. Sejnowski,et al.  Enhanced detection of artifacts in EEG data using higher-order statistics and independent component analysis , 2007, NeuroImage.

[6]  B Fyntanidou,et al.  Bispectral Index Scale variations in patients diagnosed with brain death. , 2012, Transplantation proceedings.

[7]  D. B. Preston.  Spectral Analysis and Time Series , 1983 .

[8]  Ying Sun,et al.  Surface boxplots , 2014, Stat.

[9]  Jianqing Fan,et al.  Automatic Local Smoothing for Spectral Density Estimation , 1998 .

[10]  V. G. Alekseev.  Nonparametric spectral density estimators , 2000 .

[11]  S. Logue,et al.  The neural and genetic basis of executive function: Attention, cognitive flexibility, and response inhibition , 2014, Pharmacology Biochemistry and Behavior.

[12]  M. Genton,et al.  Functional Boxplots , 2011 .

[13]  S. Jain,et al.  Parametric modeling of brain signals , 2004, Proceedings. Eighth International Conference on Information Visualisation.

[14]  D. Brillinger.  Time series - data analysis and theory , 1981, Classics in applied mathematics.

[15]  Tzyy-Ping Jung,et al.  Independent Component Analysis of Electroencephalographic Data , 1995, NIPS.

[16]  Steven C. Cramer,et al.  Resting-state cortical connectivity predicts motor skill acquisition , 2014, NeuroImage.

[17]  M. West,et al.  New methods of time series analysis of non-stationary EEG data: eigenstructure decompositions of time varying autoregressions , 1999, Clinical Neurophysiology.

[18]  Hong-ye Gao.  Choice of thresholds for wavelet shrinkage estimate of the spectrum , 1997 .

[19]  Richard A. Davis,et al.  Time Series: Theory and Methods , 2013 .

[20]  David S. Stoffer,et al.  Time series analysis and its applications , 2000 .

[21]  S. Makeig.  Auditory event-related dynamics of the EEG spectrum and effects of exposure to tones. , 1993, Electroencephalography and clinical neurophysiology.

[22]  Kristin Potter,et al.  Surface boxplots. , 2014, Stat.

[23]  J. Fermaglich.  Electric Fields of the Brain: The Neurophysics of EEG , 1982 .

[24]  A. Engel,et al.  Beta-band oscillations—signalling the status quo? , 2010, Current Opinion in Neurobiology.

[25]  J. Romo,et al.  On the Concept of Depth for Functional Data , 2009 .

[26]  A. Wennberg,et al.  Computer analysis of EEG signals with parametric models , 1981, Proceedings of the IEEE.

[27]  M. V. van Putten,et al.  Continuous EEG Monitoring During Thrombolysis in Acute Hemispheric Stroke Patients Using the Brain Symmetry Index , 2008, Journal of clinical neurophysiology : official publication of the American Electroencephalographic Society.

[28]  Thomas C. M. Lee.  A simple span selector for periodogram smoothing , 1997 .

[29]  P. Nunez,et al.  Electric fields of the brain , 1981 .

[30]  Catie Chang,et al.  Time–frequency dynamics of resting-state brain connectivity measured with fMRI , 2010, NeuroImage.

[31]  G. Wahba.  Automatic Smoothing of the Log Periodogram , 1980 .

[32]  Miles A Whittington,et al.  A beta2-frequency (20–30 Hz) oscillation in nonsynaptic networks of somatosensory cortex , 2006, Proceedings of the National Academy of Sciences.

[33]  Richard M. Leahy,et al.  General Linear Modeling of Magnetoencephalography Data , 2012 .

[34]  T. Bohlin.  Comparison of two methods of modeling stationary EEG signals , 1973 .

[35]  Marc G. Genton,et al.  Adjusted functional boxplots for spatio‐temporal data visualization and outlier detection , 2012 .

[36]  Tipu Z. Aziz,et al.  Driving Oscillatory Activity in the Human Cortex Enhances Motor Performance , 2012, Current Biology.

[37]  J. Raz,et al.  A simple generalised crossvalidation method of span selection for periodogram smoothing , 2001 .

[38]  Marc G. Genton,et al.  Functional Median Polish , 2012, Journal of Agricultural, Biological, and Environmental Statistics.

[39]  D. Nychka,et al.  Exact fast computation of band depth for large functional datasets: How quickly can one million curves be ranked? , 2012 .

[40]  Eric R. Ziegel,et al.  Generalized Linear Models , 2002, Technometrics.

[41]  Donald B. Percival,et al.  Spectral Analysis for Physical Applications , 1993 .