Statistical analyses of motion-corrupted MRI relaxometry data computed from multiple scans

https://doi.org/10.1016/j.jneumeth.2023.109950Get rights and content
Under a Creative Commons license
open access

Highlights

  • Motion violates the assumption of uniform noise level in statistical analyses.

  • QUIQI restores the validity and maximizes the sensitivity of statistical analyses.

  • We extend QUIQI to quantitative MRI maps computed from multiple scans.

  • We provide a roadmap to optimize the method’s performance in future applications.

  • QUIQI is available in the hMRI toolbox https://hmri-group.github.io/hMRI-toolbox/.

Abstract

Background

Consistent noise variance across data points (i.e. homoscedasticity) is required to ensure the validity of statistical analyses of MRI data conducted using linear regression methods. However, head motion leads to degradation of image quality, introducing noise heteroscedasticity into ordinary-least square analyses.

New method

The recently introduced QUIQI method restores noise homoscedasticity by means of weighted least square analyses in which the weights, specific for each dataset of an analysis, are computed from an index of motion-induced image quality degradation. QUIQI was first demonstrated in the context of brain maps of the MRI parameter R2 * , which were computed from a single set of images with variable echo time. Here, we extend this framework to quantitative maps of the MRI parameters R1, R2 * , and MTsat, computed from multiple sets of images.

Results

QUIQI restores homoscedasticity in analyses of quantitative MRI data computed from multiple scans. QUIQI allows for optimization of the noise model by using metrics quantifying heteroscedasticity and free energy.

Comparison with existing methods

QUIQI restores homoscedasticity more effectively than insertion of an image quality index in the analysis design and yields higher sensitivity than simply removing the datasets most corrupted by head motion from the analysis.

Conclusion

QUIQI provides an optimal approach to group-wise analyses of a range of quantitative MRI parameter maps that is robust to inherent homoscedasticity.

Keywords

Heteroscedasticity
Statistical group analysis
Quantitative MRI
Motion corruption

Data Availability

Material related to this paper is available here (Raynaud et al., 2023): https://zenodo.org/record/7692074#. ZACk9HaZOUk. This material includes maps of R2 *(2), R2 *(3), R1 and MTsat from 123 participants, and the results of the analyses presented in this manuscript conducted on this data.

Cited by (0)