ModelArray: a memory-efficient R package for statistical analysis of fixel data

Diffusion MRI is the dominant non-invasive imaging method used to characterize white matter organization in health and disease. Increasingly, fiber-specific properties within a voxel are analyzed using fixels. While tools for conducting statistical analyses of fixel data exist, currently available tools are memory intensive, difficult to scale to large datasets, and support only a limited number of statistical models. Here we introduce ModelArray, a memory-efficient R package for mass-univariate statistical analysis of fixel data. With only several lines of code, even large fixel datasets can be analyzed using a standard personal computer. At present, ModelArray supports linear models as well as generalized additive models (GAMs), which are particularly useful for studying nonlinear effects in lifespan data. Detailed memory profiling revealed that ModelArray required only limited memory even for large datasets. As an example, we applied ModelArray to fixel data derived from diffusion images acquired as part of the Philadelphia Neurodevelopmental Cohort (n=938). ModelArray required far less memory than existing tools and revealed anticipated nonlinear developmental effects in white matter. Moving forward, ModelArray is supported by an open-source software development model that can incorporate additional statistical models and other imaging data types. Taken together, ModelArray provides an efficient and flexible platform for statistical analysis of fixel data. HIGHLIGHTSO_LIModelArray is an R package for mass-univariate statistical analysis of fixel data C_LIO_LIModelArray is memory-efficient even for large-scale datasets C_LIO_LIModelArray supports linear and nonlinear modeling and is extensible to more models C_LIO_LIModelArray facilitates easy statistical analysis of large-scale fixel data C_LI Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=80 SRC="FIGDIR/small/499631v1_ufig1.gif" ALT="Figure 1"> View larger version (17K): org.highwire.dtl.DTLVardef@1946469org.highwire.dtl.DTLVardef@14c3076org.highwire.dtl.DTLVardef@10197beorg.highwire.dtl.DTLVardef@19170cd_HPS_FORMAT_FIGEXP M_FIG C_FIG