@article {819, title = {Analytic programming with FMRI data: a quick-start guide for statisticians using R.}, journal = {PLoS One}, volume = {9}, year = {2014}, month = {2014}, pages = {e89470}, abstract = {

Functional magnetic resonance imaging (fMRI) is a thriving field that plays an important role in medical imaging analysis, biological and neuroscience research and practice. This manuscript gives a didactic introduction to the statistical analysis of fMRI data using the R project, along with the relevant R code. The goal is to give statisticians who would like to pursue research in this area a quick tutorial for programming with fMRI data. References of relevant packages and papers are provided for those interested in more advanced analysis.

}, keywords = {Brain, brain mapping, Humans, Magnetic Resonance Imaging, Programming Languages}, issn = {1932-6203}, doi = {10.1371/journal.pone.0089470}, author = {Eloyan, Ani and Li, Shanshan and Muschelli, John and Pekar, Jim J and Mostofsky, Stewart H and Caffo, Brian S} } @article {822, title = {An evaluation of independent component analyses with an application to resting-state fMRI.}, journal = {Biometrics}, volume = {70}, year = {2014}, month = {2014 Mar}, pages = {224-36}, abstract = {

We examine differences between independent component analyses (ICAs) arising from different assumptions, measures of dependence, and starting points of the algorithms. ICA is a popular method with diverse applications including artifact removal in electrophysiology data, feature extraction in microarray data, and identifying brain networks in functional magnetic resonance imaging (fMRI). ICA can be viewed as a generalization of principal component analysis (PCA) that takes into account higher-order cross-correlations. Whereas the PCA solution is unique, there are many ICA methods-whose solutions may differ. Infomax, FastICA, and JADE are commonly applied to fMRI studies, with FastICA being arguably the most popular. Hastie and Tibshirani (2003) demonstrated that ProDenICA outperformed FastICA in simulations with two components. We introduce the application of ProDenICA to simulations with more components and to fMRI data. ProDenICA was more accurate in simulations, and we identified differences between biologically meaningful ICs from ProDenICA versus other methods in the fMRI analysis. ICA methods require nonconvex optimization, yet current practices do not recognize the importance of, nor adequately address sensitivity to, initial values. We found that local optima led to dramatically different estimates in both simulations and group ICA of fMRI, and we provide evidence that the global optimum from ProDenICA is the best estimate. We applied a modification of the Hungarian (Kuhn-Munkres) algorithm to match ICs from multiple estimates, thereby gaining novel insights into how brain networks vary in their sensitivity to initial values and ICA method.

}, keywords = {Adolescent, Algorithms, Attention Deficit Disorder with Hyperactivity, brain mapping, Child, Child, Preschool, Computer Simulation, Humans, Magnetic Resonance Imaging, Models, Statistical, Principal Component Analysis}, issn = {1541-0420}, doi = {10.1111/biom.12111}, author = {Risk, Benjamin B and Matteson, David S and Ruppert, David and Eloyan, Ani and Caffo, Brian S} } @article {790, title = {Health effects of lesion localization in multiple sclerosis: spatial registration and confounding adjustment.}, journal = {PLoS One}, volume = {9}, year = {2014}, month = {2014}, pages = {e107263}, abstract = {

Brain lesion localization in multiple sclerosis (MS) is thought to be associated with the type and severity of adverse health effects. However, several factors hinder statistical analyses of such associations using large MRI datasets: 1) spatial registration algorithms developed for healthy individuals may be less effective on diseased brains and lead to different spatial distributions of lesions; 2) interpretation of results requires the careful selection of confounders; and 3) most approaches have focused on voxel-wise regression approaches. In this paper, we evaluated the performance of five registration algorithms and observed that conclusions regarding lesion localization can vary substantially with the choice of registration algorithm. Methods for dealing with confounding factors due to differences in disease duration and local lesion volume are introduced. Voxel-wise regression is then extended by the introduction of a metric that measures the distance between a patient-specific lesion mask and the population prevalence map.

}, issn = {1932-6203}, doi = {10.1371/journal.pone.0107263}, author = {Eloyan, Ani and Shou, Haochang and Shinohara, Russell T and Sweeney, Elizabeth M and Nebel, Mary Beth and Cuzzocreo, Jennifer L and Calabresi, Peter A and Reich, Daniel S and Lindquist, Martin A and Crainiceanu, Ciprian M} } @article {792, title = {Shrinkage prediction of seed-voxel brain connectivity using resting state fMRI.}, journal = {Neuroimage}, volume = {102 Pt 2}, year = {2014}, month = {2014 Nov 15}, pages = {938-44}, abstract = {

Resting-state functional magnetic resonance imaging (rs-fMRI) is used to investigate synchronous activations in spatially distinct regions of the brain, which are thought to reflect functional systems supporting cognitive processes. Analyses are often performed using seed-based correlation analysis, allowing researchers to explore functional connectivity between data in a seed region and the rest of the brain. Using scan-rescan rs-fMRI data, we investigate how well the subject-specific seed-based correlation map from the second replication of the study can be predicted using data from the first replication. We show that one can dramatically improve prediction of subject-specific connectivity by borrowing strength from the group correlation map computed using all other subjects in the study. Even more surprisingly, we found that the group correlation map provided a better prediction of a subject{\textquoteright}s connectivity than the individual{\textquoteright}s own data. While further discussion and experimentation are required to understand how this can be used in practice, results indicate that shrinkage-based methods that borrow strength from the population mean should play a role in rs-fMRI data analysis.

}, issn = {1095-9572}, doi = {10.1016/j.neuroimage.2014.05.043}, author = {Shou, Haochang and Eloyan, Ani and Nebel, Mary Beth and Mejia, Amanda and Pekar, James J and Mostofsky, Stewart and Caffo, Brian and Lindquist, Martin A and Crainiceanu, Ciprian M} } @article {812, title = {Likelihood-based population independent component analysis.}, journal = {Biostatistics}, volume = {14}, year = {2013}, month = {2013 Jul}, pages = {514-27}, abstract = {

Independent component analysis (ICA) is a widely used technique for blind source separation, used heavily in several scientific research areas including acoustics, electrophysiology, and functional neuroimaging. We propose a scalable two-stage iterative true group ICA methodology for analyzing population level functional magnetic resonance imaging (fMRI) data where the number of subjects is very large. The method is based on likelihood estimators of the underlying source densities and the mixing matrix. As opposed to many commonly used group ICA algorithms, the proposed method does not require significant data reduction by a 2-fold singular value decomposition. In addition, the method can be applied to a large group of subjects since the memory requirements are not restrictive. The performance of our approach is compared with a commonly used group ICA algorithm via simulation studies. Furthermore, the proposed method is applied to a large collection of resting state fMRI datasets. The results show that established brain networks are well recovered by the proposed algorithm.

}, keywords = {Adult, Algorithms, Attention Deficit Disorder with Hyperactivity, Biostatistics, Brain, Child, Computer Simulation, Databases, Factual, Humans, Imaging, Three-Dimensional, Likelihood Functions, Magnetic Resonance Imaging, Principal Component Analysis, Statistics, Nonparametric}, issn = {1468-4357}, doi = {10.1093/biostatistics/kxs055}, author = {Eloyan, Ani and Crainiceanu, Ciprian M and Caffo, Brian S} } @article {827, title = {Analysis of group ICA-based connectivity measures from fMRI: application to Alzheimer{\textquoteright}s disease.}, journal = {PLoS One}, volume = {7}, year = {2012}, month = {2012}, pages = {e49340}, abstract = {

Functional magnetic resonance imaging (fMRI) is a powerful tool for the in vivo study of the pathophysiology of brain disorders and disease. In this manuscript, we propose an analysis stream for fMRI functional connectivity data and apply it to a novel study of Alzheimer{\textquoteright}s disease. In the first stage, spatial independent component analysis is applied to group fMRI data to obtain common brain networks (spatial maps) and subject-specific mixing matrices (time courses). In the second stage, functional principal component analysis is utilized to decompose the mixing matrices into population-level eigenvectors and subject-specific loadings. Inference is performed using permutation-based exact logistic regression for matched pairs data. The method is applied to a novel fMRI study of Alzheimer{\textquoteright}s disease risk under a verbal paired associates task. We found empirical evidence of alternative ICA-based metrics of connectivity when comparing subjects evidencing mild cognitive impairment relative to carefully matched controls.

}, keywords = {Aged, Algorithms, Alzheimer Disease, Brain, brain mapping, Case-Control Studies, Female, Humans, Image Processing, Computer-Assisted, Magnetic Resonance Imaging, Male, Middle Aged, Mild Cognitive Impairment, Nerve Net, Principal Component Analysis, Task Performance and Analysis, Time Factors}, issn = {1932-6203}, doi = {10.1371/journal.pone.0049340}, author = {Li, Shanshan and Eloyan, Ani and Joel, Suresh and Mostofsky, Stewart and Pekar, James and Bassett, Susan Spear and Caffo, Brian} }