Dataset Information

ModelArray: An R package for statistical analysis of fixel-wise data.

ABSTRACT: 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-wise data exist, currently available tools support only a limited number of statistical models. Here we introduce ModelArray, an R package for mass-univariate statistical analysis of fixel-wise data. At present, ModelArray supports linear models as well as generalized additive models (GAMs), which are particularly useful for studying nonlinear effects in lifespan data. In addition, ModelArray also aims for scalable analysis. With only several lines of code, even large fixel-wise datasets can be analyzed using a standard personal computer. Detailed memory profiling revealed that ModelArray required only limited memory even for large datasets. As an example, we applied ModelArray to fixel-wise data derived from diffusion images acquired as part of the Philadelphia Neurodevelopmental Cohort (n = 938). ModelArray 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 a flexible and efficient platform for statistical analysis of fixel-wise data.

SUBMITTER: Zhao C

PROVIDER: S-EPMC10119782 | biostudies-literature | 2023 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

ModelArray: An R package for statistical analysis of fixel-wise data.

Zhao Chenying C Tapera Tinashe M TM Bagautdinova Joëlle J Bourque Josiane J Covitz Sydney S Gur Raquel E RE Gur Ruben C RC Larsen Bart B Mehta Kahini K Meisler Steven L SL Murtha Kristin K Muschelli John J Roalf David R DR Sydnor Valerie J VJ Valcarcel Alessandra M AM Shinohara Russell T RT Cieslak Matthew M Satterthwaite Theodore D TD

NeuroImage 20230315

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-wise data exist, currently available tools support only a limited number of statistical models. Here we introduce ModelArray, an R package for mass-univariate statistical analysis of fixel-wise data. At present, ModelArray supports l ...[more]

PMID: 36931330

Similar Datasets

Project description:In healthy individuals, behavioral outcomes are highly associated with the variability on brain regional structure or neurochemical phenotypes. Similarly, in the context of neurodegenerative conditions, neuroimaging reveals that cognitive decline is linked to the magnitude of atrophy, neurochemical declines, or concentrations of abnormal protein aggregates across brain regions. However, modeling the effects of multiple regional abnormalities as determinants of cognitive decline at the voxel level remains largely unexplored by multimodal imaging research, given the high computational cost of estimating regression models for every single voxel from various imaging modalities. VoxelStats is a voxel-wise computational framework to overcome these computational limitations and to perform statistical operations on multiple scalar variables and imaging modalities at the voxel level. VoxelStats package has been developed in Matlab(®) and supports imaging formats such as Nifti-1, ANALYZE, and MINC v2. Prebuilt functions in VoxelStats enable the user to perform voxel-wise general and generalized linear models and mixed effect models with multiple volumetric covariates. Importantly, VoxelStats can recognize scalar values or image volumes as response variables and can accommodate volumetric statistical covariates as well as their interaction effects with other variables. Furthermore, this package includes built-in functionality to perform voxel-wise receiver operating characteristic analysis and paired and unpaired group contrast analysis. Validation of VoxelStats was conducted by comparing the linear regression functionality with existing toolboxes such as glim_image and RMINC. The validation results were identical to existing methods and the additional functionality was demonstrated by generating feature case assessments (t-statistics, odds ratio, and true positive rate maps). In summary, VoxelStats expands the current methods for multimodal imaging analysis by allowing the estimation of advanced regional association metrics at the voxel level.

Project description:ObjectivesMyocardial injury assessment from delayed enhancement magnetic resonance images is routinely limited to global descriptors such as size and transmurality. Statistical tools from computational anatomy can drastically improve this characterization, and refine the assessment of therapeutic procedures aiming at infarct size reduction. Based on these techniques, we propose a new characterization of myocardial injury up to the pixel resolution. We demonstrate it on the imaging data from the Minimalist Immediate Mechanical Intervention randomized clinical trial (MIMI: NCT01360242), which aimed at comparing immediate and delayed stenting in acute ST-Elevation Myocardial Infarction (STEMI) patients.MethodsWe analyzed 123 patients from the MIMI trial (62 ± 12 years, 98 male, 65 immediate 58 delayed stenting). Early and late enhancement images were transported onto a common geometry using techniques inspired by statistical atlases, allowing pixel-wise comparisons across population subgroups. A practical visualization of lesion patterns against specific clinical and therapeutic characteristics was also proposed using state-of-the-art dimensionality reduction.ResultsInfarct patterns were roughly comparable between the two treatments across the whole myocardium. Subtle but significant local differences were observed for the LCX and RCA territories with higher transmurality for delayed stenting at lateral and inferior/inferoseptal locations, respectively (15% and 23% of myocardial locations with a p-value <0.05, mainly in these regions). In contrast, global measurements were comparable for all territories (no statistically significant differences for all-except-one measurements before standardization / for all after standardization), although immediate stenting resulted in more subjects without reperfusion injury.ConclusionOur approach substantially empowers the analysis of lesion patterns with standardized comparisons up to the pixel resolution, and may reveal subtle differences not accessible with global observations. On the MIMI trial data as illustrative case, it confirmed its general conclusions regarding the lack of benefit of delayed stenting, but revealed subgroups differences thanks to the standardized and finer analysis scale.

Dataset Information

ModelArray: An R package for statistical analysis of fixel-wise data.

Publications

ModelArray: An R package for statistical analysis of fixel-wise data.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets