@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 {810, title = {Combining hidden Markov models for comparing the dynamics of multiple sleep electroencephalograms.}, journal = {Stat Med}, volume = {32}, year = {2013}, month = {2013 Aug 30}, pages = {3342-56}, abstract = {

In this manuscript, we consider methods for the analysis of populations of electroencephalogram signals during sleep for the study of sleep disorders using hidden Markov models (HMMs). Notably, we propose an easily implemented method for simultaneously modeling multiple time series that involve large amounts of data. We apply these methods to study sleep-disordered breathing (SDB) in the Sleep Heart Health Study (SHHS), a landmark study of SDB and cardiovascular consequences. We use the entire, longitudinally collected, SHHS cohort to develop HMM population parameters, which we then apply to obtain subject-specific Markovian predictions. From these predictions, we create several indices of interest, such as transition frequencies between latent states. Our HMM analysis of electroencephalogram signals uncovers interesting findings regarding differences in brain activity during sleep between those with and without SDB. These findings include stability of the percent time spent in HMM latent states across matched diseased and non-diseased groups and differences in the rate of transitioning.

}, keywords = {Data Interpretation, Statistical, Electroencephalography, Female, Humans, Longitudinal Studies, Male, Markov Chains, Middle Aged, Models, Statistical, Sleep Apnea Syndromes}, issn = {1097-0258}, doi = {10.1002/sim.5747}, author = {Langrock, Roland and Swihart, Bruce J and Caffo, Brian S and Punjabi, Naresh M and Crainiceanu, Ciprian M} } @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 {814, title = {Longitudinal scalar-on-functions regression with application to tractography data.}, journal = {Biostatistics}, volume = {14}, year = {2013}, month = {2013 Jul}, pages = {447-61}, abstract = {

We propose a class of estimation techniques for scalar-on-function regression where both outcomes and functional predictors may be observed at multiple visits. Our methods are motivated by a longitudinal brain diffusion tensor imaging tractography study. One of the study{\textquoteright}s primary goals is to evaluate the contemporaneous association between human function and brain imaging over time. The complexity of the study requires the development of methods that can simultaneously incorporate: (1) multiple functional (and scalar) regressors; (2) longitudinal outcome and predictor measurements per patient; (3) Gaussian or non-Gaussian outcomes; and (4) missing values within functional predictors. We propose two versions of a new method, longitudinal functional principal components regression (PCR). These methods extend the well-known functional PCR and allow for different effects of subject-specific trends in curves and of visit-specific deviations from that trend. The new methods are compared with existing approaches, and the most promising techniques are used for analyzing the tractography data.

}, keywords = {Anisotropy, Biostatistics, Brain, Diffusion Tensor Imaging, Humans, Linear Models, Models, Statistical, Multiple Sclerosis, Principal Component Analysis, Regression Analysis}, issn = {1468-4357}, doi = {10.1093/biostatistics/kxs051}, author = {Gertheiss, Jan and Goldsmith, Jeff and Crainiceanu, Ciprian and Greven, Sonja} }