Dataset Information

Independent Component Analysis Involving Autocorrelated Sources With an Application to Functional Magnetic Resonance Imaging.

ABSTRACT: Independent component analysis (ICA) is an effective data-driven method for blind source separation. It has been successfully applied to separate source signals of interest from their mixtures. Most existing ICA procedures are carried out by relying solely on the estimation of the marginal density functions, either parametrically or nonparametrically. In many applications, correlation structures within each source also play an important role besides the marginal distributions. One important example is functional magnetic resonance imaging (fMRI) analysis where the brain-function-related signals are temporally correlated. In this article, we consider a novel approach to ICA that fully exploits the correlation structures within the source signals. Specifically, we propose to estimate the spectral density functions of the source signals instead of their marginal density functions. This is made possible by virtue of the intrinsic relationship between the (unobserved) sources and the (observed) mixed signals. Our methodology is described and implemented using spectral density functions from frequently used time series models such as autoregressive moving average (ARMA) processes. The time series parameters and the mixing matrix are estimated via maximizing the Whittle likelihood function. We illustrate the performance of the proposed method through extensive simulation studies and a real fMRI application. The numerical results indicate that our approach outperforms several popular methods including the most widely used fastICA algorithm. This article has supplementary material online.

SUBMITTER: Lee S

PROVIDER: S-EPMC4979078 | biostudies-literature | 2011

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Independent Component Analysis Involving Autocorrelated Sources With an Application to Functional Magnetic Resonance Imaging.

Lee Seonjoo S Shen Haipeng H Truong Young Y Lewis Mechelle M Huang Xuemei X

Journal of the American Statistical Association 20110101 495

Independent component analysis (ICA) is an effective data-driven method for blind source separation. It has been successfully applied to separate source signals of interest from their mixtures. Most existing ICA procedures are carried out by relying solely on the estimation of the marginal density functions, either parametrically or nonparametrically. In many applications, correlation structures within each source also play an important role besides the marginal distributions. One important exam ...[more]

PMID: 27524847

Similar Datasets

Project description:ObjectivesThe aim of this study was to develop and verify a method to obtain good temporal resolution T2-weighted 4-dimensional (4D-T2w) magnetic resonance imaging (MRI) by using motion information from T1-weighted 4D (4D-T1w) MRI, to support treatment planning in MR-guided radiotherapy.Materials and methodsTen patients with primary non-small cell lung cancer were scanned at 1.5 T axially with a volumetric T2-weighted turbo spin echo sequence gated to exhalation and a volumetric T1-weighted stack-of-stars spoiled gradient echo sequence with golden angle spacing acquired in free breathing. From the latter, 20 respiratory phases were reconstructed using the recently developed 4D joint MoCo-HDTV algorithm based on the self-gating signal obtained from the k-space center. Motion vector fields describing the respiratory cycle were obtained by deformable image registration between the respiratory phases and projected onto the T2-weighted image volume. The resulting 4D-T2w volumes were verified against the 4D-T1w volumes: an edge-detection method was used to measure the diaphragm positions; the locations of anatomical landmarks delineated by a radiation oncologist were compared and normalized mutual information was calculated to evaluate volumetric image similarity.ResultsHigh-resolution 4D-T2w MRI was obtained. Respiratory motion was preserved on calculated 4D-T2w MRI, with median diaphragm positions being consistent with less than 6.6 mm (2 voxels) for all patients and less than 3.3 mm (1 voxel) for 9 of 10 patients. Geometrical positions were coherent between 4D-T1w and 4D-T2w MRI as Euclidean distances between all corresponding anatomical landmarks agreed to within 7.6 mm (Euclidean distance of 2 voxels) and were below 3.8 mm (Euclidean distance of 1 voxel) for 355 of 470 pairs of anatomical landmarks. Volumetric image similarity was commensurate between 4D-T1w and 4D-T2w MRI, as mean percentage differences in normalized mutual information (calculated over all respiratory phases and patients), between corresponding respiratory phases of 4D-T1w and 4D-T2w MRI and the tie-phase of 4D-T1w and 3-dimensional T2w MRI, were consistent to 0.41% ± 0.37%. Four-dimensional T2w MRI displayed tumor extent, structure, and position more clearly than corresponding 4D-T1w MRI, especially when mobile tumor sites were adjacent to organs at risk.ConclusionsA methodology to obtain 4D-T2w MRI that retrospectively applies the motion information from 4D-T1w MRI to 3-dimensional T2w MRI was developed and verified. Four-dimensional T2w MRI can assist clinicians in delineating mobile lesions that are difficult to define on 4D-T1w MRI, because of poor tumor-tissue contrast.

Project description:The General Linear Model (GLM) used in task-fMRI relates activated brain areas to extrinsic task conditions. The translation of resulting neural activation into a hemodynamic response is commonly approximated with a linear convolution model using a hemodynamic response function (HRF). There are two major limitations in GLM analysis. Firstly, the GLM assumes that neural activation is either on or off and matches the exact stimulus duration in the corresponding task timings. Secondly, brain networks observed in resting-state fMRI experiments present also during task experiments, but the GLM approach models these task-unrelated brain activity as noise. A novel kernel matrix factorization approach, called hemodynamic matrix factorization (HMF), is therefore proposed that addresses both limitations by assuming that task-related and task-unrelated brain activity can be modeled with the same convolution model as in GLM analysis. By contrast to the GLM, the proposed HMF is a blind source separation (BSS) technique, which decomposes fMRI data into modes. Each mode comprises of a neural activation time course and a spatial mapping. Two versions of HMF are proposed in which the neural activation time course of each mode is convolved with either the canonical HRF or predetermined subject-specific HRFs. Firstly, HMF with the canonical HRF is applied to two open-source cohorts. These cohorts comprise of several task experiments including motor, incidental memory, spatial coherence discrimination, verbal discrimination task and a very short localization task, engaging multiple parts of the eloquent cortex. HMF modes were obtained whose neural activation time course followed original task timings and whose corresponding spatial map matched cortical areas known to be involved in the respective task processing. Secondly, the alignment of these neural activation time courses to task timings were further improved by replacing the canonical HRF with subject-specific HRFs during HMF mode computation. In addition to task-related modes, HMF also produced seemingly task-unrelated modes whose spatial maps matched known resting-state networks. The validity of a fMRI task experiment relies on the assumption that the exposure to a stimulus for a given time causes an imminent increase in neural activation of equal duration. The proposed HMF is an attempt to falsify this assumption and allows to identify subject task participation that does not comply with the experiment instructions.

Dataset Information

Independent Component Analysis Involving Autocorrelated Sources With an Application to Functional Magnetic Resonance Imaging.

Publications

Independent Component Analysis Involving Autocorrelated Sources With an Application to Functional Magnetic Resonance Imaging.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure