QDECR: A Flexible, Extensible Vertex-Wise Analysis Framework in R

The cerebral cortex is fundamental to the functioning of the mind and body. In vivo cortical morphology can be studied through magnetic resonance imaging in several ways, including reconstructing surface-based models of the cortex. However, existing software for surface-based statistical analyses cannot accommodate “big data” or commonly used statistical methods such as the imputation of missing data, extensive bias correction, and non-linear modeling. To address these shortcomings, we developed the QDECR package, a flexible and extensible R package for group-level statistical analysis of cortical morphology. QDECR was written with large population-based epidemiological studies in mind and was designed to fully utilize the extensive modeling options in R. QDECR currently supports vertex-wise linear regression. Design matrix generation can be done through simple, familiar R formula specification, and includes user-friendly extensions for R options such as polynomials, splines, interactions and other terms. QDECR can handle unimputed and imputed datasets with thousands of participants. QDECR has a modular design, and new statistical models can be implemented which utilize several aspects from other generic modules which comprise QDECR. In summary, QDECR provides a framework for vertex-wise surface-based analyses that enables flexible statistical modeling and features commonly used in population-based and clinical studies, which have until now been largely absent from neuroimaging research.

[1]  Gary King,et al.  Amelia II: A Program for Missing Data , 2011 .

[2]  Tomáš Paus,et al.  Population neuroscience: Why and how , 2010, Human brain mapping.

[3]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[4]  A. Zeileis,et al.  Extended Model Formulas in R : Multiple Parts and Multiple Responses , 2010 .

[5]  KJ Worsley,et al.  SurfStat: A Matlab toolbox for the statistical analysis of univariate and multivariate surface and volumetric data using linear mixed effects models and random field theory , 2009, NeuroImage.

[6]  Brian Avants,et al.  Neuroconductor: an R platform for medical imaging analysis , 2019, Biostatistics.

[7]  Bruce Fischl,et al.  FreeSurfer , 2012, NeuroImage.

[8]  Stef van Buuren,et al.  MICE: Multivariate Imputation by Chained Equations in R , 2011 .

[9]  Karl J. Friston,et al.  A Voxel-Based Method for the Statistical Analysis of Gray and White Matter Density Applied to Schizophrenia , 1995, NeuroImage.

[10]  A. Dale,et al.  Improved Localizadon of Cortical Activity by Combining EEG and MEG with MRI Cortical Surface Reconstruction: A Linear Approach , 1993, Journal of Cognitive Neuroscience.

[11]  Alan C. Evans,et al.  Amygdalar reactivity is associated with prefrontal cortical thickness in a large population-based sample of adolescents , 2019, PloS one.

[12]  Thomas E. Nichols,et al.  Statistical Challenges in “Big Data” Human Neuroimaging , 2018, Neuron.

[13]  A. Gelman,et al.  Multiple Imputation with Diagnostics (mi) in R: Opening Windows into the Black Box , 2011 .

[14]  Anders M. Dale,et al.  The Adolescent Brain Cognitive Development (ABCD) study: Imaging acquisition across 21 sites , 2018, Developmental Cognitive Neuroscience.

[15]  Lena Osterhagen,et al.  Multiple Imputation For Nonresponse In Surveys , 2016 .

[16]  Henry Markram,et al.  A Brief History of Simulation Neuroscience , 2019, Front. Neuroinform..

[17]  Ludovica Griffanti,et al.  Image processing and Quality Control for the first 10,000 brain imaging datasets from UK Biobank , 2017, NeuroImage.

[18]  D. Rubin INFERENCE AND MISSING DATA , 1975 .

[19]  Ayse Pinar Saygin,et al.  Smoothing and cluster thresholding for cortical surface-based group analysis of fMRI data , 2006, NeuroImage.

[20]  Raphael Gottardo,et al.  Orchestrating high-throughput genomic analysis with Bioconductor , 2015, Nature Methods.

[21]  Peter Bühlmann,et al.  MissForest - non-parametric missing value imputation for mixed-type data , 2011, Bioinform..

[22]  Bruce Fischl,et al.  False positive rates in surface-based anatomical analysis , 2018, NeuroImage.