Independent component analysis (ICA) on group-level voxel-based morphometry (VBM) produces the coefficient matrix and the component matrix. The former contains variability among multiple subjects for further statistical analysis, and the latter reveals spatial maps common for all subjects. ICA algorithms converge to local optimization points in practice and the mostly applied stability investigation approach examines the stability of the extracted components. We found that the practically stable components do not guarantee to produce the practically stable coefficients of ICA decomposition for the further statistical analysis. Consequently, we proposed a novel approach including two steps: 1), the stability index for the coefficient matrix and the stability index for the component matrix were examined, respectively; 2) the two indices were multiplied to analyze the stability of ICA decomposition. The proposed approach was used to study the sMRI data of Type II diabetes mellitus group (DM) and the healthy control group (HC). Group differences in VBM were found in the superior temporal gyrus. Besides, it was revealed that the VBMs of the region of the HC group were significantly correlated with Montreal Cognitive Assessment (MoCA) describing the level of cognitive disorder. In contrast to the widely applied approach to investigating the stability of the extracted components for ICA decomposition, we proposed to examine the stability of ICA decomposition by fusion the stability of both coefficient matrix and the component matrix. Therefore, the proposed approach can examine the stability of ICA decomposition sufficiently.