Improved estimation of subject-level functional connectivity using full and partial correlation with empirical Bayes shrinkage

&NA; Reliability of subject‐level resting‐state functional connectivity (FC) is determined in part by the statistical techniques employed in its estimation. Methods that pool information across subjects to inform estimation of subject‐level effects (e.g., Bayesian approaches) have been shown to enhance reliability of subject‐level FC. However, fully Bayesian approaches are computationally demanding, while empirical Bayesian approaches typically rely on using repeated measures to estimate the variance components in the model. Here, we avoid the need for repeated measures by proposing a novel measurement error model for FC describing the different sources of variance and error, which we use to perform empirical Bayes shrinkage of subject‐level FC towards the group average. In addition, since the traditional intra‐class correlation coefficient (ICC) is inappropriate for biased estimates, we propose a new reliability measure denoted the mean squared error intra‐class correlation coefficient (ICCMSE) to properly assess the reliability of the resulting (biased) estimates. We apply the proposed techniques to test‐retest resting‐state fMRI data on 461 subjects from the Human Connectome Project to estimate connectivity between 100 regions identified through independent components analysis (ICA). We consider both correlation and partial correlation as the measure of FC and assess the benefit of shrinkage for each measure, as well as the effects of scan duration. We find that shrinkage estimates of subject‐level FC exhibit substantially greater reliability than traditional estimates across various scan durations, even for the most reliable connections and regardless of connectivity measure. Additionally, we find partial correlation reliability to be highly sensitive to the choice of penalty term, and to be generally worse than that of full correlations except for certain connections and a narrow range of penalty values. This suggests that the penalty needs to be chosen carefully when using partial correlations. HighlightsEmpirical Bayes shrinkage methods for functional connectivity are proposed.A novel reliability measure analogous to intraclass correlation is proposed.Shrinkage significantly improves reliability of full and partial correlations.Partial correlation reliability is highly sensitive to ridge regression penalty.Partial correlation reliability is worse overall but better for some connections.

[1]  M. Chun,et al.  Functional connectome fingerprinting: Identifying individuals based on patterns of brain connectivity , 2015, Nature Neuroscience.

[2]  S. Choi,et al.  Individual parcellation of resting fMRI with a group functional connectivity prior , 2017, NeuroImage.

[3]  E. Wagenmakers,et al.  Scientific rigor and the art of motorcycle maintenance , 2014, Nature Biotechnology.

[4]  Martin A. Lindquist,et al.  Shrinkage prediction of seed-voxel brain connectivity using resting state fMRI , 2014, NeuroImage.

[5]  Bradley P. Carlin,et al.  Markov Chain Monte Carlo in Practice: A Roundtable Discussion , 1998 .

[6]  Gaël Varoquaux,et al.  Learning and comparing functional connectomes across subjects , 2013, NeuroImage.

[7]  B. Biswal,et al.  The resting brain: unconstrained yet reliable. , 2009, Cerebral cortex.

[8]  Dustin Scheinost,et al.  Influences on the Test–Retest Reliability of Functional Connectivity MRI and its Relationship with Behavioral Utility , 2017, Cerebral cortex.

[9]  Mark Jenkinson,et al.  The minimal preprocessing pipelines for the Human Connectome Project , 2013, NeuroImage.

[10]  Eswar Damaraju,et al.  Tracking whole-brain connectivity dynamics in the resting state. , 2014, Cerebral cortex.

[11]  Daniel P. Kennedy,et al.  The Autism Brain Imaging Data Exchange: Towards Large-Scale Evaluation of the Intrinsic Brain Architecture in Autism , 2013, Molecular Psychiatry.

[12]  Evan M. Gordon,et al.  Functional System and Areal Organization of a Highly Sampled Individual Human Brain , 2015, Neuron.

[13]  Kuncheng Li,et al.  Reliability correction for functional connectivity: Theory and implementation , 2015, Human brain mapping.

[14]  M. Chun,et al.  Functional connectome fingerprinting: Identifying individuals based on patterns of brain connectivity , 2015, Nature Neuroscience.

[15]  C. Stein,et al.  Estimation with Quadratic Loss , 1992 .

[16]  Ludovica Griffanti,et al.  Automatic denoising of functional MRI data: Combining independent component analysis and hierarchical fusion of classifiers , 2014, NeuroImage.

[17]  Steen Moeller,et al.  Pushing spatial and temporal resolution for functional and diffusion MRI in the Human Connectome Project , 2013, NeuroImage.

[18]  Evan M. Gordon,et al.  Precision Functional Mapping of Individual Human Brains , 2017, Neuron.

[19]  Stephen M. Smith,et al.  Probabilistic independent component analysis for functional magnetic resonance imaging , 2004, IEEE Transactions on Medical Imaging.

[20]  Jean-Baptiste Poline,et al.  Brain covariance selection: better individual functional connectivity models using population prior , 2010, NIPS.

[21]  Michael C. Frank,et al.  Estimating the reproducibility of psychological science , 2015, Science.

[22]  Nora D. Volkow,et al.  Temporal Evolution of Brain Functional Connectivity Metrics: Could 7 Min of Rest be Enough? , 2016, Cerebral cortex.

[23]  Steen Moeller,et al.  Multiband multislice GE‐EPI at 7 tesla, with 16‐fold acceleration using partial parallel imaging with application to high spatial and temporal whole‐brain fMRI , 2010, Magnetic resonance in medicine.

[24]  Ying Guo,et al.  Predicting individual brain functional connectivity using a Bayesian hierarchical model , 2017, NeuroImage.

[25]  Archana Venkataraman,et al.  Intrinsic functional connectivity as a tool for human connectomics: theory, properties, and optimization. , 2010, Journal of neurophysiology.

[26]  Stephen M. Smith,et al.  Multiplexed Echo Planar Imaging for Sub-Second Whole Brain FMRI and Fast Diffusion Imaging , 2010, PloS one.

[27]  Mary E. Meyerand,et al.  The effect of scan length on the reliability of resting-state fMRI connectivity estimates , 2013, NeuroImage.

[28]  B. Efron,et al.  Data Analysis Using Stein's Estimator and its Generalizations , 1975 .

[29]  Madeleine B. Thompson,et al.  Slice Sampling with Multivariate Steps , 2011 .

[30]  Kevin Murphy,et al.  How long to scan? The relationship between fMRI temporal signal to noise ratio and necessary scan duration , 2007, NeuroImage.

[31]  Gaël Varoquaux,et al.  Population-Shrinkage of Covariance to Estimate Better Brain Functional Connectivity , 2017, MICCAI.

[32]  Yufeng Zang,et al.  Toward reliable characterization of functional homogeneity in the human brain: Preprocessing, scan duration, imaging resolution and computational space , 2013, NeuroImage.

[33]  Per Christian Hansen,et al.  Analysis of Discrete Ill-Posed Problems by Means of the L-Curve , 1992, SIAM Rev..

[34]  J. Polimeni,et al.  Blipped‐controlled aliasing in parallel imaging for simultaneous multislice echo planar imaging with reduced g‐factor penalty , 2012, Magnetic resonance in medicine.

[35]  Yufeng Zang,et al.  Functional brain hubs and their test–retest reliability: A multiband resting-state functional MRI study , 2013, NeuroImage.

[36]  M. B. Nebel,et al.  Quantifying the reliability of image replication studies: The image intraclass correlation coefficient (I2C2) , 2013, Cognitive, affective & behavioral neuroscience.

[37]  James M. Flegal,et al.  A Practical Sequential Stopping Rule for High-Dimensional Markov Chain Monte Carlo , 2016 .

[38]  Amir Shmuel,et al.  Evaluation and calibration of functional network modeling methods based on known anatomical connections , 2012, NeuroImage.

[39]  Mark W. Woolrich,et al.  Network modelling methods for FMRI , 2011, NeuroImage.

[40]  Thomas E. Nichols,et al.  A positive-negative mode of population covariation links brain connectivity, demographics and behavior , 2015, Nature Neuroscience.

[41]  Yanmei Tie,et al.  Defining language networks from resting‐state fMRI for surgical planning—a feasibility study , 2014, Human brain mapping.

[42]  N. Filippini,et al.  Group comparison of resting-state FMRI data using multi-subject ICA and dual regression , 2009, NeuroImage.

[43]  Dimitris Samaras,et al.  Deriving reproducible biomarkers from multi-site resting-state data: An Autism-based example , 2016, NeuroImage.

[44]  Gene H. Golub,et al.  Generalized cross-validation as a method for choosing a good ridge parameter , 1979, Milestones in Matrix Computation.

[45]  Brian Caffo,et al.  Comparing test-retest reliability of dynamic functional connectivity methods , 2017, NeuroImage.

[46]  C. Morris Parametric Empirical Bayes Inference: Theory and Applications , 1983 .

[47]  V. A. Morozov Regular Methods for Solving Linear and Nonlinear Ill-Posed Problems , 1984 .

[48]  Raag D. Airan,et al.  Factors affecting characterization and localization of interindividual differences in functional connectivity using MRI , 2016, Human brain mapping.

[49]  D. Yurgelun-Todd,et al.  Reproducibility of Single-Subject Functional Connectivity Measurements , 2011, American Journal of Neuroradiology.

[50]  Steen Moeller,et al.  ICA-based artefact removal and accelerated fMRI acquisition for improved resting state network imaging , 2014, NeuroImage.

[51]  Wei Sun,et al.  Partial correlation matrix estimation using ridge penalty followed by thresholding and re‐estimation , 2014, Biometrics.

[52]  Essa Yacoub,et al.  The WU-Minn Human Connectome Project: An overview , 2013, NeuroImage.

[53]  Yikai Wang,et al.  An Efficient and Reliable Statistical Method for Estimating Functional Connectivity in Large Scale Brain Networks Using Partial Correlation , 2016, Front. Neurosci..

[54]  Martin A. Lindquist,et al.  Improving reliability of subject-level resting-state fMRI parcellation with shrinkage estimators , 2014, NeuroImage.

[55]  Marina Vannucci,et al.  A Bayesian Approach for Estimating Dynamic Functional Network Connectivity in fMRI Data , 2018, Journal of the American Statistical Association.

[56]  Brian A. Nosek,et al.  Power failure: why small sample size undermines the reliability of neuroscience , 2013, Nature Reviews Neuroscience.

[57]  Theo G. M. van Erp,et al.  Multisite reliability of MR-based functional connectivity , 2017, NeuroImage.

[58]  Susan Spear Bassett,et al.  Modified test statistics by inter-voxel variance shrinkage with an application to f MRI. , 2009, Biostatistics.