1 Department of Physics, Concordia University, Montreal, Canada.
2 School of Health, Concordia University, Montreal, Canada.
3 EPIC Centre, Montreal Heart Institute, Montreal, Canada.
4 Department of Psychology, Concordia University, Montreal, Canada.
5 Department of Electrical Engineering, Concordia University, Montreal, Canada.
6 ViTAA medical solutions, Montreal, Canada.
7 Biospective Inc., Montreal, Canada.
8 Neurology and Neurosurgery Department, Montreal Neurological Institute, McGill University, Montreal, Canada.
9 McConnell Brain Imaging Center, Montreal Neurological Institute, Montreal, Canada.
10 Ludmer Center for NeuroInformatics and Mental Health, Montreal, Canada.
11 Department of Neurology, Max Planck Institute for Human Cognitive and Brain Sciences, Leipzig, Germany.

Multivariate approaches have recently gained in popularity to address the physiological unspecificity of neuroimaging metrics and to better characterize the complexity of biological processes underlying behavior. However, commonly used approaches are biased by the intrinsic associations between variables, or they are computationally expensive and may be more complicated to implement than standard univariate approaches. Here, we propose using the Mahalanobis distance (D2), an individual-level measure of deviation relative to a reference distribution that accounts for covariance between metrics. To facilitate its use, we introduce an open-source python-based tool for computing D2 relative to a reference group or within a single individual: the MultiVariate Comparison (MVComp) toolbox. The toolbox allows different levels of analysis (i.e., group- or subject-level), resolutions (e.g., voxel-wise, ROI-wise) and dimensions considered (e.g., combining MRI metrics or WM tracts). Several example cases are presented to showcase the wide range of possible applications of MVComp and to demonstrate the functionality of the toolbox. The D2 framework was applied to the assessment of white matter (WM) microstructure at 1) the group-level, where D2 can be computed between a subject and a reference group to yield an individualized measure of deviation. We observed that clustering applied to D2 in the corpus callosum yields parcellations that highly resemble known topography based on neuroanatomy, suggesting that D2 provides an integrative index that meaningfully reflects the underlying microstructure. 2) At the subject level, D2 was computed between voxels to obtain a measure of (dis)similarity. The loadings of each MRI metric (i.e., its relative contribution to D2) were then extracted in voxels of interest to showcase a useful option of the MVComp toolbox. These relative contributions can provide important insights into the physiological underpinnings of differences observed. Integrative multivariate models are crucial to expand our understanding of the complex brain-behavior relationships and the multiple factors underlying disease development and progression. Our toolbox facilitates the implementation of a useful multivariate method, making it more widely accessible.