@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 {843, title = {Cloak and DAG: a response to the comments on our comment.}, journal = {Neuroimage}, volume = {76}, year = {2013}, month = {2013 Aug 1}, pages = {446-9}, abstract = {

Our original comment (Lindquist and Sobel, 2011) made explicit the types of assumptions neuroimaging researchers are making when directed graphical models (DGMs), which include certain types of structural equation models (SEMs), are used to estimate causal effects. When these assumptions, which many researchers are not aware of, are not met, parameters of these models should not be interpreted as effects. Thus it is imperative that neuroimaging researchers interested in issues involving causation, for example, effective connectivity, consider the plausibility of these assumptions for their particular problem before using SEMs. In cases where these additional assumptions are not met, researchers may be able to use other methods and/or design experimental studies where the use of unrealistic assumptions can be avoided. Pearl does not disagree with anything we stated. However, he takes exception to our use of potential outcomes{\textquoteright} notation, which is the standard notation used in the statistical literature on causal inference, and his comment is devoted to promoting his alternative conventions. Glymour{\textquoteright}s comment is based on three claims that he inappropriately attributes to us. Glymour is also more optimistic than us about the potential of using directed graphical models (DGMs) to discover causal relations in neuroimaging research; we briefly address this issue toward the end of our rejoinder.

}, keywords = {Artifacts, Brain, Computer Simulation, Humans, Image Interpretation, Computer-Assisted, Meta-Analysis as Topic}, issn = {1095-9572}, doi = {10.1016/j.neuroimage.2011.11.027}, author = {Lindquist, Martin A and Sobel, Michael E} } @article {684, title = {Corrected confidence bands for functional data using principal components.}, journal = {Biometrics}, volume = {69}, year = {2013}, month = {2013 Mar}, pages = {41-51}, abstract = {

Functional principal components (FPC) analysis is widely used to decompose and express functional observations. Curve estimates implicitly condition on basis functions and other quantities derived from FPC decompositions; however these objects are unknown in practice. In this article, we propose a method for obtaining correct curve estimates by accounting for uncertainty in FPC decompositions. Additionally, pointwise and simultaneous confidence intervals that account for both model- and decomposition-based variability are constructed. Standard mixed model representations of functional expansions are used to construct curve estimates and variances conditional on a specific decomposition. Iterated expectation and variance formulas combine model-based conditional estimates across the distribution of decompositions. A bootstrap procedure is implemented to understand the uncertainty in principal component decomposition quantities. Our method compares favorably to competing approaches in simulation studies that include both densely and sparsely observed functions. We apply our method to sparse observations of CD4 cell counts and to dense white-matter tract profiles. Code for the analyses and simulations is publicly available, and our method is implemented in the R package refund on CRAN.

}, keywords = {Brain, CD4 Lymphocyte Count, Computer Simulation, Confidence Intervals, HIV, HIV Infections, Humans, Magnetic Resonance Imaging, Models, Statistical, Multiple Sclerosis, Principal Component Analysis}, issn = {1541-0420}, doi = {10.1111/j.1541-0420.2012.01808.x}, author = {Goldsmith, J and Greven, S and Crainiceanu, C} } @article {824, title = {A decision-theory approach to interpretable set analysis for high-dimensional data.}, journal = {Biometrics}, volume = {69}, year = {2013}, month = {2013 Sep}, pages = {614-23}, abstract = {

A key problem in high-dimensional significance analysis is to find pre-defined sets that show enrichment for a statistical signal of interest; the classic example is the enrichment of gene sets for differentially expressed genes. Here, we propose a new decision-theory approach to the analysis of gene sets which focuses on estimating the fraction of non-null variables in a set. We introduce the idea of "atoms," non-overlapping sets based on the original pre-defined set annotations. Our approach focuses on finding the union of atoms that minimizes a weighted average of the number of false discoveries and missed discoveries. We introduce a new false discovery rate for sets, called the atomic false discovery rate (afdr), and prove that the optimal estimator in our decision-theory framework is to threshold the afdr. These results provide a coherent and interpretable framework for the analysis of sets that addresses the key issues of overlapping annotations and difficulty in interpreting p values in both competitive and self-contained tests. We illustrate our method and compare it to a popular existing method using simulated examples, as well as gene-set and brain ROI data analyses.

}, keywords = {Algorithms, Bayes Theorem, Biometry, Brain, Computer Simulation, Data Interpretation, Statistical, Decision Theory, Functional Neuroimaging, Gene Expression Profiling, Genomics, Humans, Magnetic Resonance Imaging, Models, Statistical, Oligonucleotide Array Sequence Analysis}, issn = {1541-0420}, doi = {10.1111/biom.12060}, author = {Boca, Simina M and Bravo, H{\'e}ctor C{\'e}orrada and Caffo, Brian and Leek, Jeffrey T and Parmigiani, Giovanni} } @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 {803, title = {Quantifying the reliability of image replication studies: the image intraclass correlation coefficient (I2C2).}, journal = {Cogn Affect Behav Neurosci}, volume = {13}, year = {2013}, month = {2013 Dec}, pages = {714-24}, abstract = {

This article proposes the image intraclass correlation (I2C2) coefficient as a global measure of reliability for imaging studies. The I2C2 generalizes the classic intraclass correlation (ICC) coefficient to the case when the data of interest are images, thereby providing a measure that is both intuitive and convenient. Drawing a connection with classical measurement error models for replication experiments, the I2C2 can be computed quickly, even in high-dimensional imaging studies. A nonparametric bootstrap procedure is introduced to quantify the variability of the I2C2 estimator. Furthermore, a Monte Carlo permutation is utilized to test reproducibility versus a zero I2C2, representing complete lack of reproducibility. Methodologies are applied to three replication studies arising from different brain imaging modalities and settings: regional analysis of volumes in normalized space imaging for characterizing brain morphology, seed-voxel brain activation maps based on resting-state functional magnetic resonance imaging (fMRI), and fractional anisotropy in an area surrounding the corpus callosum via diffusion tensor imaging. Notably, resting-state fMRI brain activation maps are found to have low reliability, ranging from .2 to .4. Software and data are available to provide easy access to the proposed methods.

}, keywords = {Adult, Brain, brain mapping, Computer Simulation, Female, Humans, Male, Models, Biological, Neuroimaging, Reproducibility of Results, Statistics as Topic}, issn = {1531-135X}, doi = {10.3758/s13415-013-0196-0}, author = {Shou, H and Eloyan, A and Lee, S and Zipunnikov, V and Crainiceanu, A N and Nebel, N B and Caffo, B and Lindquist, M A and Crainiceanu, C M} } @article {842, title = {Dynamic connectivity regression: determining state-related changes in brain connectivity.}, journal = {Neuroimage}, volume = {61}, year = {2012}, month = {2012 Jul 16}, pages = {907-20}, abstract = {

Most statistical analyses of fMRI data assume that the nature, timing and duration of the psychological processes being studied are known. However, often it is hard to specify this information a priori. In this work we introduce a data-driven technique for partitioning the experimental time course into distinct temporal intervals with different multivariate functional connectivity patterns between a set of regions of interest (ROIs). The technique, called Dynamic Connectivity Regression (DCR), detects temporal change points in functional connectivity and estimates a graph, or set of relationships between ROIs, for data in the temporal partition that falls between pairs of change points. Hence, DCR allows for estimation of both the time of change in connectivity and the connectivity graph for each partition, without requiring prior knowledge of the nature of the experimental design. Permutation and bootstrapping methods are used to perform inference on the change points. The method is applied to various simulated data sets as well as to an fMRI data set from a study (N=26) of a state anxiety induction using a socially evaluative threat challenge. The results illustrate the method{\textquoteright}s ability to observe how the networks between different brain regions changed with subjects{\textquoteright} emotional state.

}, keywords = {Algorithms, Brain, brain mapping, Computer Simulation, Emotions, Humans, Image Interpretation, Computer-Assisted, Magnetic Resonance Imaging, Neural Pathways}, issn = {1095-9572}, doi = {10.1016/j.neuroimage.2012.03.070}, author = {Cribben, Ivor and Haraldsdottir, Ragnheidur and Atlas, Lauren Y and Wager, Tor D and Lindquist, Martin A} } @article {845, title = {Graphical models, potential outcomes and causal inference: comment on Ramsey, Spirtes and Glymour.}, journal = {Neuroimage}, volume = {57}, year = {2011}, month = {2011 Jul 15}, pages = {334-6}, abstract = {

Ramsey, Spirtes and Glymour (RSG) critique a method proposed by Neumann et al. (2010) for the discovery of functional networks from fMRI meta-analysis data. We concur with this critique, but are unconvinced that directed graphical models (DGMs) are generally useful for estimating causal effects. We express our reservations using the "potential outcomes" framework for causal inference widely used in statistics.

}, keywords = {Artifacts, Brain, Computer Simulation, Humans, Image Interpretation, Computer-Assisted, Meta-Analysis as Topic}, issn = {1095-9572}, doi = {10.1016/j.neuroimage.2010.10.020}, author = {Lindquist, Martin A and Sobel, Michael E} } @article {Paninski2010, title = {A new look at state-space models for neural data.}, journal = {Journal of computational neuroscience}, volume = {29}, number = {1-2}, year = {2010}, pages = {107{\textendash}26}, abstract = {

State space methods have proven indispensable in neural data analysis. However, common methods for performing inference in state-space models with non-Gaussian observations rely on certain approximations which are not always accurate. Here we review direct optimization methods that avoid these approximations, but that nonetheless retain the computational efficiency of the approximate methods. We discuss a variety of examples, applying these direct optimization techniques to problems in spike train smoothing, stimulus decoding, parameter estimation, and inference of synaptic properties. Along the way, we point out connections to some related standard statistical methods, including spline smoothing and isotonic regression. Finally, we note that the computational methods reviewed here do not in fact depend on the state-space setting at all; instead, the key property we are exploiting involves the bandedness of certain matrices. We close by discussing some applications of this more general point of view, including Markov chain Monte Carlo methods for neural decoding and efficient estimation of spatially-varying firing rates.

}, keywords = {1 introduction, Action Potentials, Action Potentials: physiology, Animals, Computer Simulation, for inference in state-space, forward-backward methods, hidden markov model, models, neural coding, Neurological, Neurons, Neurons: physiology, Retinal Ganglion Cells, Retinal Ganglion Cells: physiology, state-space models, Statistical, Synapses, Synapses: physiology, tridiagonal matrix}, issn = {1573-6873}, doi = {10.1007/s10827-009-0179-x}, url = {http://www.ncbi.nlm.nih.gov/pubmed/19649698}, author = {Paninski, Liam and Ahmadian, Yashar and Gil, Daniel and Shinsuke, Ferreira and Rahnama, Kamiar and Michael, Rad and Vogelstein, Joshua T and Wu, Wei and Ferreira, Daniel Gil and Koyama, Shinsuke and {Rahnama Rad}, Kamiar and Vidne, Michael} }