# Improved Multiclassification of Schizophrenia Based on Xgboost and Information Fusion for Small Datasets

Authors: Wenjing Zhu; Shoufeng Shen; Zhijun Zhang
Journal: Computational and Mathematical Methods in Medicine
 (2022)

---

## Abstract

To improve the performance in multiclass classification for small datasets, a new approach for schizophrenic classification is proposed in the present study. Firstly, the Xgboost classifier is introduced to discriminate the two subtypes of schizophrenia from health controls by analyzing the functional magnetic resonance imaging (fMRI) data, while the gray matter volume (GMV) and amplitude of low-frequency fluctuations (ALFF) are extracted as the features of classifiers. Then, the D-S combination rule of evidence is used to achieve fusion to determine the basic probability assignment based on the output of different classifiers. Finally, the algorithm is applied to classify 38 healthy controls, 16 deficit schizophrenic patients, and 31 nondeficit schizophrenic patients. 10-folds cross-validation method is used to assess classification performance. The results show the proposed algorithm with a sensitivity of 73.89%, which is higher than other classification algorithms, such as supported vector machine (SVM), logistic regression (LR),K-nearest neighbor (KNN) algorithm, random forest (RF), BP neural network (NN), classification and regression tree (CART), naive Bayes classifier (NB), extreme gradient boosting (Xgboost), and deep neural network (DNN). The accuracy of the fusion algorithm is higher than that of classifier based on the GMV or ALFF in the small datasets. The accuracy rate of the improved multiclassification method based on Xgboost and fusion algorithm is higher than that of other machine learning methods, which can further assist the diagnosis of clinical schizophrenia.


---

## Body

## 1. Introduction

Schizophrenia (SZ) is a serious mental illness that interferes with a person’s ability to think clearly, manage emotions, make decisions, and relate to others [1]. The positive symptoms, such as hallucinations and delusions, can lead to suicidal or aggressive behavior, while negative symptoms and cognitive impairment lead to a decline in quality of life and social function; all these symptoms will cause tremendous human suffering and economic burden [2]. However, SZ is diagnosed on the basis of clinical evaluation of symptoms and functional, no objective diagnostic biomarker set. In addition, there are two main types of schizophrenia, which are called deficit schizophrenia (DS) and nondeficit schizophrenia (NDS) [3]. Deficit syndrome of schizophrenia, also called negative symptoms of schizophrenia, includes social withdrawal, loss of motivation, poverty of speech, and blunting of affect. Compared with NDS, DS has greater cognitive impairment, worse long-term prognosis, and lower recovery rates, which persist or are found even during psychotic remissions [4, 5]. Therefore, it is quite important to diagnose SZ accurately and discriminate the two subtypes of SZ from healthy control (HC), particularly the discrimination between the DS and NDS.Classification is a machine learning algorithm where we get the labeled data as input and we need to predict the output into a class [6]. If there are two classes, then it is called binary classification [7]. If there are more than two classes, then it is called multiclass classification. Nowadays, classification algorithm is widely used in the medical diagnosis, especially in the field of mental disorders. Early published research applied support vector machines with radial basis function kernel method to classify 15 schizophrenic patients and 15 HCs based on the structural image of the hippocampal complex only with 63% classification accuracy [8]. With the rapid development of computational psychiatry, a growing body of classification approaches are applied to discriminate SZ in recent years, such as logistic regression (LR), support vector machine (SVM), neural network (NN), random forest (RF), extreme gradient boosting (XGBoost), and deep learning [9]. In the early studies, many researchers focused on the binary classification problem in SZ. Greenstein et al. proposed the logistic regress classifier to discriminate 99 SZ patients from 99 HCs with 73.7% accuracy [10]. Nieuwenhuis et al. proposed the SVM classifier to discriminate 128 SZ patients from 111 HCs and achieved an accuracy of 71.4% [11]. Thereby, in the past years, feature reduction approaches were discussed and applied to improve the performance of classification. Ershad and Hashemi proposed the dispelling reduction approaches [12], Juneja and his colleagues obtained the discriminative features by using SVD model and a novel multivariate feature selection algorithm [13]. However, the accuracy was not high enough based on the classical classification algorithms, which was usually less than 75%. Then, many new classifiers and feature selection approaches are proposed to improve the classification performance. Up to 2018, Wang et al. developed the SVM model to discriminate SZ from HCs and achieved an accuracy of 92.4% [14]. In 2020, Kim et al. proposed the feature reduction method when there are redundant or correlated features based on the FDR value and achieved an accuracy of 96.2% [15]. Patel et al. proposed a classification algorithm to discriminate SZ versus HCs by busing deep learning in fMRI, and the accuracy was 92% [16]. Nowadays, the binary classification is not a hard work in the field of mental disorders. However, there is not only one type of SZ, such as DS and NDS, which is more difficult to discriminate each other or from HCs. In order to solve this problem, multiple classification methods for schizophrenic subtypes are necessary.Few previous literatures have reported multiclass classifications for different types of psychiatric patients, and most of these classifications used traditional machine learning methods such as SVM and LR. For examples, Zhu et al. proposed the SVM model to classify first-episode, drug-naive SZ, ultrahigh risk for psychosis and HC with the global balanced accuracy only 73.37% and the sensitivity only 68.42%, using the fivefold cross-validation method [17]. Soon afterwards, multiple classification methods have been explored, such as three SVM models to classify the SZ, bipolar disorder, and HC [18] and classify depression, bipolar disorder, and HC [19]; an SVM combined with recursive feature elimination was used to classify first-episode SZ, chronic SZ, and HC [20]. Unfortunately, almost all the accuracy rate of the multiclass classification is less than 70%. In addition, because of the poor coordination in psychiatric patients, the amount of imaging data is generally small, which leads to an accuracy far below 70%.In recent years, deep learning is widely used in the pattern recognition. Zeng et al. proposed a deep discriminant autoencoder network to learn imaging site-shared functional connectivity features to discriminate SZ from normal subjects. In their work, the accuracy of 85% is achieved [21]. Oh et al. collected 873 structural MRI datasets and discriminate the SZ from normal subjects by using a deep convolutional neural network [22]. Srinivasagopalan et al. proposed a deep learning algorithm for diagnosing SZ [23]. These deep learning algorithms are based on the original image data, such as CNN and DNN. However, the original image data is usually quite difficult to obtain. Many deep learning algorithms on small dataset is usually overfitting. Therefore, the machine learning algorithms are more suitable to improve the performance of the classical classifiers.To the best of our knowledge, there is no study to implement the multiclass classification of DS, NDS, and HC based on multimodal imaging data of schizophrenia. Therefore, in order to achieve multiclass classification of schizophrenia and obtain higher classification results in small data, a new classification algorithm is proposed in this paper. In this algorithm, GMV and ALFF are selected as the features to construct multiclassifier based on Xgboost, respectively. Then, the fused model is built to improve the accuracy for the small datasets. The D-S fusion model is used to combine the output from different classifiers to determine the probability assignment for different subtypes and HC. The rest of this paper is structured as follows. In Section2, the Xgboost classifier is proposed to discriminate DS and NDS from HC and the fusion model is introduced to combine information of output. The results obtained by applying our model are shown in Section 3. In Section 4, the main contribution of this paper is summarized.


## 2. The Fused Classification Algorithm Based on Xgboost for Three-Class Classification

In this section, a fused classification algorithm is proposed to improve the accuracy for the small datasets. This algorithm is applied to discriminate the subtypes of SZ from HC. There are three labels (DS, NDS, and HC) that should be assigned to each collected subject. To solve this problem, an improved multiclassification algorithm is introduced. Firstly, Xgboost algorithm is applied to classify them, which is one of the most widely used machine learning algorithms in classification problems [24]. The classifiers are constructed based on the features of GMV or ALFF, which is extracted from fMRI data. Then, the fusion model is used to combine the output information of the different classifier to determine the probability assignment of each class. Finally, the test subject will be classified into the class with the maximum probability. The flow chart of the proposed algorithm is shown in Figure 1.Figure 1
Flow chat of improved multiclassification algorithm.


### 2.1. Classifier Based on the Xgboost Algorithm

The Xgboost algorithm is composed of many weak classification and regression trees (CARTs). Taking theith dataset (xi,yi) as an example, xi is the input variable with several attributes of fMRI data and yi is the real value of the given subject. For example, yi=1,0,0 means the ith subject is HC, yi=0,1,0 means the ith subject is DS, and yi=0,0,1 means the ith subject is NDS.Then, an Xgboost model can be mathematically expressed in the following form [25]:
(1)y^i=∑k=1Kfkxi,where K is the number of the CARTs, fk is the predicted value of each independent CART, and y^i is the predicted value with respect to input xi.The additive training model of Xgboost can be expressed as(2)y^i0=0,y^ik=∑t=1kftxi=y^kk−1+fkxi,where y^ik is the predicted value of the kth CART.The objective function of Xgboost includes a loss function and regularization term, which is expressed as(3)obj=∑i=1nlyi,y^i+∑k=1KΩfk,where lyi,y^i can be used to measure the error between predicted value and real value, n is the number of the subjects, and Ω is the regularization item to avoid overfitting.The specific form ofΩfk of the kth CART is given as
(4)Ωfk=γT+12λw2,where γ and λ present the penalty coefficients, T is the number of leaf nodes, and w is the weight of the leaf nodes.Then, the objective function of thetth step objt can be calculated by Equation (5) based on the previous step objt−1 based on the Equations (2) and (3).
(5)objt=∑i=1nlyi,y^it−1+ftxi+∑k=1t−1Ωfk+Ωft.By applying the second-order Taylor expansion to above equation, the objective function can be transformed into(6)objt=∑i=1nlyi,y^it−1+giftxi+12hift2xi+Ωft+Const.In the above expression,Const is a constant term at the step t; the parameters pair gi and hi can be calculated as
(7)gi=∂lyi,y^it−1∂y^it−1,hi=∂2lyi,y^it−1∂y∧it−12.Becauselyi,y^it−1is the constant item, the objective function can be rewritten as
(8)objt=∑i=1ngiftxi+12hift2xi+γT+12λ∑j=1Twj2+const′,where const′ is a new constant item at step t.According to the definition offk, fk can be written in the following form as
(9)fkx=wqx.Then, Equation (8) can be rewritten in the following form as
(10)objt=∑j=1T∑i∈Ijgiwj+12∑i∈Ijhi+λwj2+γT+const′,denoted by∑i∈Ijgi=Gjand∑i∈Ijhi=Hj; then, the objective function is expressed as
(11)objt=∑j=1TGjwj+12Hj+λwj2+γT+const′.The leaf nodes of thetth CART are each independent; Gj and Hj are the determined items. Then, minimizing the function equation (11), the optimal parameter wj can be calculated as
(12)wj∗=−GjHj+λ.Therefore, the final objective function is shown in the following form as(13)obj∗=−12∑j=1TGj2Hj+λ+γT+const′.The splitting algorithm [26] based on the above function is used to find the best split in Xgboost by
(14)Gain=12GL2HL+λ+GR2HR+λ−GL+GR2HL+HR+λ−γ.The gain function has four terms: the first two terms are the profits of left and right parts of a node, whereGL and GR are the left and right parts of Gj and HL and HR are the left and right parts of Hj; and the third item is the total profit of that node. The last item is the regularization item for preventing overfitting. The greedy algorithm determines whether a node obtains the maximum gain. Thus far, the optimal tree structure that maximizes the gain can be generated.The above description leads to the split finding algorithm for Xgboost presented as Algorithm1.Algorithm 1: Greedy algorithm for split finding.
Require:I, instance set of current node; d, feature dimensionEnsure: Split with max score1.Gain⟵0, G⟵∑i∈Ijgi, H⟵∑i∈Ijhi2.fork=1 to mdo3.GL⟵0, HL⟵04.for alli in sorted I by xjkdo5.GL⟵GL+gj, HL⟵HL+hj6.GR⟵G−GL, HR⟵H−HL7.score⟵maxscore,GL2/HL+λ+GR2/HR+λ−G2/H+λ8.end for9. end for


### 2.2. The Fusion Model Based on the D-S Evidence Theory

Analyses of the amplitude of low-frequency fluctuations (ALFF) and gray matter volume (GMV) are two important methods used in fMRI studies. Selecting GMV as the feature to construct the classifier based on the Xgboost algorithm, the predicted value of theith subject can be obtained and denoted as y^i1, while the other predicted value y^i2can be calculated by selecting ALFF as the feature,y^i1≠y^i2usually. Taking the ith subject as an example, the output information of the classifiers is y^i1=0.7,0.3,0.9 and y^i2=0.7,0.5,0.6. According to the output of classifier based on the feature of GMV, this subject should be classified into NDS, while the classifier based on the feature of ALFF will classify this subject into HC.To overcome the conflict of different classifiers, the D-S evidence model is used to fuse the information ofy^i1 and y^i2 [27]. The softmax function [28] is used in many machine learning applications for multiclass classifications to assign probability for each subject. The softmax function is expressed in the following form as Equation (15), which is used to calculate the probability assignment of each class.
(15)y^ij1∗=ey∧ij1∑j=13ey∧ij1,y^ij2∗=ey∧ij2∑j=13ey∧ij2,where y^ij1 is the element in the set y^i1 and j=1,2,3 represents HC, DS, and NDS, respectively. For example,y^i11∗is the probability that theith subject is HC according to the GMV feature, andy^i12∗is the probability that theith subject is HC according to the ALFF feature.Though softmax function, the probability assignment of theith subject can be obtained as
(16)y^i1∗=0.3458,0.2318,0.4224,y^i2∗=0.3672,0.3006,0.3322.Then, the combined strategy based on the D-S evidence is expressed as(17)y^i1∗⊕y^i2∗HC=1K×y^i11∗×y^i12∗,y^i1∗⊕y^i2∗DS=1K×y^i11∗×y^i12∗,y^i1∗⊕y^i2∗NDS=1K×y^i11∗×y^i12∗.The probability of subject into each groupy^i1∗⊕y^i2∗∈0,1, where K reflects the conflict level of evidences and can be represented as
(18)K=y^i11∗×y^i12∗+y^i11∗×y^i12∗+y^i11∗×y^i12∗.The result of the above example isy^i1∗⊕y^i2∗=0.3768,0.2068,0.4164 by Equation (17). According to the output of the fusion information, this subject will be classified into the class of NDS.


## 2.1. Classifier Based on the Xgboost Algorithm

The Xgboost algorithm is composed of many weak classification and regression trees (CARTs). Taking theith dataset (xi,yi) as an example, xi is the input variable with several attributes of fMRI data and yi is the real value of the given subject. For example, yi=1,0,0 means the ith subject is HC, yi=0,1,0 means the ith subject is DS, and yi=0,0,1 means the ith subject is NDS.Then, an Xgboost model can be mathematically expressed in the following form [25]:
(1)y^i=∑k=1Kfkxi,where K is the number of the CARTs, fk is the predicted value of each independent CART, and y^i is the predicted value with respect to input xi.The additive training model of Xgboost can be expressed as(2)y^i0=0,y^ik=∑t=1kftxi=y^kk−1+fkxi,where y^ik is the predicted value of the kth CART.The objective function of Xgboost includes a loss function and regularization term, which is expressed as(3)obj=∑i=1nlyi,y^i+∑k=1KΩfk,where lyi,y^i can be used to measure the error between predicted value and real value, n is the number of the subjects, and Ω is the regularization item to avoid overfitting.The specific form ofΩfk of the kth CART is given as
(4)Ωfk=γT+12λw2,where γ and λ present the penalty coefficients, T is the number of leaf nodes, and w is the weight of the leaf nodes.Then, the objective function of thetth step objt can be calculated by Equation (5) based on the previous step objt−1 based on the Equations (2) and (3).
(5)objt=∑i=1nlyi,y^it−1+ftxi+∑k=1t−1Ωfk+Ωft.By applying the second-order Taylor expansion to above equation, the objective function can be transformed into(6)objt=∑i=1nlyi,y^it−1+giftxi+12hift2xi+Ωft+Const.In the above expression,Const is a constant term at the step t; the parameters pair gi and hi can be calculated as
(7)gi=∂lyi,y^it−1∂y^it−1,hi=∂2lyi,y^it−1∂y∧it−12.Becauselyi,y^it−1is the constant item, the objective function can be rewritten as
(8)objt=∑i=1ngiftxi+12hift2xi+γT+12λ∑j=1Twj2+const′,where const′ is a new constant item at step t.According to the definition offk, fk can be written in the following form as
(9)fkx=wqx.Then, Equation (8) can be rewritten in the following form as
(10)objt=∑j=1T∑i∈Ijgiwj+12∑i∈Ijhi+λwj2+γT+const′,denoted by∑i∈Ijgi=Gjand∑i∈Ijhi=Hj; then, the objective function is expressed as
(11)objt=∑j=1TGjwj+12Hj+λwj2+γT+const′.The leaf nodes of thetth CART are each independent; Gj and Hj are the determined items. Then, minimizing the function equation (11), the optimal parameter wj can be calculated as
(12)wj∗=−GjHj+λ.Therefore, the final objective function is shown in the following form as(13)obj∗=−12∑j=1TGj2Hj+λ+γT+const′.The splitting algorithm [26] based on the above function is used to find the best split in Xgboost by
(14)Gain=12GL2HL+λ+GR2HR+λ−GL+GR2HL+HR+λ−γ.The gain function has four terms: the first two terms are the profits of left and right parts of a node, whereGL and GR are the left and right parts of Gj and HL and HR are the left and right parts of Hj; and the third item is the total profit of that node. The last item is the regularization item for preventing overfitting. The greedy algorithm determines whether a node obtains the maximum gain. Thus far, the optimal tree structure that maximizes the gain can be generated.The above description leads to the split finding algorithm for Xgboost presented as Algorithm1.Algorithm 1: Greedy algorithm for split finding.
Require:I, instance set of current node; d, feature dimensionEnsure: Split with max score1.Gain⟵0, G⟵∑i∈Ijgi, H⟵∑i∈Ijhi2.fork=1 to mdo3.GL⟵0, HL⟵04.for alli in sorted I by xjkdo5.GL⟵GL+gj, HL⟵HL+hj6.GR⟵G−GL, HR⟵H−HL7.score⟵maxscore,GL2/HL+λ+GR2/HR+λ−G2/H+λ8.end for9. end for


## 2.2. The Fusion Model Based on the D-S Evidence Theory

Analyses of the amplitude of low-frequency fluctuations (ALFF) and gray matter volume (GMV) are two important methods used in fMRI studies. Selecting GMV as the feature to construct the classifier based on the Xgboost algorithm, the predicted value of theith subject can be obtained and denoted as y^i1, while the other predicted value y^i2can be calculated by selecting ALFF as the feature,y^i1≠y^i2usually. Taking the ith subject as an example, the output information of the classifiers is y^i1=0.7,0.3,0.9 and y^i2=0.7,0.5,0.6. According to the output of classifier based on the feature of GMV, this subject should be classified into NDS, while the classifier based on the feature of ALFF will classify this subject into HC.To overcome the conflict of different classifiers, the D-S evidence model is used to fuse the information ofy^i1 and y^i2 [27]. The softmax function [28] is used in many machine learning applications for multiclass classifications to assign probability for each subject. The softmax function is expressed in the following form as Equation (15), which is used to calculate the probability assignment of each class.
(15)y^ij1∗=ey∧ij1∑j=13ey∧ij1,y^ij2∗=ey∧ij2∑j=13ey∧ij2,where y^ij1 is the element in the set y^i1 and j=1,2,3 represents HC, DS, and NDS, respectively. For example,y^i11∗is the probability that theith subject is HC according to the GMV feature, andy^i12∗is the probability that theith subject is HC according to the ALFF feature.Though softmax function, the probability assignment of theith subject can be obtained as
(16)y^i1∗=0.3458,0.2318,0.4224,y^i2∗=0.3672,0.3006,0.3322.Then, the combined strategy based on the D-S evidence is expressed as(17)y^i1∗⊕y^i2∗HC=1K×y^i11∗×y^i12∗,y^i1∗⊕y^i2∗DS=1K×y^i11∗×y^i12∗,y^i1∗⊕y^i2∗NDS=1K×y^i11∗×y^i12∗.The probability of subject into each groupy^i1∗⊕y^i2∗∈0,1, where K reflects the conflict level of evidences and can be represented as
(18)K=y^i11∗×y^i12∗+y^i11∗×y^i12∗+y^i11∗×y^i12∗.The result of the above example isy^i1∗⊕y^i2∗=0.3768,0.2068,0.4164 by Equation (17). According to the output of the fusion information, this subject will be classified into the class of NDS.


## 3. Experiments and Results

To evaluate performance of the proposed classification method, we scan 246 brain regions of 85 subjects by fMRI to extract the features of GMV and ALFF, including 16 DS, 31 NDS, and 38 HC. Obviously, it is a hard task to classify them in these small datasets. We apply the proposed algorithm to classify all the subjects.The present study classifies all subjects into three classes by applying linear regression (LR), supported vector classifier (SVC),K-nearest Neighbor (KNN), neural network (NN), naive Bayes (NB), classification and regression tree (CART), random forest (RF), extreme gradient boosting (Xgboost), deep neural network (DNN), and the proposed fusion algorithm. The flow of 10-fold cross-validation is shown in Figure 2. The results of 10-fold cross-validation of different classifiers are shown as below.Figure 2
Flow chat of 10-fold cross-validation.The receiver operating characteristic (ROC) curve is considered to evaluate the performance of classifiers. For different classification thresholds, the true-positive rate (TPR) (Equation (19)) is plotted against the false-positive rate (FPR) (Equation (20)). The area under the ROC curve (AUC) indicates the classifier’s ability to distinguish between classes. The value of the AUC is in the range [0,1]. AUC is 1 for a perfect classifier. In this work, the ROC curve is plotted for each class, as this is a multiclass problem. The microaverage and macroaverage are also computed by summing the individual values for true positive (TP), true negative (TN), false positive (FP), and false negative (FN). Then, the accuracy (Equation (21)), recall (Equation (22)), precision (Equation (23)), and F1-score (Equation (24)) are selected as the important metrics to evaluate the performance of different classifiers [29].
(19)TPR=TPTP+FN,(20)FPR=FPFP+TN,(21)Accuracy=TP+TNTP+TN+FP+FN,(22)Precision=TPTP+FP,(23)Recall=TPTP+FN,(24)F1‐score=2×precision×recallprecision+recall.The ROC of the above classifiers are shown in Figures3–8. In these figures, DS is the class 0, HC is the class 1, and NDS is the class 2. The key classification metrics are extracted from the results listed in Table 1.Figure 3
Results of logistic classifier (GMV).Figure 4
Results of logistic classifier (ALFF).Figure 5
Results of SVM with kernel of RBF (GMV).Figure 6
Results of SVM with kernel of RBF (ALFF).Figure 7
Results of Xgboost (GMV).Figure 8
Results of Xgboost (ALFF).Table 1
Classification metrics of different methods.
MethodsAccuracyPrecisionRecallF1-scoresAUCLR (GMV)64.8611%65.0509%64.8611%0.62370.8160LR (ALFF)65.9722%65.5394%65.9722%0.63880.7743SVC (GMV)61.5278%51.5602%61.5277%0.55080.7403SVC (ALFF)67.0833%59.4147%67.0833%0.61140.6114KNN (GMV)63.4722%61.4259%63.4722%0.58990.7174KNN (ALFF)58.6111%67.3175%58.6111%0.56010.6731NN (GMV)55.4167%56.4749%47.2222%0.54250.8066NN (ALFF)64.4444%67.0764%60.8333%0.62660.7618NB (GMV)53.8889%55.8657%53.8889%0.52160.7389NB (ALFF)68.8889%74.0139%68.8889%0.67810.8308CART (GMV)45.8333%47.7870%51.9444%0.47210.6375CART (ALFF)60%57.6075%56.6667%0.56670.6542RF (GMV)62.6389%60.0883%62.6389%0.58740.7722RF (ALFF)65.6944%57.1042%65.6944%0.59430.7969Xgboost (GMV)62.6389%60.4445%62.6389%0.58880.7958Xgboost (ALFF)68.0556%62.9610%68.0556%0.63640.8351DNN62.7083%40.3380%62.7083%0.51520.6723Our classifier (fusion)73.8889%65.4242%73.8889%0.67460.8524Due to the limited number of the DS, the AUC of class 0 is lower than other classes in the above figures. Therefore, it is hard to discriminate DS from NDS and HC. Taking Figure5 as an example, the macroaverage AUC score is 0.81 and the microaverage AUC score is 0.85 based on the feature of GMV by using the SVM classifier, which are better than the performance based on the feature of ALFF. The AUC score of HC is 0.85 based on the feature of GMV, which is better than the AUC scores of DS and NDS. The Xgboost classifier showed better performance than the SVM and logistic classifier. The microaverage AUC score is 0.90 and 0.80, respectively, on the feature of GMV and ALFF. The macroaverage AUC score is 0.90 and 0.84, respectively, on the feature of GMV and ALFF. The ROC of the fusion algorithm is shown in Figure 9. In this figure, we find the AUC score of the class 0 is 1, which means it is a perfect classifier. This classifier combines the advantages of the classifier based on feature of GMV and ALFF. Therefore, the performance of this fusion classifier is much better than others.Figure 9
Results of Xgboost (fusion).From Table1, 10-fold cross-validation showed our algorithm with an accuracy of 73.89%, which is higher than other classifiers. Many metrics of the proposed classifier are better than other classifiers. In many present studies, the accuracy of the present classifier is usually less than 70% when the datasets are small. In this paper, the proposed fusion classifier will improve the performance effectively for the small datasets by combining the advantages of each feature. The accuracy of different classifiers is shown in Figure 10.Figure 10
Accuracy of different classifiers.


## 4. Conclusion

In this paper, a new multiple classifier method was proposed for the small datasets and applied to discriminate the two subtypes of schizophrenia and health controls based on the fMRI data. Due to the limitation data and indexes, this study constructed the Xgboost algorithm based on the different features. To improve the accuracy, the fusion model was used to combine the information from different classifiers. Finally, the subject would be classified into the class with the maximum probability. This method was applied to classify 38 healthy controls, 16 deficit schizophrenic patients, and 31 nondeficit schizophrenic patients. 10-fold cross-validation showed our algorithm with a sensitivity of 73.89%, which was much higher than other classification algorithms when the datasets were small. In addition, the proposed algorithm can be used to discriminate different classes for the large datasets. The performance of the proposed algorithm would be effective than other algorithms when the datasets are small. It will bring the better performance in diagnosing subtypes of schizophrenia. Although the findings in our study are rigorous, there are some limitations: (1) relatively small sample size; (2) interference caused by antipsychotic drugs during experiment; and (3) limitations of the algorithm itself. In the future work, more subjects will be collected in the project, including different subtypes and HC. The original image data should be obtained in the processing of experiments, and more deep learning approaches will be proposed to solve this multiclass classification problem.



---

*Source: 1581958-2022-07-19.xml*