Study of local EEG specificities in children with mental development disorders using independent component analysis

Baseline EEGs in the frequency range of 3–13 Hz in children with mental disorders of perinatal origin during wakefulness with the eyes open were analyzed using independent component analysis. In cases of severe mental retardation, a significant increase in the power density spectra of the θ band was revealed in the left-sided frontotemporal and right-sided temporal cortices, which allows us to consider these regions to be putative sources of slow activity and markers for a lesion or immaturity in the fronto-thalamic system, as well as for the temporal areas responsible for the auditory analysis and synthesis of speech signals and the integration of audio-visual information.     

Mining HIV dynamics using independent component analysis

MOTIVATION: We implement a data mining technique based on the method of Independent Component Analysis (ICA) to generate reliable independent data sets for different HIV therapies. We show that this technique takes advantage of the ICA power to eliminate the noise generated by artificial interaction of HIV system dynamics. Moreover, the incorporation of the actual laboratory data sets into the analysis phase offers a powerful advantage when compared with other mathematical procedures that consider the general behavior of HIV dynamics. RESULTS: The ICA algorithm has been used to generate different patterns of the HIV dynamics under different therapy conditions. The Kohonen Map has been used to eliminate redundant noise in each pattern to produce a reliable data set for the simulation phase. We show that under potent antiretroviral drugs, the value of the CD4+ cells in infected persons decreases gradually by about 11% every 100 days and the levels of the CD8+ cells increase gradually by about 2% every 100 days. AVAILABILITY: Executable code and data libraries are available by contacting the corresponding author. IMPLEMENTATION: Mathematica 4 has been used to simulate the suggested model. A Pentium III or higher platform is recommended.     

Fine-scale genetic mapping using independent component analysis

The aim of genetic mapping is to locate the loci responsible for specific traits such as complex diseases. These traits are normally caused by mutations at multiple loci of unknown locations and interactions. In this work, we model the biological system that relates DNA polymorphisms with complex traits as a linear mixing process. Given this model, we propose a new fine-scale genetic mapping method based on independent component analysis. The proposed method outputs both independent associated groups of SNPs in addition to specific associated SNPs with the phenotype. It is applied to a clinical data set for the Schizophrenia disease with 368 individuals and 42 SNPs. It is also applied to a simulation study to investigate in more depth its performance. The obtained results demonstrate the novel characteristics of the proposed method compared to other genetic mapping methods. Finally, we study the robustness of the proposed method with missing genotype values and limited sample sizes.     

Model-free functional MRI analysis using topographic independent component analysis

Data-driven fMRI analysis techniques include independent component analysis (ICA) and different types of clustering in the temporal domain. Since each of these methods has its particular strengths, it is natural to look for an approach that unifies Kohonen's self-organizing map and ICA. This is given by the topographic independent component analysis. While achieved by a slight modification of the ICA model, it can be at the same time used to define a topographic order (clusters) between the components, and thus has the usual computational advantages associated with topographic maps. In this contribution, we can show that when applied to fMRI analysis it outperforms FastICA.     

Investigations into resting-state connectivity using independent component analysis

Inferring resting-state connectivity patterns from functional magnetic resonance imaging (fMRI) data is a challenging task for any analytical technique. In this paper, we review a probabilistic independent component analysis (PICA) approach, optimized for the analysis of fMRI data, and discuss the role which this exploratory technique can take in scientific investigations into the structure of these effects. We apply PICA to fMRI data acquired at rest, in order to characterize the spatio-temporal structure of such data, and demonstrate that this is an effective and robust tool for the identification of low-frequency resting-state patterns from data acquired at various different spatial and temporal resolutions. We show that these networks exhibit high spatial consistency across subjects and closely resemble discrete cortical functional networks such as visual cortical areas or sensory-motor cortex.     

Application of independent component analysis to microarrays

We apply linear and nonlinear independent component analysis (ICA) to project microarray data into statistically independent components that correspond to putative biological processes, and to cluster genes according to over- or under-expression in each component. We test the statistical significance of enrichment of gene annotations within clusters. ICA outperforms other leading methods, such as principal component analysis, k-means clustering and the Plaid model, in constructing functionally coherent clusters on microarray datasets from Saccharomyces cerevisiae, Caenorhabditis elegans and human.     

The influence of visual imagery experience on EEG spatial organization

The comparison of EEG spatial organization between groups of 23 students of graphic arts department ("professional" subjects) and 39 subjects of another specialization ("non-professional" subjects) was made in order to find EEG correlates of visual imagery experience. Changes in the spatial organization of biopotentials (spatial synchronization and spatial disorder, spectral power and coherence) were analyzed while subjects mentally composed visual images from two simple elements, right angle and oblique line. The total number of elements presented for the image composition increased with each subsequent task (in total, four tasks were presented) from the number adequate to simultaneous perception and conscious processing (less than 7 +/- 2) to a much higher number. Intergroup differences, especially, in the degree of the spatial disorder (non-linear processes), were most evident under conditions when the subjects operated with a greater number of elements (tasks 3 and 4). This parameter increased more rapidly in "professionals" than in "non-professionals". These changes were most pronounced in the right anterior cortex. In "non-professional" subjects, spatial synchronization (linear processes) increased in the right posterior area. In "professional" subjects, coherence and spectral power increased in a greater number of narrow EEG frequency subbands than in "non-professional" subjects. The findings suggest that the imagery performance in subjects with visual imagery experience involves complicated neurodynamic processes such as non-linear dynamics and numerous EEG spatial resonance systems.     

Gene expression data classification using consensus independent component analysis

We propose a new method for tumor classification from gene expression data, which mainly contains three steps. Firstly, the original DNA microarray gene expression data are modeled by independent component analysis （ICA）. Secondly, the most discriminant eigenassays extracted by ICA are selected by the sequential floating forward selection technique. Finally, support vector machine is used to classify the modeling data. To show the validity of the proposed method, we applied it to classify three DNA microarray datasets involving various human normal and tumor tissue samples. The experimental results show that the method is efficient and feasible.     

Nonlinear dynamics of the EEG separated by independent component analysis after sound and light stimulation

 The electroencephalogram (EEG) is a multiscaled signal consisting of several time-series components each with different dominant frequency ranges and different origins. Nonlinear measures of the EEG reflect the complexity of the overall EEG, but not of each component in it. The aim of this study is to examine effect of the sound and light (SL) stimulation on the complexity of each component of the EEG. We used independent component analysis to obtain independent components of the EEG. The first positive Lyapunov exponent (L1) was estimated as a nonlinear measure of complexity. Twelve subjects were administered photic and auditory stimuli with a frequency of 10 Hz, which corresponded to the alpha frequency of the EEG, by a sound and light entrainment device. We compared the L1 values of the EEGs and their independent components between baseline and after the SL stimulation. We detected that the L1 values of the EEG decreased after the SL stimulation in all channels except C3 and F4, indicating that the complexity of the EEG decreased. We showed that alpha components increased in proportion but decreased in complexity after the SL stimulation. The beta independent components were found to decrease in proportion and complexity. These results suggest that decreased complexity of the EEG after the SL stimulation may be principally caused by decreased complexity and increased proportion of the alpha independent components. We showed also that theta components increased in complexity after the SL stimulation. We propose that nonlinear dynamical analysis combined with independent component analysis may be helpful in understanding the temporal characteristics of the EEG, which cannot be detected by conventional linear or nonlinear methods. Received: 12 March 2001 / Accepted in revised form: 27 November 2001     

Simulation of retinal ganglion cell response using fast independent component analysis

Advances in neurobiology suggest that neuronal response of the primary visual cortex to natural stimuli may be attributed to sparse approximation of images, encoding stimuli to activate specific neurons although the underlying mechanisms are still unclear. The responses of retinal ganglion cells (RGCs) to natural and random checkerboard stimuli were simulated using fast independent component analysis. The neuronal response to stimuli was measured using kurtosis and Treves–Rolls sparseness, and the kurtosis, lifetime and population sparseness were analyzed. RGCs exhibited significant lifetime sparseness in response to natural stimuli and random checkerboard stimuli. About 65 and 72% of RGCs do not fire all the time in response to natural and random checkerboard stimuli, respectively. Both kurtosis of single neurons and lifetime response of single neurons values were larger in the case of natural than in random checkerboard stimuli. The population of RGCs fire much less in response to random checkerboard stimuli than natural stimuli. However, kurtosis of population sparseness and population response of the entire neurons were larger with natural than random checkerboard stimuli. RGCs fire more sparsely in response to natural stimuli. Individual neurons fire at a low rate, while the occasional “burst” of neuronal population transmits information efficiently.     

Bilinear neural network with 3-D attention for brain decoding of motor imagery movements from the human EEG

Deep learning has achieved great success in areas such as computer vision and natural language processing. In the past, some work used convolutional networks to process EEG signals and reached or exceeded traditional machine learning methods. We propose a novel network structure and call it QNet. It contains a newly designed attention module: 3D-AM, which is used to learn the attention weights of EEG channels, time points, and feature maps. It provides a way to automatically learn the electrode and time selection. QNet uses a dual branch structure to fuse bilinear vectors for classification. It performs four, three, and two classes on the EEG Motor Movement/Imagery Dataset. The average cross-validation accuracy of 65.82%, 74.75%, and 82.88% was obtained, which are 7.24%, 4.93%, and 2.45% outperforms than the state-of-the-art, respectively. The article also visualizes the attention weights learned by QNet and shows its possible application for electrode channel selection.     

Spectral clustering independent component analysis for tissue classification from brain MRI

A spectral angle based feature extraction method, Spectral Clustering Independent Component Analysis (SC-ICA), is proposed in this work to improve the brain tissue classification from Magnetic Resonance Images (MRI). SC-ICA provides equal priority to global and local features; thereby it tries to resolve the inefficiency of conventional approaches in abnormal tissue extraction. First, input multispectral MRI is divided into different clusters by a spectral distance based clustering. Then, Independent Component Analysis (ICA) is applied on the clustered data, in conjunction with Support Vector Machines (SVM) for brain tissue analysis. Normal and abnormal datasets, consisting of real and synthetic T1-weighted, T2-weighted and proton density/fluid-attenuated inversion recovery images, were used to evaluate the performance of the new method. Comparative analysis with ICA based SVM and other conventional classifiers established the stability and efficiency of SC-ICA based classification, especially in reproduction of small abnormalities. Clinical abnormal case analysis demonstrated it through the highest Tanimoto Index/accuracy values, 0.75/98.8%, observed against ICA based SVM results, 0.17/96.1%, for reproduced lesions. Experimental results recommend the proposed method as a promising approach in clinical and pathological studies of brain diseases.     

Bayesian independent component analysis recovers pathway signatures from blood metabolomics data

Interpreting the complex interplay of metabolites in heterogeneous biosamples still poses a challenging task. In this study, we propose independent component analysis (ICA) as a multivariate analysis tool for the interpretation of large-scale metabolomics data. In particular, we employ a Bayesian ICA method based on a mean-field approach, which allows us to statistically infer the number of independent components to be reconstructed. The advantage of ICA over correlation-based methods like principal component analysis (PCA) is the utilization of higher order statistical dependencies, which not only yield additional information but also allow a more meaningful representation of the data with fewer components. We performed the described ICA approach on a large-scale metabolomics data set of human serum samples, comprising a total of 1764 study probands with 218 measured metabolites. Inspecting the source matrix of statistically independent metabolite profiles using a weighted enrichment algorithm, we observe strong enrichment of specific metabolic pathways in all components. This includes signatures from amino acid metabolism, energy-related processes, carbohydrate metabolism, and lipid metabolism. Our results imply that the human blood metabolome is composed of a distinct set of overlaying, statistically independent signals. ICA furthermore produces a mixing matrix, describing the strength of each independent component for each of the study probands. Correlating these values with plasma high-density lipoprotein (HDL) levels, we establish a novel association between HDL plasma levels and the branched-chain amino acid pathway. We conclude that the Bayesian ICA methodology has the power and flexibility to replace many of the nowadays common PCA and clustering-based analyses common in the research field.     

Integrated analysis of gene expression and copy number data on gene shaving using independent component analysis

DNA microarray gene expression and microarray-based comparative genomic hybridization (aCGH) have been widely used for biomedical discovery. Because of the large number of genes and the complex nature of biological networks, various analysis methods have been proposed. One such method is "gene shaving," a procedure which identifies subsets of the genes with coherent expression patterns and large variation across samples. Since combining genomic information from multiple sources can improve classification and prediction of diseases, in this paper we proposed a new method, "ICA gene shaving" (ICA, independent component analysis), for jointly analyzing gene expression and copy number data. First we used ICA to analyze joint measurements, gene expression and copy number, of a biological system and project the data onto statistically independent biological processes. Next, we used these results to identify patterns of variation in the data and then applied an iterative shaving method. We investigated the properties of our proposed method by analyzing both simulated and real data. We demonstrated that the robustness of our method to noise using simulated data. Using breast cancer data, we showed that our method is superior to the Generalized Singular Value Decomposition (GSVD) gene shaving method for identifying genes associated with breast cancer.     

Resting‐state networks of the neonate brain identified using independent component analysis

Resting‐state functional magnetic resonance imaging (rs‐fMRI) has been successfully used to probe the intrinsic functional organization of the brain and to study brain development. Here, we implemented a combination of individual and group independent component analysis (ICA) of FSL on a 6‐min resting‐state data set acquired from 21 naturally sleeping term‐born (age 26 ± 6.7 d), healthy neonates to investigate the emerging functional resting‐state networks (RSNs). In line with the previous literature, we found evidence of sensorimotor, auditory/language, visual, cerebellar, thalmic, parietal, prefrontal, anterior cingulate as well as dorsal and ventral aspects of the default‐mode‐network. Additionally, we identified RSNs in frontal, parietal, and temporal regions that have not been previously described in this age group and correspond to the canonical RSNs established in adults. Importantly, we found that careful ICA‐based denoising of fMRI data increased the number of networks identified with group‐ICA, whereas the degree of spatial smoothing did not change the number of identified networks. Our results show that the infant brain has an established set of RSNs soon after birth.     

Generalized optimal spatial filtering using a kernel approach with application to EEG classification

Common spatial patterns (CSP) has been widely used for finding the linear spatial filters which are able to extract the discriminative brain activities between two different mental tasks. However, the CSP is difficult to capture the nonlinearly clustered structure from the non-stationary EEG signals. To relax the presumption of strictly linear patterns in the CSP, in this paper, a generalized CSP (GCSP) based on generalized singular value decomposition (GSVD) and kernel method is proposed. Our method is able to find the nonlinear spatial filters which are formulated in the feature space defined by a nonlinear mapping through kernel functions. Furthermore, in order to overcome the overfitting problem, the regularized GCSP is developed by adding the regularized parameters. The experimental results demonstrate that our method is an effective nonlinear spatial filtering method.     

Estimation of leaf water potential by thermal imagery and spatial analysis

Canopy temperature has long been recognized as an indicator of plant water status and as a potential tool for irrigation scheduling. In the present study, the potential of using thermal images for an in-field estimation of the water status of cotton under a range of irrigation regimes was investigated. Thermal images were taken with a radiometric infrared video camera. Specific leaves that appeared in the camera field of view were sampled, their LWP was measured and their temperature was calculated from the images. Regression models were built in order to predict LWP according to the crop canopy temperature and to the empirical formulation of the crop water stress index (CWSI). Statistical analysis revealed that the relationship between CWSI and LWP was more stable and had slightly higher correlation coefficients than that between canopy temperature and LWP. The regression models of LWP against CWSI and against leaf temperatures were used to create LWP maps. The classified LWP maps showed that there was spatial variability in each treatment, some of which may be attributed to the difference between sunlit and shaded leaves. The distribution of LWP in the maps showed that irrigation treatments were better distinguished from each other when the maps were calculated from CWSI than from leaf temperature alone. Furthermore, the inclusion of the spatial pattern in the classification enhanced the differences between the treatments and was better matched to irrigation amounts. Optimal determination of the water status from thermal images should be based on an overall view of the physical status as well as on the analysis of the spatial structure. Future study will involve investigating the robustness of the models and the potential of using water status maps, derived from aerial thermal images, for irrigation scheduling and variable management in commercial fields.