Blind separation and localization of correlated P300 subcomponents from single trial recordings using extended PARAFAC2 tensor model

A novel mathematical model based on multi-way data construction and analysis with the goal of simultaneously separating and localizing the brain sources specially the subcomponents of event related potentials (ERPs) is introduced. We represent multi-channel EEG data using a third-order tensor with modes: space (channels), time samples, and number of segments. Then, a multi-way technique, in particular, generalized version of PARAFAC2 method, is developed to blindly separate and localize mutually/temporally correlated P3a and P3b sources as subcomponents of P300 signal. In this paper the non-orthogonality of the ERP subcomponents is defined within the tensor model. In order to obtain essentially unique estimation of the signal components one parametric and one structural constraint are defined and imposed. The method is applied to both simulated and real data and has been shown to perform very well even in low signal to noise ratio situations. In addition, the method is compared with spatial principal component analysis (sPCA) and its superiority is demonstrated by using simulated signals.

[1]  Saeid Sanei,et al.  Simultaneous localization and separation of biomedical signals by tensor factorization , 2009, 2009 IEEE/SP 15th Workshop on Statistical Signal Processing.

[2]  N. Higham Computing real square roots of a real matrix , 1987 .

[3]  Florian Roemer,et al.  Multi-dimensional PARAFAC2 component analysis of multi-channel EEG data including temporal tracking , 2010, 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology.

[4]  Rasmus Bro,et al.  MULTI-WAY ANALYSIS IN THE FOOD INDUSTRY Models, Algorithms & Applications , 1998 .

[5]  J. Polich Clinical application of the P300 event-related brain potential. , 2004, Physical medicine and rehabilitation clinics of North America.

[6]  Joseph Dien,et al.  The ERP PCA Toolkit: An open source program for advanced statistical analysis of event-related potential data , 2010, Journal of Neuroscience Methods.

[7]  R. Bro,et al.  A fast non‐negativity‐constrained least squares algorithm , 1997 .

[8]  Joseph Dien,et al.  Evaluation of PCA and ICA of simulated ERPs: Promax vs. infomax rotations , 2007, Human brain mapping.

[9]  Bülent Yener,et al.  Unsupervised Multiway Data Analysis: A Literature Survey , 2009, IEEE Transactions on Knowledge and Data Engineering.

[10]  José Carlos Príncipe,et al.  A New Spatiotemporal Filtering Method for Single-Trial Estimation of Correlated ERP Subcomponents , 2011, IEEE Transactions on Biomedical Engineering.

[11]  R. Bro,et al.  PARAFAC2—Part I. A direct fitting algorithm for the PARAFAC2 model , 1999 .