Fast and Scalable Multiway Analysis of Neural Data

Analysis of neural data with multiple modes and high density has recently become a trend with the advances in neuroscience research and practices. There exists a pressing need for an approach to accurately and uniquely capture the features without loss or destruction of the interactions amongst the modes (typically) of space, time, and frequency. Moreover, the approach must be able to quickly analyze the neural data of exponentially growing scales and sizes, in tens or even hundreds of channels, so that timely conclusions and decisions may be made. A salient approach to multi-way data analysis is the parallel factor analysis (PARAFAC) that manifests its effectiveness in the decomposition of the electroencephalography (EEG). However, the conventional PARAFAC is only suited for offline data analysis due to the high complexity, which computes to be O(n) with the increasing data size. In this study, a large-scale PARAFAC method has been developed, which is supported by general-purpose computing on the graphics processing unit (GPGPU). Comparing to the PARAFAC running on conventional CPU-based platform, the new approach dramatically excels by >360 times in run-time performance, and effectively scales by >400 times in all dimensions. Moreover, the proposed approach forms the basis of a model for the analysis of electrocochleography (ECoG) recordings obtained from epilepsy patients, which proves to be effective in the epilepsy state detection. The time evolutions of the proposed model are well correlated with the clinical observations. Moreover, the frequency signature is stable and high in the ictal phase. Furthermore, the spatial signature explicitly identifies the propagation of neural activities among various brain regions. The model supports real-time analysis of ECoG in >1000 channels on an inexpensive and available cyber-infrastructure.

[1]  Jean Gotman,et al.  Removing Eye‐movement Artifacts from the EEG during the Intracarotid Amobarbital Procedure , 2005, Epilepsia.

[2]  H. Law Research methods for multimode data analysis , 1984 .

[3]  Jack Dongarra,et al.  An Improved MAGMA GEMM for Fermi GPUs , 2010 .

[4]  Fumikazu Miwakeichi,et al.  Decomposing EEG data into space–time–frequency components using Parallel Factor Analysis , 2004, NeuroImage.

[5]  R. Bro PARAFAC. Tutorial and applications , 1997 .

[6]  Patrick Dupont,et al.  Canonical decomposition of ictal scalp EEG reliably detects the seizure onset zone , 2007, NeuroImage.

[7]  R. Bro,et al.  A new efficient method for determining the number of components in PARAFAC models , 2003 .

[8]  Didier G. Leibovici,et al.  Multi-way modelling of high-dimensionality electroencephalographic data , 2001 .

[9]  John D. Owens,et al.  GPU Computing , 2008, Proceedings of the IEEE.

[10]  G. Rabinovici,et al.  Distinct clinical and metabolic deficits in PCA and AD are not related to amyloid distribution , 2011, Neurology.

[11]  Konrad P Kording,et al.  How advances in neural recording affect data analysis , 2011, Nature Neuroscience.

[12]  G. Fox,et al.  High Performance Multidimensional Scaling for Large High-Dimensional Data Visualization , 2012 .

[13]  B. Kowalski,et al.  Tensorial resolution: A direct trilinear decomposition , 1990 .

[14]  Numaan Ashraf CUDA accelerated face recognition , 2010 .

[15]  R. Harshman,et al.  PARAFAC: parallel factor analysis , 1994 .

[16]  Rasmus Bro,et al.  Multiway analysis of epilepsy tensors , 2007, ISMB/ECCB.

[17]  Karl Pearson F.R.S. LIII. On lines and planes of closest fit to systems of points in space , 1901 .

[18]  Andrzej Cichocki,et al.  PARAFAC algorithms for large-scale problems , 2011, Neurocomputing.

[19]  P. J. Narayanan,et al.  Singular value decomposition on GPU using CUDA , 2009, 2009 IEEE International Symposium on Parallel & Distributed Processing.

[20]  Emmanuel Agullo,et al.  QR Factorization on a Multicore Node Enhanced with Multiple GPU Accelerators , 2011, 2011 IEEE International Parallel & Distributed Processing Symposium.

[21]  Rasmus Bro,et al.  Improving the speed of multi-way algorithms:: Part I. Tucker3 , 1998 .

[22]  Andrzej Cichocki,et al.  Multiway array decomposition analysis of EEGs in Alzheimer's disease , 2012, Journal of Neuroscience Methods.

[23]  Lars Kai Hansen,et al.  Parallel Factor Analysis as an exploratory tool for wavelet transformed event-related EEG , 2006, NeuroImage.

[24]  Dan Brandt,et al.  Investigation of GPGPU for use in processing of EEG in real-time , 2010 .

[25]  Chengbiao Lu,et al.  Characteristics of Evoked Potential Multiple EEG Recordings in Patients with Chronic Pain by Means of Parallel Factor Analysis , 2012, Comput. Math. Methods Medicine.

[26]  J. G. F. Francis,et al.  The QR Transformation - Part 2 , 1962, Comput. J..

[27]  Gerrit Kateman,et al.  Generalized rank annihilation method. I: Derivation of eigenvalue problems , 1994 .

[28]  Pierre Comon,et al.  Independent component analysis, A new concept? , 1994, Signal Process..

[29]  Tamara G. Kolda,et al.  Tensor Decompositions and Applications , 2009, SIAM Rev..

[30]  J. Chang,et al.  Analysis of individual differences in multidimensional scaling via an n-way generalization of “Eckart-Young” decomposition , 1970 .

[31]  Hong Bao,et al.  GPGPU-Aided Ensemble Empirical-Mode Decomposition for EEG Analysis During Anesthesia , 2010, IEEE Transactions on Information Technology in Biomedicine.

[32]  Xiaoli Li,et al.  Massively parallel non-stationary EEG data processing on GPGPU platforms with Morlet continuous wavelet transform , 2012, Journal of Internet Services and Applications.

[33]  D. Burdick An introduction to tensor products with applications to multiway data analysis , 1995 .