Current analytical techniques applied to functional magnetic resonance imaging (fMRI) data require a priori knowledge or specific assumptions about the time courses of processes contributing to the measured signals. Here we describe a new method for analyzing fMRI data based on the independent component analysis (ICA) algorithm of Bell and Sejnowski ([1995]: Neural Comput 7:1129-1159). We decomposed eight fMRI data sets from 4 normal subjects performing Stroop color-naming, the Brown and Peterson work/number task, and control tasks into spatially independent components. Each component consisted of voxel values at fixed three-dimensional locations (a component "map"), and a unique associated time course of activation. Given data from 144 time points collected during a 6-min trial, ICA extracted an equal number of spatially independent components. In all eight trials, ICA derived one and only one component with a time course closely matching the time course of 40-sec alternations between experimental and control tasks. The regions of maximum activity in these consistently task-related components generally overlapped active regions detected by standard correlational analysis, but included frontal regions not detected by correlation. Time courses of other ICA components were transiently task-related, quasiperiodic, or slowly varying. By utilizing higher-order statistics to enforce successively stricter criteria for spatial independence between component maps, both the ICA algorithm and a related fourth-order decomposition technique (Comon [1994]: Signal Processing 36:11-20) were superior to principal component analysis (PCA) in determining the spatial and temporal extent of task-related activation. For each subject, the time courses and active regions of the task-related ICA components were consistent across trials and were robust to the addition of simulated noise. Simulated movement artifact and simulated task-related activations added to actual fMRI data were clearly separated by the algorithm. ICA can be used to distinguish between nontask-related signal components, movements, and other artifacts, as well as consistently or transiently task-related fMRI activations, based on only weak assumptions about their spatial distributions and without a priori assumptions about their time courses. ICA appears to be a highly promising method for the analysis of fMRI data from normal and clinical populations, especially for uncovering unpredictable transient patterns of brain activity associated with performance of psychomotor tasks.

Objectives: Electrical potentials produced by blinks and eye movements present serious problems for electroencephalographic (EEG) and event-related potential (ERP) data interpretation and analysis, particularly for analysis of data from some clinical populations. Often, all epochs contaminated by large eye artifacts are rejected as unusable, though this may prove unacceptable when blinks and eye movements occur frequently.Methods: Frontal channels are often used as reference signals to regress out eye artifacts, but inevitably portions of relevant EEG signals also appearing in EOG channels are thereby eliminated or mixed into other scalp channels. A generally applicable adaptive method for removing artifacts from EEG records based on blind source separation by independent component analysis (ICA) (Neural Computation 7 (1995) 1129; Neural Computation 10(8) (1998) 2103; Neural Computation 11(2) (1999) 606) overcomes these limitations.Results: Results on EEG data collected from 28 normal controls and 22 clinical subjects performing a visual selective attention task show that ICA can be used to effectively detect, separate and remove ocular artifacts from even strongly contaminated EEG recordings. The results compare favorably to those obtained using rejection or regression methods.Conclusions: The ICA method can preserve ERP contributions from all of the recorded trials and all the recorded data channels, even when none of the single trials are artifact-free. q

Averaged event-related potential (ERP) data recorded from the human scalp reveal electroencephalographic (EEG) activity that is reliably time-locked and phaselocked to experimental events. We report here the application of a method based on information theory that decomposes one or more ERPs recorded at multiple scalp sensors into a sum of components with fixed scalp distributions and sparsely activated, maximally independent time courses. Independent component analysis (ICA) decomposes ERP data into a number of components equal to the number of sensors. The derived components have distinct but not necessarily orthogonal scalp projections. Unlike dipole-fitting methods, the algorithm does not model the locations of their generators in the head. Unlike methods that remove second-order correlations, such as principal component analysis (PCA), ICA also minimizes higherorder dependencies. Applied to detected-and undetectedtarget ERPs from an auditory vigilance experiment, the algorithm derived ten components that decomposed each of the major response peaks into one or more ICA components with relatively simple scalp distributions. Three of these components were active only when the subject detected the targets, three other components only when the target went undetected, and one in both cases. Three additional components accounted for the steady-state brain response to a 39-Hz background click train. Major features of the decomposition proved robust across sessions and changes in sensor number and placement. This method of ERP analysis can be used to compare responses from multiple stimuli, task conditions, and subject states.Although the locations of the brain areas generating eventrelated potentials (ERPs) cannot be uniquely determined by scalp recordings from any number of channels (1), several methods have been proposed for decomposing evoked responses into activations of distinct neural sources. Most of these also attempt to locate the active areas, by assuming either that they have a known or simple spatial configuration (2) or that generators are restricted to a small subset of possible locations and orientations (3). Other methods based on rotations of principal components use optimization criteria not directly related to brain anatomy and physiology. These methods may assume that each response component has the same time course of activation in every experimental condition (4). All these methods use second-order spatiotemporal correlations to perform the decomposition.Here we report a statistical method for decomposing one or more event-related brain responses into a sum of components with spatially fixed scalp distributions and maximally independent (though possibly overlapping) time courses. Independence requires the absence of higher-order as well as secondorder correlations between the time courses. Independence, therefore, is a stronger condition than decorrelation and, in particular, is not satisfied by decomposition into principal components by principal component analysis (PCA).Although the ...

Abstract-Recent demand and interest in wireless, mobile-based healthcare has driven significant interest towards developing alternative biopotential electrodes for patient physiological monitoring. The conventional wet adhesive Ag/AgCl electrodes used almost universally in clinical applications today provide an excellent signal but are cumbersome and irritating for mobile use. While electrodes that operate without gels, adhesives and even skin contact have been known for many decades, they have yet to achieve any acceptance for medical use. In addition, detailed knowledge and comparisons between different electrodes are not well known in the literature. In this paper, we explore the use of dry/noncontact electrodes for clinical use by first explaining the electrical models for dry, insulated and noncontact electrodes and show the performance limits, along with measured data. The theory and data show that the common practice of minimizing electrode resistance may not always be necessary and actually lead to increased noise depending on coupling capacitance. Theoretical analysis is followed by an extensive review of the latest dry electrode developments in the literature. The paper concludes with highlighting some of the novel systems that dry electrode technology has enabled for cardiac and neural monitoring followed by a discussion of the current challenges and a roadmap going forward.

Eye movements, eye blinks, cardiac signals, muscle noise, and line noise present serious problems for electroencephalographic (EEG) interpretation and analysis when rejecting contaminated EEG segments results in an unacceptable data loss. Many methods have been proposed to remove artifacts from EEG recordings, especially those arising from eye movements and blinks. Often regression in the time or frequency domain is performed on parallel EEG and electrooculographic (EOG) recordings to derive parameters characterizing the appearance and spread of EOG artifacts in the EEG channels. Because EEG and ocular activity mix bidirectionally, regressing out eye artifacts inevitably involves subtracting relevant EEG signals from each record as well. Regression methods become even more problematic when a good regressing channel is not available for each artifact source, as in the case of muscle artifacts. Use of principal component analysis (PCA) has been proposed to remove eye artifacts from multichannel EEG. However, PCA cannot completely separate eye artifacts from brain signals, especially when they have comparable amplitudes. Here, we propose a new and generally applicable method for removing a wide variety of artifacts from EEG records based on blind source separation by independent component analysis (ICA). Our results on EEG data collected from normal and autistic subjects show that ICA can effectively detect, separate, and remove contamination from a wide variety of artifactual sources in EEG records with results comparing favorably with those obtained using regression and PCA methods. ICA can also be used to analyze blink‐related brain activity.

The analysis of electroencephalographic (EEG) and magnetoencephalographic (MEG)

Human event-related potentials (ERPs) were recorded from 10 subjects presented with visual target and nontarget stimuli at five screen locations and responding to targets presented at one of the locations. The late positive response complexes of 25-75 ERP average waveforms from the two task conditions were simultaneously analyzed with Independent Component Analysis, a new computational method for blindly separating linearly mixed signals. Three spatially fixed, temporally independent, behaviorally relevant, and physiologically plausible components were identified without reference to peaks in single-channel waveforms. A novel frontoparietal component (P3f) began at approximately 140 msec and peaked, in faster responders, at the onset of the motor command. The scalp distribution of P3f appeared consistent with brain regions activated during spatial orienting in functional imaging experiments. A longer-latency large component (P3b), positive over parietal cortex, was followed by a postmotor potential (Pmp) component that peaked 200 msec after the button press and reversed polarity near the central sulcus. A fourth component associated with a left frontocentral nontarget positivity (Pnt) was evoked primarily by target-like distractors presented in the attended location. When no distractors were presented, responses of five faster-responding subjects contained largest P3f and smallest Pmp components; when distractors were included, a Pmp component appeared only in responses of the five slower-responding subjects. Direct relationships between component amplitudes, latencies, and behavioral responses, plus similarities between component scalp distributions and regional activations reported in functional brain imaging experiments suggest that P3f, Pmp, and Pnt measure the time course and strength of functionally distinct brain processes.

A method is given for determining the time course and spatial extent of consistently and transiently taskrelated activations from other physiological and artifactual components that contribute to functional MRI (fMRI) recordings. Independent component analysis (ICA) was used to analyze two fMRI data sets from a subject performing 6-min trials composed of alternating 40-sec Stroop color-naming and control task blocks. Each component consisted of a fixed threedimensional spatial distribution of brain voxel values (a ''map'') and an associated time course of activation. For each trial, the algorithm detected, without a priori knowledge of their spatial or temporal structure, one consistently task-related component activated during each Stroop task block, plus several transiently task-related components activated at the onset of one or two of the Stroop task blocks only. Activation patterns occurring during only part of the fMRI trial are not observed with other techniques, because their time courses cannot easily be known in advance. Other ICA components were related to physiological pulsations, head movements, or machine noise. By using higherorder statistics to specify stricter criteria for spatial independence between component maps, ICA produced improved estimates of the temporal and spatial extent of task-related activation in our data compared with principal component analysis (PCA). ICA appears to be a promising tool for exploratory analysis of fMRI data, particularly when the time courses of activation are not known in advance.Univariate methods for the analysis of functional MRI (fMRI) data typically examine each brain volume element or voxel individually, to determine whether the activity level at that voxel reaches a prespecified criterion for task-related activity. A common criterion is a predetermined level of significance for a statistic, such as the Student t (1) or KolmogorovSmirnov (2) statistic, under the null hypothesis that the distribution of a voxel's values during the behavioral control task is identical to that during performance of the experimental task(s). Correlational analysis (3) determines whether the similarity between a voxel's time course and a prediction of the task-related modulation, the reference function, exceeds a specified threshold. These methods then assemble individually selected (or ''active'') voxels, ignoring statistical relationships between voxels, to create a spatially distributed map demonstrating areas of significant activation.To enhance the statistical power of standard analysis techniques based on correlation or univariate statistical tests, fMRI experimenters often use alternating task-block designs in which the subject performs two or more tasks successively in alternating 20-to 40-sec blocks. By averaging over a number of task-block cycles, small consistently task-related (CTR) differences in hemodynamic activation can be detected. Isolated stimulus paradigms, such as that employed by Buckner et al. (4), avoid overlapping hemodynamic responses produc...