Subjects
The structural sMRI data used in our study were acquired from the Soochow University, which is composed of 68 undergraduates. The study was approved by the Ethics Committee of the Third Affiliated Hospital of Soochow University. Written informed consents was obtained from all subjects. All subjects did not receive stimulants or hypnotics before acquisition in order to keep them awake and let the brain work normally. All participants' vision was normal or corrected to normal, and they were right-handed. After the test, each participant will receive a small gift or financial reward. All subjects are required to perform Rosenberg Self-esteem Scale (RSES) test. The RSES is originally developed by Rosenberg in 1965 to assess the overall feelings of undergraduates about self-worth and self-acceptance. It is the most used self-esteem measurement tool in the psychology community [18]. We ranked the RSES test scores from highest to lowest, and then divided them into two groups: high self-esteem group and low self-esteem group. Table 5 provides detailed information of all participants.
Imaging acquisition and preprocessing
All images were collected on a 3T Siemens Medical Systems equipment. The acquisition parameters are set as: echo time (TE) = 2.98 ms, repetition time (TR) = 2300 ms, flip angle (FA) = 9 deg, voxel size = 1 × 1 × 1 mm3, slice thickness = 1 mm, field of view (FoV) = 256mm.
We use an automatic pipeline for sMRI image processing. Firstly, we adjusted the image orientation (axial, coronal, and sagittal) to match the template image, and performed offset field correction to remove the gray-scale unevenness of the image [19]. Secondly, the brain image was extracted by removing the skull and cerebellum [20]. Thirdly, gray matter (GM), white matter (WM) and cerebrospinal fluid (CSF) were segmented from the background [21]. Fourth, the segmented image was registered to the template labeled with the Automated Anatomical Labeling (AAL) template [22]. Fifth, in order to calculate the morphological features based on the cortex, the middle layer of the cerebral cortex was depicted [23]. After the whole processing, the morphological measurements of GM volume, WM volume, CSF volume, cortical thickness, and cortical surface area of each ROI were obtained for each subject. It should be noted that we removed 12 subcortical ROIs from AAL template considering that the cerebral cortex contains more neurons.
Classification framework
The framework of the proposed classification algorithm based on multi-resolution ROI brain network is shown in Figure 5, mainly including multiple anatomical network construction, feature selection, and classification. Multi-resolution ROI based multiple anatomical brain network were constructed based on morphological features (volume of different brain tissue, cortical thickness, and cortical surface area). Feature selection can reduce the dimensionality of high-dimensional brain network features, only retaining the features that can maximize the specificity of the subjects. The optimal feature subset can be trained by the classifier as neuroimaging markers representing different self-esteem levels.
Construction of multiple anatomical networks
Feature selection
In order to reduce the feature dimension and filter out the most discriminative features, we adopted several feature selection methods. First, we preliminarily select the features by comparing the statistics of different features. The statistical t-test (p <0.05) is adopted to remove features with small differences (the features with small differences are difficult to distinguish the two groups). Then, another filter-based feature selection method called minimum redundancy and maximum correlation (mRMR) is used to remove the redundant features [24]. The core idea of mRMR is to maximize the correlation between features and classification variables, and minimize the correlation between different features. After the above two filter-based feature selections, the machine learning recursive feature elimination (SVM-RFE) method [25] is used to further reduce the feature dimension. SVM-RFE is proposed in classification of cancer, and has good performance and strong generalization ability. It is the combination of SVM and subsequent search strategy. It trains samples through the model, and then ranks the scores of each feature to remove the feature with the smallest score, and then trains the model again with the remaining features for the next iteration, and finally selects the number of features that are needed. After completing the entire feature selection steps, the optimal feature subset is obtained.
Classification using multi-kernel SVM
There are various types of features in the multiple brain network, one is the high-resolution ROI features in the fourth layer, and the other is the brain network features corresponding to different layers. Multi-kernel machine learning method can integrate these various types of features into a final classifier. Firstly, a Gaussian Radial Basis Function (RBF) kernel function is used to construct a kernel matrix for each type of feature. Secondly, the two kernel matrices are integrated into the multi-kernel matrix through appropriate weight coefficients [25]. Comparing the results of using linear kernel function and using RBF function (non-linear), we discover that the RBF kernel can significantly improve the classification performance. Therefore, we choose the RBF kernel function to construct the multi-kernel classifier. Finally, the optimal features subset can be obtained.
Cross-validation
The nested cross-validation method has been applied in our previous research [26]. In the inner loop, the training set are used to determine the parameters of the classifier. In the outer loop, the testing set is used to evaluate the generalization ability of the classifier. It should be noted that at the beginning of the experiment, the entire data set was randomly divided into two parts, one for training and the other one for testing. The training set and testing set can be exchanged throughout the verification process, while the processing steps remain unchanged.