Statistical Instability of TBSS Analysis Based on DTI Fitting Algorithm

Ivan I Maximov,Heike Thönneßen,Kerstin Konrad,Laura Amort,Irene Neuner,N Jon Shah,Ivan I. Maximov,N. Jon Shah
DOI: https://doi.org/10.1111/jon.12215
2015-02-14
Journal of Neuroimaging
Abstract:Voxel-based DTI analysis is an important approach in the comparison of subject groups by detecting and localizing gray and white matter changes in the brain. One of the principal problems for intersubject comparison is the absence of a "gold standard" processing pipeline. As a result, contradictory results may be obtained from identical data using different data processing pipelines, for example, in the data normalization or smoothing procedures. Tract-based spatial statistics (TBSS) shows potential to overcome this problem by automatic detection of white matter changes and decreasing variation in the performed analysis. However, skeleton projection approaches, such as TBSS, critically depend on the accuracy of the diffusion scalar metric estimations. In this work, we demonstrate that the agreement and reliability of TBSS results depend on the applied DTI data processing algorithm. Statistical tests have been performed using two in vivo measured datasets and compared with different implementations of the least squares algorithm. As a result, we recommend repeating TBSS analysis using different fitting algorithms, in particular, using on iteratively-assessed robust estimators, as accurate and more reliable approach in voxel-based analysis, particularly, for TBSS. Repeating TBSS analysis allows one to detect and localize suspicious regions in white matter which were estimated as the regions with significant difference. Finally, we did not find a favorite fitting algorithm (or class of them) which can be marked as more reliable for group comparison.
radiology, nuclear medicine & medical imaging,clinical neurology,neuroimaging
What problem does this paper attempt to address?
### What problem does this paper attempt to solve? This paper aims to explore whether the results of Tract - based Spatial Statistics (TBSS) analysis based on Diffusion Tensor Imaging (DTI) depend on the DTI data processing algorithms used. Specifically, the author hopes to reveal the statistical instability and reproducibility problems of TBSS analysis by studying the influence of different fitting algorithms on TBSS results. #### Background and Motivation 1. **Introduction to DTI and TBSS**: - Diffusion Tensor Imaging (DTI) is a powerful technique for in - vivo detection of the white matter structure of the human brain. - TBSS is a commonly used method that reduces the variation in inter - group comparisons by projecting the FA images of all subjects onto a "skeleton" and automatically detects white matter changes. 2. **Existing Problems**: - Different research teams use different DTI data processing pipelines, which may lead to contradictory results even from the same data. - The results of TBSS analysis are highly dependent on the accuracy of diffusion scalar metric estimations, and these estimations depend on the fitting algorithms used. #### Research Objectives - **Verify the Dependence of TBSS Results**: Study the influence of different DTI fitting algorithms (such as ordinary least squares, weighted least squares, nonlinear least squares, RESTORE, etc.) on the results of TBSS analysis. - **Improve the Reliability of TBSS**: It is recommended to use multiple fitting algorithms in TBSS analysis, especially robust estimators with iterative evaluation, to improve the accuracy and reliability of the results. - **Identify Suspicious Areas**: Detect and locate white matter suspicious areas due to statistical instability by repeating TBSS analysis. #### Main Findings - **Influence of Different Algorithms**: There are significant differences in FA/MD values produced by different fitting algorithms, especially between similar methods implemented in different software packages. - **Variability of TBSS Results**: For the same data set, TBSS results obtained by different fitting algorithms may be completely different, and even show differences within the same group. - **Advantages of Robust Estimators**: The robust estimator based on the modified least trimmed squares (MLTS) shows better consistency and stability in both data sets. #### Conclusion The author recommends that when performing TBSS analysis, multiple fitting algorithms, especially robust estimators, should be considered to ensure the reliability and reproducibility of the results. In addition, they emphasize the importance of choosing appropriate fitting algorithms in multi - center clinical studies to reduce the variability of results and improve the credibility of the study. ### Formula Explanation The formulas involved in this paper are mainly related to Diffusion Tensor Imaging (DTI), such as the calculation of diffusion tensors and the estimation of scalar metrics (such as FA, MD). The following is the Markdown - format display of several key formulas: 1. **Definition of Diffusion Tensor**: \[ D=\begin{pmatrix} D_{xx}&D_{xy}&D_{xz}\\ D_{yx}&D_{yy}&D_{yz}\\ D_{zx}&D_{zy}&D_{zz} \end{pmatrix} \] 2. **Fractional Anisotropy (FA)**: \[ FA = \sqrt{\frac{3}{2}}\cdot\frac{\sqrt{(\lambda_1-\bar{\lambda})^2+(\lambda_2 - \bar{\lambda})^2+(\lambda_3-\bar{\lambda})^2}}{\sqrt{\lambda_1^2+\lambda_2^2+\lambda_3^2}} \] where \(\lambda_1,\lambda_2,\lambda_3\) are the eigenvalues of the diffusion tensor, and \(\bar{\lambda}\) is their average. 3. **Mean Diffusivity (MD)**: \[ MD=\frac{\lambda_1+