About this Journal Submit a Manuscript Table of Contents
Computational Intelligence and Neuroscience
Volume 2011 (2011), Article ID 129365, 9 pages
http://dx.doi.org/10.1155/2011/129365
Research Article

EEGIFT: Group Independent Component Analysis for Event-Related EEG Data

1Department of Biological and Medical Psychology, University of Bergen, Jonas Lies Vei 91, 5011 Bergen, Norway
2Mind Research Network, 1101 Yale Boulevard, N.E, Albuquerque, NM 87131, New Mexico, USA
3Department of Electrical and Computer Engineering, University of New Mexico, Albuquerque, USA

Received 1 October 2010; Revised 4 March 2011; Accepted 9 April 2011

Academic Editor: Sylvain Baillet

Copyright © 2011 Tom Eichele et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Independent component analysis (ICA) is a powerful method for source separation and has been used for decomposition of EEG, MRI, and concurrent EEG-fMRI data. ICA is not naturally suited to draw group inferences since it is a non-trivial problem to identify and order components across individuals. One solution to this problem is to create aggregate data containing observations from all subjects, estimate a single set of components and then back-reconstruct this in the individual data. Here, we describe such a group-level temporal ICA model for event related EEG. When used for EEG time series analysis, the accuracy of component detection and back-reconstruction with a group model is dependent on the degree of intra- and interindividual time and phase-locking of event related EEG processes. We illustrate this dependency in a group analysis of hybrid data consisting of three simulated event-related sources with varying degrees of latency jitter and variable topographies. Reconstruction accuracy was tested for temporal jitter 1, 2 and 3 times the FWHM of the sources for a number of algorithms. The results indicate that group ICA is adequate for decomposition of single trials with physiological jitter, and reconstructs event related sources with high accuracy.

1. Introduction

Event related brain responses to simple cognitive tasks are composed of multiple dynamic, temporally and regionally overlapping, functionally separable sub-processes which add to existing oscillatory background activity [15]. In other words, event-related processes are spatially and temporally mixed across the brain, and the scalp EEG samples a volume-conducted, spatially degraded version of the responses, where the potential at any location and latency can be considered a mixture of multiple independent sources that stem from large-scale synchronous field potentials [6, 7]. The exploration of the trial-to-trial variability in these responses provides important clues about the dynamics and adaptability of cognitive processes [2, 6, 810]

One powerful and increasingly popular method that allows for decomposition of EEG data and assessment of single trial variability is blind source separation with independent component analysis (ICA). ICA algorithms solve a two-dimensional linear mixing problem of spatially, and/or temporally independent sources [11, 12]. ICA models spatio-temporal data as a linear combination of maps and timecourses while attempting to maximize the statistical independence between either the maps (spatial ICA, sICA) or the time courses (temporal ICA, tICA). The method has general applicability to Gaussian mixtures, regarding psychophysiological and neuro imaging data it has been successfully used with averaged ERPs [13], single trial EEG [6, 7, 14], structural and functional MRI [15], and recently also in EEG-fMRI integration [1622]. Tools for data analysis with ICA are implemented for example in the academic freeware toolboxes EEGLAB [23], ICALAB [24] and GIFT [15, 25], both running in Matlab, as well as the stand-alone package FSL-MELODIC [26].

The basic ICA model applies to single subject data, thus one inherent limitation to the use of ICA in typical multi-subject/session EEG studies is that this method is not naturally suited to generalize results from a group of subjects. This is because ICAs from separate runs or participants will generate different sets of components with different order and scaling that need to be matched across datasets to allow group inferences. This is in contradistinction to the straightforward way of making group inferences from, for example, ERP component averages from selected channels and latencies in the general linear model [27, 28]. Therefore, a method combining individual components is desired where group inferences are straightforward. There are two strategies to allow for matching of independent components across individuals: one is to combine individual ICs across subjects with clustering techniques [7, 2931]. Clustering usually involves selection of suitable algorithms and features of interest, that is, topography, timecourse, spectrum and so forth by which between-subject correspondences of components are identified. This requires additional assumptions about the data and expert user input, and differences in algorithm and feature selection, as well as user bias can then create equivocal results. Alternatively, a more parsimonious approach is to create aggregate data containing observations from all subjects, directly estimate components that are consistently expressed in the population in a single set of ICs and then back-reconstruct estimated components to the individual data. This approach has so far predominantly been used for spatial ICA of fMRI [25, 26, 32]. We have recently adopted a group ICA method for parallel and joint decomposition of concurrent EEG-fMRI recordings [21, 22, 33]. Here, we present a group-level temporal ICA model based on the rationale proposed by Calhoun et al. [25] for single-trial analysis of event related EEG timeseries that we also implement in the toolbox EEGIFT, which is available from the Mind Research Networks’ medical image analysis lab webpages at http://icatb.sourceforge.net/gift/eegift_startup.php along with documentation and tutorial datasets. EEGIFT runs in Matlab (The Mathworks, Natick, MA), and employs pre-processed data from EEGLAB [23], a popular free toolbox for EEG processing which can be downloaded from http://sccn.ucsd.edu/eeglab/.

In order to make the component estimation computationally feasible, we employ a data reduction using principal component analysis (PCA). The number of components estimated from the data can be based on minimum description length [34] principle (MDL) or on other estimates of dimensionality. Due to aggregation and data reduction with PCA preceding component estimation, group ICA of EEG time-domain data is preferentially suited to the detection of components that contribute to event-related potentials. Processes that are not time/phase-locked within and across subjects, such as background rhythms cannot be satisfyingly reconstructed, in these cases the transformation of the data into the frequency domain prior to ICA decomposition is useful [35]. For time domain data it follows that the accuracy of component detection and back-reconstruction with a group model is dependent on the degree of intra- and inter-individual time and phase-locking of event related EEG processes. Similar to findings in early studies of PCA decomposition of ERP averages [36, 37], excess latency jitter results in splitting of a single source into two (or more) independent components representing the source and its approximate time derivative [21]. Here, we illustrate this dependency in a group analysis of 20 hybrid datasets consisting of three simulated event-related sources with varying degrees of latency jitter, mixed with real EEG data from 20 subjects that participated in a passive listening experiment. It is not trivial to recover reliable estimates about ERP latency jitter from real EEG data, we assume from a literature search and unpublished observations in our lab that single trial peak latencies of larger components such as the auditory N1 and P3 vary approximately 20–40 ms2, roughly corresponding to the full width at half maximum (FWHM) of these components (for comparison, see, e.g., [8]). We tested the reconstruction accuracy (RA), expressed as the variance of the source accounted for by the reconstructed IC (in terms of 𝑅2) for latency jitter 1, 2 and 3 times the average FWHM of the sources for the Infomax, ERICA, JADE, fastICA, and SIMBEC algorithms that are implemented among others in GIFT and EEGIFT.

2. Method

2.1. Group tICA

The group ICA model is divided into the underlying data generation and mixing process, recording, pre-processing, reduction, component estimation and back-reconstruction (schematically illustrated in Figure 1). We assume that the scalp EEG signal is a gaussian mixture containing statistically independent non-gaussian source timeseries 𝑠(𝑡)=[𝑠1(𝑡),𝑠2(𝑡),,𝑠𝑁(𝑡)]𝑇 indicated by 𝑠𝑖(𝑡) at time 𝑡 for the 𝑖th source from 𝑁 sources. The sources have weights that specify the contribution to each timepoint. The weights are multiplied by each source’s fixed topography. Secondly, it is assumed that the 𝑁 sources are linearly mixed so that a given time point contains a weighted mixture of the sources. The linear combination of sources is represented by the unknown mixing system 𝐴 where 𝐴𝑠(𝑡)=𝑢(𝑡) where 𝑢(𝑡)=[𝑢1(𝑡),𝑢2(𝑡),,𝑢𝑁(𝑡)]𝑇 and represents 𝑁 ideal samples of the signal 𝑢𝑖(𝑡) at time 𝑡, for the 𝑖th source in the brain. The sampling of the electric activity on the scalp results in 𝑦(𝑡)=[𝑦1(𝑡),𝑦2(𝑡),,𝑦𝐾(𝑡)]𝑇  where the EEG is sampled at 𝐾 timepoints where 𝑡{1,2,,𝐾}. A set of possible transformations during pre-processing, such as downsampling and filtering determine the effective sampling such that 𝑦(𝑗)=[𝑦1(𝑗),𝑦2(𝑗),,𝑦𝐾(𝑗)]𝑇 where the effective temporal sampling is indexed by 𝑗=1,2,,𝐾.

129365.fig.001
Figure 1: Group ICA. In the group ICA model, we assume that the EEG is a linear mixture of temporally independent sources in each subject 𝑠(𝑡). The linear combination of sources is represented by the unknown mixing matrix 𝐴, and yields the ideal samples of brain activity u(t), and the signals recorded with the EEG amplifier (𝐸). Transformations (𝑇) during preprocessing contain filtering, epoching, artefact rejection, individual ICA for additional artefact reduction and so forth, altering the effective temporal sampling and dimensionality of the data 𝑦(𝑖). For each individual separately, the pre-processed single trial data are pre-whitened and reduced to 𝑅 via PCA. Group data is generated by concatenating individual principal components in the aggregate data set 𝐺. Temporal ICA is performed in this set, estimating aggregate components (𝐶). From the aggregate components, the individual data are reconstructed (see text for details).
2.2. Data Reduction

For each individual separately, the pre-processed single trial data 𝑦(𝑗) are pre-whitened and reduced via PCA (Figure 1,𝑅11,,𝑅𝑀1) containing the major proportion of variance in the N uncorrelated timecourses of 𝑥(𝑗)=[𝑥1(𝑗),𝑥2(𝑗),,𝑥𝑁(𝑗)]𝑇. PCA whitening preconditions the data and simplifies ICA estimation due to the orthogonal projection, reduction of complexity, and de-noising, as well as compressing the data and thus reducing the computational load. Group data is generated by concatenating individual principal components in the aggregate data set 𝐺. In detail, let 𝑋𝑖=𝑅𝑖1𝑌𝑖 be the 𝐿-𝑏𝑦-𝑉 reduced data matrix from subject 𝑖 where 𝑌𝑖 is the 𝑄-𝑏𝑦-𝑉 data matrix containing preprocessed EEG epochs from all channels, 𝑅𝑖1 is the 𝐿-𝑏𝑦-𝑄 reducing matrix from the principal component decomposition, 𝑉 is the number of timepoints (samples per epoch * trials), 𝑄 is the number of scalp channels, and 𝐿 is the size of the channel dimension following reduction. The next step is to concatenate the reduced data from all subjects into a matrix and reduce this matrix to 𝑁, the number of components to be estimated. The 𝑁-𝑏𝑦-𝑉 reduced, concatenated matrix for the 𝑀 subjects is𝑋=𝐺1𝑅11𝑌1𝑅𝑀1𝑌𝑀,(1) where 𝐺1 is an 𝑁-𝑏𝑦-𝐿𝑀 reducing matrix from a second PCA decomposition and is multiplied on the right by the 𝐿𝑀-𝑏𝑦-𝑉 concatenated data matrix for the 𝑀 subjects.

2.3. ICA Estimation

The idea is to find the mixing matrix 𝐴 and compute the 𝑠𝑖 sources for the group. After concatenation of individual principal components in the aggregate data set 𝐺, this matrix is decomposed by 𝐼𝐶𝐴, estimating the optimal inverse of the mixing matrix 𝐴, and a single set of source timecourses (𝑆). Following 𝐼𝐶𝐴 estimation, we can write 𝐴𝑆𝑋=, where 𝐴 is the 𝑁-𝑏𝑦-𝑁 mixing matrix and 𝑆 are the 𝑁-𝑏𝑦-𝑉 component timecourses. Substituting this expression for 𝑋 into (1) and multiplying both sides by 𝐺 results in𝐺𝐴𝑅𝑆=11𝑌1𝑅𝑀1𝑌𝑀.(2)

2.4. Partitioning and Single Subject Reconstruction

Partitioning the matrix 𝐺 by subject provides the following expression𝐺1𝐺𝑀𝐴𝑅𝑆=11𝑌1𝑅𝑀1𝑌𝑀.(3) We then write the equation for subject 𝑖 by working only with the elements in partition 𝑖 of the above matrices such that 𝐺𝑖𝐴𝑆𝑖=𝑅𝑖1𝑌𝑖.(4) The matrix 𝑆𝑖 in (4) contains the single subject component timecourses for subject 𝑖, calculated from the following equation𝑆𝑖=𝐺𝑖𝐴1𝑅𝑖1𝑌𝑖.(5) We now multiply both sides of (4) by 𝑅𝑖 and write𝑌𝑖𝐹𝑖𝐺𝑖𝐴𝑆𝑖(6) yielding the ICA decomposition of the data from subject 𝑖 contained in the matrix 𝑌𝑖. The 𝑁-𝑏𝑦-𝑉 matrix 𝑆𝑖contains the 𝑁 component timecourses, and the 𝑄-𝑏𝑦-𝑁 matrix 𝐹𝑖𝐺𝑖𝐴 is the single subject mixing matrix, yielding the scalp maps for 𝑁 components.

2.5. Generation of Hybrid Data

In this simulation, 20 mutually uncorrelated hybrid EEG datasets were generated containing 63 channels, 256 timepoints and 500 trials. Three event-related sources (S1–S3) with variable topographies across datasets were mixed with real EEG data from 20 participants from an unrelated study. For each single trial, an event related response (ERR) was simulated with two Gaussians (7). ERR𝑥=𝑎1𝑒((𝑥𝑏)/3𝑐1)223𝑎2𝑒((𝑥𝑏)/2𝑐2)2.(7) The amplitudes 𝑎1 and 𝑎2 were varied randomly and independently from 0.5–2.5 and the widths 𝑐1 and 𝑐2 from 0.5–1.5, introducing additional jitter of the ERR amplitude and shape. Latency jitter is a relevant source of variability in single trials, affecting the accuracy of component estimation [21, 36, 37]. The three sources simulated here had non-overlapping peak latencies, and latency variability (within-“subject”) in 𝑏 was in a range of 20 samples in S1, corresponding to the FWHM of the ERR, 40 samples (2 FWHM) in S2, and 60 samples (3 FWHM) in S3. Sine waves with random phase and amplitude modulation were additionally entered as background activity into each source. Across individual datasets the average peak latency of each source dataset was varied by 20 samples (between-“subjects”). For each source, the scalp distributions were generated as dipolar maps covering six channels (of 63), with 50% overlap between S1-S2, and S2-S3, respectively. Across datasets, the location of each source was systematically varied. These sources were normalized to unit variance, and mixed with normalized real EEG data with the same dimensions from 20 participants. The resulting hybrid data are shown in the top half of Figure 2 for two datasets.

129365.fig.002
Figure 2: Hybrid Data. The figure shows two hybrid datasets with three source topographies (𝐴) and timecourses (𝑠) on top, the mixture of the sources with real EEG (𝑋) in (𝑏), and the reconstruction of the sources after group ICA (𝑊,𝑦) in (𝑐) section.
2.6. Independent Component Analysis

In order to generate a reference value for the performance of group ICA we computed individual ICA solutions in EEGLAB for each of the datasets, employing the Infomax algorithm [38]. For the group ICA, all subjects were analyzed at once, and principal component analysis (PCA) was used for compression to allow the datasets to be processed together. The number of components is estimated by doing singular value decomposition on the data and the resulting eigenvalues are passed to MDL method [34]. Here we selected 20 components as the top 20 components usually explain more than 95% of the variance in the data. In our experience, the exact choice of the number of components does not critically affect the results as long as this number is not much smaller than the true number of sources. In the PCA steps, data from each dataset was reduced over the spatial dimension, that is. from the number of channels to 20 principal components, concatenated across subjects, and again reduced to 20 components. Temporal ICA was then performed using the Infomax algorithm [38] with subsequent back-reconstruction into single subjects.In order to assess reconstruction accuracy of group ICA for different numbers of estimated components, we estimated the solutions for 10, 20, 30, 40 and 50 components using Infomax. For comparison between algorithms, we also estimated solutions using the fastICA [39], JADE [40], SIMBEC [41] and ERICA [42] algorithms. The reconstruction accuracy of group ICA was expressed as the variance of the simulated sources accounted for by the reconstructed ICs averaged across the 20 datasets (𝑅2), separately for the entire single trial images, the amplitude modulation across trials around the component peak latency (averaged in a 20 sample window), the component average timeseries, and topographies (Table 1). Results for two hybrid datasets with variable topographic and temporal representations of the three sources are illustrated in Figure 2.

tab1
Table 1: Reconstruction Accuracy. The table summarizes the reconstruction accuracy (RA, mean across datasets ± S.E.M.) of different group ICA models. RA stands for proportion of the source variability accounted for by the ICs averaged across the 20 datasets, and was computed separately for the entire single trial images, the amplitude modulation across trials around the component peak latency, the component average timeseries, and topographies.
2.7. Application to Real Data

In addition to the quantification of the model performance, we illustrate group ICA in the context of typical recording conditions and pre-processing steps that we employ with a decomposition of an auditory oddball dataset. 32 healthy participants took part in the experiment after providing a written statement of informed consent. Participants were sitting in an electro-magnetically shielded and sound-attenuated testing chamber (Rainford EMC Systems, Wigan, UK) and were fitted with 61 Ag/AgCl scalp electrodes mounted in an elastic cap (EasiCap, Falk Minow Services, Breitenbrunn, Germany) and two additional channels monitoring eye movements. All channels were referenced to the nose, and impedances were kept below 10 kΩ. EEGs were recorded continuously at 500 Hz sampling frequency with a band-pass from  .01–250 Hz with BrainAmp DC amplifiers (BrainProducts, Munich, Germany). The experiment consisted of detecting an infrequent target sound within a series of frequent regular sounds and participants were asked to respond as quickly as possible by pressing a button with their right index finger. The standard stimulus was a center-panned 500 Hz tone, the target stimuli were left or right-panned location deviants, respectively. Targets occurred with a probability of 0.1 for each location. Stimulus duration was 75 ms and the inter-stimulus interval was 700 ms, with no immediate repetitions of targets. All stimuli were presented at approximately 65 decibels above threshold. EEGs were down-sampled to 250 Hz, filtered with a zero-phase Butterworth filter from 1–45 Hz, and re-referenced to common average reference. In the example application, the data were segmented from −800 to 1200 ms around target stimuli and thus contain a sequence of standard-target-standard sounds, to dissociate between obligatory stimulus-related and target-related components, respectively. Trials with amplitudes exceeding ±150μV on any of the channels were excluded from further analysis. Concatenated single sweeps around target onset were subjected to single subject Infomax ICA in EEGLAB [23] running in MATLAB. Components with topographies and timecourses attributeable artifacts were identified and removed from the data [43]. For each dataset, we extract 20 components from the data based on the MDL method. Missing trials were padded with the mean from surrounding trials because there are gradual changes across trials. Single-trials were additionally denoised with a wavelet filter [44]. Hereafter, group ICA was computed, estimating 20 components.

3. Results

Reconstruction accuracy (RA) for all analyses is summarized in Table 1. Individual ICA reconstructs the source timecourses with near-perfect accuracy and independent of latency jitter, the topographies are recovered with an accuracy of around 0.9. Group ICA models yield overall lower RA than individual ICA, coming closest to individual ICA in S1, where the entire source timecourse is reconstructed with 0.86 on the average, Infomax with 20 components yields the best performance at 0.88. RA for the amplitude modulation around the peak latency of S1 is 0.9 on the average, for the average timecourses RA is 0.99. These two latter features are typically most relevant for making inferences about electrophysiological data, that is. the overall shape of the waveform and the (single trial) peak amplitudes of components. The topographies of S1 are reconstructed with an accuracy of 0.87. The group ICA result for all features falls with increasing latency jitter regardless of algorithm choice. For the entire timecourse, average RA across algorithms is S1—0.86, S2—0.76, S3—0.50, with the most pronounced drop-off of all features. The peak RA for the peak is S1—0.90, S2—0.87, S3—0.65. RA for the timecourse average is S1—0.95, S2—0.95, S3—0.82, the topographies are reconstructed with S1—0.87, S2—0.80, S3—0.77 overall a less pronounced drop-off. This effect is more pronounced with increasing the number of estimated components as can be seen in the RA of all features of S3. In summary, all three sources were recovered with sufficient reconstruction accuracy of all four features.

The decomposition of the real data yielded a number of event related components that showed differential responses between standard and target sounds. We related these in terms of the topography and peak latency to the ERP components N1, T-complex, P2, N2 and subcomponents of the P3. In Figure 3, we show the group mean of one independent component that represents the auditory N1, together with the reconstructions for three representative subjects. The group mean of the N1 component is representative for the group level average, and as is expected somewhat smaller in amplitude than compared to a single subject’s mean amplitude because of inter-individual latency differences. Parts of this dataset accompany the toolbox as tutorial material and the entire dataset is available upon request.

129365.fig.003
Figure 3: Real Data Decomposition. This figure is an illustration of one component from a group ICA decomposition of an auditory oddball dataset (n=32 participants), and shows the group mean of one independent component in the leftmost column, and reconstructions for three subject s in the other columns. Note that the group mean component peaks have smaller amplitude due to latency jitter. Top: topography, Middle: single trial image, bottom: component event related average. Vertical lines in the lower left indicate the onset of the preceding standard, target, and succeeding standard, respectively. Topography and latency identify this component as the auditory N1, and typical N1-enhancement is clearly visible in response to targets.

4. Discussion

This work presents an approach to perform a temporal independent component analysis on single-trial time domain EEG data for multiple subjects simultaneously. Our model uses a combination of principal component analysis for data reduction, subsequent independent component analysis on the aggregate data, and back-reconstruction of the aggregate mixing matrix in individual Subjects [20, 25]. This method is implemented in the freeware toolbox EEGIFT that runs in Matlab (R13 and newer) and is downloadable from http://mialab.mrn.org/ or http://icatb.sourceforge.net/. EEGIFT has a graphical user interface (GUI) that allow import of EEG data from multiple participants pre-processed in EEGLAB (http://sccn.ucsd.edu/eeglab/). Another GUI allows the user to specify of analysis details, such as PCA data reduction, model order, choice of ICA algorithm, and the respective parameters. EEGIFT also allow robust estimation with ICASSO [45]. The analysis output is stored after back-reconstruction as individual timecourses and topographies in Matlab format which can be used for specifying between-condition and/or between-groups statistical tests. Individual data, group averages, and population statistics can also be visualized in a GUI as topographies, grand mean timecourses and single trial images. Analyses can also be batch-scripted for convenience. A full documentation of functions in EEGIFT, including a tutorial walkthrough, accompanies the download package (http://icatb.sourceforge.net/gift/eegift_startup.php). User support is provided through the GIFT mailing list (Icatb-discuss).

In EEGIFT, the back-reconstructed timecourses and topographies are a function primarily of the variability within subjects, as opposed to representing a representation of the average across subjects. A simulation affords to create a situation of a known ground truth with sources and noise parameters, and is useful to illustrate of the concept of the group model. Here, we generated hybrid data with realistic spatial, temporal and amplitude variability. The reconstruction of the solutions with different algorithms and numbers of components demonstrates that this approach offers a straightforward and computationally tractable solution to the problem of multi-subject analysis with ICA. The results from the decomposition of hybrid dataset illustrate that this group model is able to recover timecourses and topographies of event related responses on single trial level with overall sufficiently high accuracy. As laid out in the introduction, the critical determinant for the success of a group model computed for time-domain data is intra- and inter-subject latency jitter. There is considerable difficulty in reliably estimating latency jitter of event related responses in real EEG data; we assume here that physiologic jitter is roughly on the order of the full width at half maximum of sources, in which case the reconstruction accuracy of the group model reaches more than 90% of the accuracy of the individual ICA estimates. Individual ICA with subsequent ordering of components across subject would achieve the highest accuracy with an ideal clustering technique. We did not attempt to directly compare the performance of group-level ICA with that of individual ICA with subsequent clustering, apart from possible computational limitations, the challenge for existing algorithms is to identify and cluster components with poor between-subject correspondence of their topographies and timecourses [7, 2931]. Group ICA avoids any such problems since the decomposition is computed for all data-sets/subjects simultaneously, estimating a single set of components with the same order across datasets. Consequently, it is straightforward to perform random-effects population tests for the timecourses as well as the topographies. This can also be done in the statistical parametric mapping framework where testing within the Gaussian random field theory with adjustments for multiple comparisons (e.g., [28, 46, 47]). Apart from being used as a primary tool for inference, we believe that there is also a utility for group ICA in the mining stage of prediction-based analysis: For example, in cases where a-priori models of the event related responses can only be poorly specified, for example, where the selection criteria for appropriate electrodes and time-windows from the EEG data are unclear, group ICA results can be used for model specification.

As noted in Section 1, the limitation of the current method is that responses with poor time/phase-locking are not satisfyingly reconstructed [21, 36, 37]. This is a consequence of the data reduction and aggregation, which inherently limits the visibility of this method to evoked activity, both within and across subjects. Sources that have a loose relation to stimulus/response onset, if extracted and identified, are usually represented in more than one component. This is in clear contradistinction to the superior performance of temporal ICA on concatenated EEG epochs from single subjects, which is insensitive to trial-to-trial phase/latency variability of sources. For detection of poorly time-locked processes in this framework one could consider time-domain data with frequency or time-frequency transformed data [35], or incorporate correction for latency jitter during pre-processing where feasible [8, 44]. We assume from testing the performance of this model extensively in simulated, hybrid, and real data that it is a very useful addition to the available ICA-toolbelt, in that it affords a straightforward possibility to perform group analysis of event related EEG responses. Since the underlying generative model is flexible and modality independent, and the software implementations are highly interrelated, a genuine advantage of this method is that the EEG components can be fused with results from diverse biomedical imaging modalities such as sMRI, fMRI, DTI and VBMas well as genetic information (SNP) within the same conceptual and computational framework [15, 18, 20, 23, 48, 49]. This then affords multimodal inferences which can bear novel insights about brain function.

Acknowledgment

The present study was financially supported with a BILATGRUNN Grant from the Research Council of Norway to Tom Eichele and by the National Institutes of Health, under Grants 1 R01 EB 000840 and 1 R01 EB 005846 to Vince Calhoun.

References

  1. E. Halgren and K. Marinkovic, “General principles for the physiology of cognition as suggested by intracranial ERPs,” in Recent Advances in Event-Related Brain Potential Research, C. Ogura, Y. Koga, and M. Shimokochi, Eds., pp. 1072–1084, Elsevier, Amsterdam, The Netherlands, 1995.
  2. A. Arieli, A. Sterkin, A. Grinvald, and AD. Aertsen, “Dynamics of ongoing activity: explanation of the large variability in evoked cortical responses,” Science, vol. 273, no. 5283, pp. 1868–1871, 1996. View at Scopus
  3. E. Halgren, K. Marinkovic, and P. Chauvel, “Generators of the late cognitive potentials in auditory and visual oddball tasks,” Electroencephalography and Clinical Neurophysiology, vol. 106, no. 2, pp. 156–164, 1998. View at Publisher · View at Google Scholar · View at Scopus
  4. T. Eichele, K. Specht, M. Moosmann et al., “Assessing the spatiotemporal evolution of neuronal activation with single-trial event-related potentials and functional MRI,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 49, pp. 17798–17803, 2005. View at Publisher · View at Google Scholar · View at Scopus
  5. K. A. Kiehl, M. C. Stevens, K. R. Laurens, G. Pearlson, V. D. Calhoun, and P. F. Liddle, “An adaptive reflexive processing model of neurocognitive function: supporting evidence from a large scale (n = 100) fMRI study of an auditory oddball task,” NeuroImage, vol. 25, no. 3, pp. 899–915, 2005. View at Publisher · View at Google Scholar · View at Scopus
  6. S. Makeig, S. Debener, J. Onton, and A. Delorme, “Mining event-related brain dynamics,” Trends in Cognitive Sciences, vol. 8, no. 5, pp. 204–210, 2004. View at Publisher · View at Google Scholar · View at Scopus
  7. J. Onton and S. Makeig, “Information-based modeling of event-related brain dynamics,” in Progress in Brain Research, C. Neuper and G. Pfurtscheller, Eds., vol. 159, pp. 99–120, Elsevier, Amsterdam, The Netherlands, 2006.
  8. K. M. Spencer, “Averaging, detection, and classification of single-trial ERPs,” in Event Related Potentials, T. C. Handy, Ed., pp. 209–228, The MIT, Cambridge, Mass, USA, 2005.
  9. S. Debener, M. Ullsperger, M. Siegel, and A. K. Engel, “Single-trial EEG-fMRI reveals the dynamics of cognitive function,” Trends in Cognitive Sciences, vol. 10, no. 12, pp. 558–563, 2006. View at Publisher · View at Google Scholar · View at Scopus
  10. H. Eichele and H. T. Juvodden, “Mal-adaptation of event-related EEG responses preceding performance errors,” Frontiers in Human Neuroscience, vol. 10, no. 4, p. 65, 2010.
  11. A. Hyvärinen and E. Oja, “Independent component analysis: algorithms and applications,” Neural Networks, vol. 13, no. 4-5, pp. 411–430, 2000. View at Publisher · View at Google Scholar · View at Scopus
  12. J. V. Stone, “Independent component analysis: an introduction,” Trends in Cognitive Sciences, vol. 6, no. 2, pp. 59–64, 2002. View at Publisher · View at Google Scholar · View at Scopus
  13. S. Makeig, T. P. Jung, A. J. Bell, D. Ghahremani, and T. J. Sejnowski, “Blind separation of auditory event-related brain responses into independent components,” Proceedings of the National Academy of Sciences of the United States of America, vol. 94, no. 20, pp. 10979–10984, 1997. View at Publisher · View at Google Scholar · View at Scopus
  14. T. P. Jung, S. Makeig, M. Westerfield, J. Townsend, E. Courchesne, and T. J. Sejnowski, “Analysis and visualization of single-trial event-related potentials,” Human Brain Mapping, vol. 14, no. 3, pp. 166–185, 2001. View at Publisher · View at Google Scholar · View at Scopus
  15. V. D. Calhoun and T. Adali, “Unmixing fMRI with independent component analysis,” IEEE Engineering in Medicine and Biology Magazine, vol. 25, no. 2, pp. 79–90, 2006. View at Publisher · View at Google Scholar · View at Scopus
  16. S. Debener, M. Ullsperger, M. Siegel, K. Fiehler, D. Y. Von Cramon, and A. K. Engel, “Trial-by-trial coupling of concurrent electroencephalogram and functional magnetic resonance imaging identifies the dynamics of performance monitoring,” Journal of Neuroscience, vol. 25, no. 50, pp. 11730–11737, 2005. View at Publisher · View at Google Scholar · View at Scopus
  17. B. Feige, K. Scheffler, F. Esposito, F. Di Salle, J. Hennig, and E. Seifritz, “Cortical and subcortical correlates of electroencephalographic alpha rhythm modulation,” Journal of Neurophysiology, vol. 93, no. 5, pp. 2864–2872, 2005. View at Publisher · View at Google Scholar · View at Scopus
  18. V. D. Calhoun, T. Adali, N. R. Giuliani, J. J. Pekar, K. A. Kiehl, and G. D. Pearlson, “Method for multimodal analysis of independent source differences in schizophrenia: combining gray matter structural and auditory oddball functional data,” Human Brain Mapping, vol. 27, no. 1, pp. 47–62, 2006. View at Publisher · View at Google Scholar · View at Scopus
  19. V. D. Calhoun, T. Adali, G. D. Pearlson, and K. A. Kiehl, “Neuronal chronometry of target detection: fusion of hemodynamic and event-related potential data,” NeuroImage, vol. 30, no. 2, pp. 544–553, 2006. View at Publisher · View at Google Scholar · View at Scopus
  20. T. Eichele, V. D. Calhoun, M. Moosmann et al., “Unmixing concurrent EEG-fMRI with parallel independent component analysis,” International Journal of Psychophysiology, vol. 67, no. 3, pp. 222–234, 2008. View at Publisher · View at Google Scholar · View at Scopus
  21. M. Moosmann, T. Eichele, H. Nordby, K. Hugdahl, and V. D. Calhoun, “Joint independent component analysis for simultaneous EEG-fMRI: principle and simulation,” International Journal of Psychophysiology, vol. 67, no. 3, pp. 212–221, 2008. View at Publisher · View at Google Scholar · View at Scopus
  22. T. Eichele, V. D. Calhoun, and S. Debener, “Mining EEG-fMRI using independent component analysis,” International Journal of Psychophysiology, vol. 73, no. 1, pp. 53–61, 2009. View at Publisher · View at Google Scholar · View at Scopus
  23. A. Delorme and S. Makeig, “EEGLAB: an open source toolbox for analysis of single-trial EEG dynamics including independent component analysis,” Journal of Neuroscience Methods, vol. 134, no. 1, pp. 9–21, 2004. View at Publisher · View at Google Scholar · View at Scopus
  24. A. Cichocki, S.-I. Amari, et al., Adaptive Blind Signal and Image Processing, Wiley, New York, NY, USA, 2002.
  25. V. D. Calhoun, T. Adali, G. D. Pearlson, and J. J. Pekar, “A method for making group inferences from functional MRI data using independent component analysis,” Human Brain Mapping, vol. 14, no. 3, pp. 140–151, 2001. View at Publisher · View at Google Scholar · View at Scopus
  26. C. F. Beckmann and S. M. Smith, “Tensorial extensions of independent component analysis for multisubject FMRI analysis,” NeuroImage, vol. 25, no. 1, pp. 294–311, 2005. View at Publisher · View at Google Scholar · View at Scopus
  27. T. W. Picton, S. Bentin, P. Berg et al., “Guidelines for using human event-related potentials to study cognition: recording standards and publication criteria,” Psychophysiology, vol. 37, no. 2, pp. 127–152, 2000. View at Publisher · View at Google Scholar · View at Scopus
  28. S. J. Kiebel and K. J. Friston, “Statistical parametric mapping for event-related potentials: I. Generic considerations,” NeuroImage, vol. 22, no. 2, pp. 492–502, 2004. View at Publisher · View at Google Scholar · View at Scopus
  29. F. Esposito, T. Scarabino, A. Hyvarinen et al., “Independent component analysis of fMRI group studies by self-organizing clustering,” NeuroImage, vol. 25, no. 1, pp. 193–205, 2005. View at Publisher · View at Google Scholar · View at Scopus
  30. J. Onton, M. Westerfield, J. Townsend, and S. Makeig, “Imaging human EEG dynamics using independent component analysis,” Neuroscience and Biobehavioral Reviews, vol. 30, no. 6, pp. 808–822, 2006. View at Publisher · View at Google Scholar · View at Scopus
  31. F. de Martino, F. Gentile, F. Esposito et al., “Classification of fMRI independent components using IC-fingerprints and support vector machine classifiers,” NeuroImage, vol. 34, no. 1, pp. 177–194, 2007. View at Publisher · View at Google Scholar · View at Scopus
  32. V. J. Schmithorst and S. K. Holland, “Comparison of three methods for generating group statistical inferences from independent component analysis of functional magnetic resonance imaging data,” Journal of Magnetic Resonance Imaging, vol. 19, no. 3, pp. 365–368, 2004. View at Publisher · View at Google Scholar · View at Scopus
  33. T. Eichele and V. Calhoun, “Parallel EEG-fMRI ICA decomposition,” in Integrating EEG and FMRI, Analysis and Application, M. Ullsperger and S. Debener, Eds., Oxford University Press, New York, NY, USA, 2010.
  34. J. Rissanen, “A universal prior for integers and estimation by minimum description length,” Annals of Statistics, vol. 11, no. 2, pp. 416–431, 1983. View at Publisher · View at Google Scholar
  35. L. Wu, T. Eichele, and V. D. Calhoun, “Reactivity of hemodynamic responses and functional connectivity to different states of alpha synchrony: a concurrent EEG-fMRI study,” NeuroImage, vol. 52, no. 4, pp. 1252–60, 2010. View at Publisher · View at Google Scholar · View at Scopus
  36. E. Donchin and E. F. Heffley, “Multivariate analysis of event-related potential data: a tutorial review,” in Multidisciplinary Perspectives in Event-Related Potential Research, D. Otto, Ed., Government Printing Office, Washington, DC, USA, 1978.
  37. J. Mocks, “The influence of latency jitter in principal component analysis of event-related potentials,” Psychophysiology, vol. 23, no. 4, pp. 480–484, 1986. View at Scopus
  38. A. J. Bell and T. J. Sejnowski, “An information-maximization approach to blind separation and blind deconvolution,” Neural Computation, vol. 7, no. 6, pp. 1129–1159, 1995. View at Scopus
  39. A. Hyvärinen and E. Oja, “A fast fixed-point algorithm for independent component analysis,” Neural Computation, vol. 9, no. 7, pp. 1483–1492, 1997. View at Scopus
  40. J. F. Cardoso and A. Souloumiac, “Blind beamforming for non Gaussian signals,” IEE-Proceedings-F, vol. 140, no. 6, pp. 362–370, 1993. View at Scopus
  41. S. Cruces, A. Cichocki, et al., “Criteria for the simultaneous extraction of arbitrary groups of sources,” in Proceedings of the International Conference on Independent Component Analysis and Blind Signal Separation, San Diego, Calif, USA, December 2001.
  42. S. Cruces, L. Castedo, A. Cichocki, et al., “Novel blind source separation algorithms using cumulants,” in Proceedings of the IEEE Interntional Conference on Acoustics, Speech, and Signal Processing, pp. 3152–3155, Istanbul, Turkey, June 2000. View at Scopus
  43. T. P. Jung, S. Makeig, C. Humphries et al., “Removing electroencephalographic artifacts by blind source separation,” Psychophysiology, vol. 37, no. 2, pp. 163–178, 2000. View at Publisher · View at Google Scholar · View at Scopus
  44. R. Quian Quiroga and H. Garcia, “Single-trial event-related potentials with wavelet denoising,” Clinical Neurophysiology, vol. 114, no. 2, pp. 376–390, 2003. View at Publisher · View at Google Scholar · View at Scopus
  45. J. Himberg, A. Hyvärinen, and F. Esposito, “Validating the independent components of neuroimaging time series via clustering and visualization,” NeuroImage, vol. 22, no. 3, pp. 1214–1222, 2004. View at Publisher · View at Google Scholar · View at Scopus
  46. K. J. Friston, A. P. Holmes, K. J. Worsley, J. P. Poline, C. D. Frith, and R. S. J. Frackowiak, “Statistical parametric maps in functional imaging: a general linear approach,” Human Brain Mapping, vol. 2, no. 4, pp. 189–210, 1995. View at Scopus
  47. J. M. Kilner, S. J. Kiebel, and K. J. Friston, “Applications of random field theory to electrophysiology,” Neuroscience Letters, vol. 374, no. 3, pp. 174–178, 2005. View at Publisher · View at Google Scholar · View at Scopus
  48. J. Liu, G. Pearlson, et al., “Combining fMRI and SNP data to investigate connections between brain function and genetics using parallel ICA,” Human Brain Mapping, vol. 30, no. 1, pp. 241–255, 2007.
  49. V. D. Calhoun, J. Liu, and T. Adali, “A review of group ICA for fMRI data and ICA for joint inference of imaging, genetic, and ERP data,” NeuroImage, vol. 45, no. 1, pp. S163–S172, 2009. View at Publisher · View at Google Scholar · View at Scopus