A Hybrid Crow Search and Grey Wolf Optimization Technique for Enhanced Medical Data Classification in Diabetes Diagnosis System

Diabetes is an extremely serious hazard to global health and its incidence is increasing vividly. In this paper, we develop an effective system to diagnose diabetes disease using a hybrid optimization-based Support Vector Machine (SVM).The proposed hybrid optimization technique integrates a Crow Search algorithm (CSA) and Binary Grey Wolf Optimizer (BGWO) for exploiting the full potential of SVM in the diabetes diagnosis system. The effectiveness of our proposed hybrid optimization-based SVM (hereafter called CS-BGWO-SVM) approach is carefully studied on the real-world databases such as UCIPima Indian standard dataset and the diabetes type dataset from the Data World repository. To evaluate the CS-BGWO-SVM technique, its performance is related to several state-of-the-arts approaches using SVM with respect to predictive accuracy, Intersection Over-Union (IoU), specificity, sensitivity, and the area under receiver operator characteristic curve (AUC). The outcomes of empirical analysis illustrate that CS-BGWO-SVM can be considered as a more efficient approach with outstanding classification accuracy. Furthermore, we perform the Wilcoxon statistical test to decide whether the proposed cohesive CS-BGWO-SVM approach offers a substantial enhancement in terms of performance measures or not. Consequently, we can conclude that CS-BGWO-SVM is the better diabetes diagnostic model as compared to modern diagnosis methods previously reported in the literature.


Introduction
Precisemedical diagnosis and effective disease management are two significant problems in life science analysis and have a constructive effect on the health condition of individuals. At present, the majority of healthcare institutions are adopting innovative technologies for their administrativeendeavor to enable an improved healthcare service in terms of techno-economic aspects [1]. Therefore, often colossal and devastating volumes of data (i.e., petabytes or terabytes) are collected from this medical industry. The potential and usefulness of big data in the medical sector, such as disease management and providing an effective treatment are realized from its impending nature of the appropriate methods to intelligently discover knowledge and build models from data [2]. Intellectual database management systems provide efficient tools to extract knowledge from this complex and massive data and convert all extant data into useful information to develop an intelligent decision support system. Diabetes has been a serious global health issue for many years. Innumerable worldwide and nationwide epidemiological studies have witnessed the increasing number of diabetic patients globally. It is a non-communicable disease categorized by hyperglycemia (i.e., high blood glucose) and associated with disorders of carbohydrate, protein, and fat. Besides, it deteriorates nearly all organs of the body including the kidney, heart, foot, eyes, skin, and so on. The diabetes managing applications demand unremitting clinical services to decrease the jeopardy of long-term problems and to thwart deadly consequences.
In the last few decades, health care professionals pay much more attention to an early diagnosis of diabetes disease and design several approaches to resolve these diagnosis issues. Exploiting machine learning techniques in the Diabetes Diagnosis System (DDS) will help the medical professional to provide timely care for patients by diagnosing the disease at the initial stage. Researchers have employed artificial intelligence techniques to design some DDS, which enhance the performance of the diabetes management system. Recently, several research works apply the SVM method to find out diabetes [3][4][5]. SVM is a classifier based on a supervised machine learning algorithm and capable of differentiating the 'intrinsic features' of different data samples for nonlinear problems. Nevertheless, the effectiveness of the SVM classification is closely associated with the feature of the selected parameters. Hence, parameter optimization techniques perform a vital role in SVM performance and affect the classification accuracy of DDS. The SVM-based classifier consists of two processes namely, training and testing. This approach recognizes the metrics for classification using known datasets in training phase. In the testing phase, it carries out classification on the basis of those metrics and labels them to the equivalent class. In the present work, we aim to achieve the diabetes classification using a hybrid optimization technique-based SVM classification.
While developing and implementing a meta-heuristic, it should be remembered that a suitable balance between exploitation and exploration must be realized effectively. By the way, one important alternate solution is to design a hybrid approach in which (at least) two algorithms are assimilated together to increase the system performance. The key objective of this study is to develop and design a cohesive optimizer aiming at the assimilation of CSA and BGWO to establish higher classification accuracy on real-world problems. CSA is a meta-heuristic optimizer that imitates the intellectual behaviour of crows. In this study, we utilize the cleverness of crows to store food in a hiding place and retrieve it whenever needed. The merit of CSA lies in the aptitude to evadethe trapping in the local optimum efficiently when coping with multimodal optimization problems in a more complex searching space. Nevertheless, the exploitation phase of CSA is not much effective.
Grey wolf optimizer (GWO) is also a meta-heuristic optimizer that imitates the social behaviour of grey wolves such as leadershiphierarchyand the group hunting. GWO has the superior abilities to circumvent algorithms getting stuck into local optimal value [6]. It also has a fast convergence speed in finding the optimum solution. Generally, GWO is good at o exploitation. But, it suffers fromthe deprived global searching ability. Therefore, in certain circumstances, it cannot handle the problem effectively and flops to calculate the global optimum solution. This study develops a binary grey wolf optimization technique for handling parameter optimization problems, and then SVM is used to achieve diagnostic processes according to the selected features recognized by the BGWO algorithm.
Considering the virtues of CSA and BGWO, these optimization techniques are perfect for hybridization. The hybrid optimizer exploits the strengths of both CSA and BGWO techniques successfully to provide promising candidate solutions and realize global optima effectively. The proposed model adopts CSA to select the better initial position and then implements BGWO to appraise the position of each searching agent for obtaining the optimum parameters with better accuracy. Hence, the proposed approach integrates these two algorithms to balance exploitation and exploration appropriately and provide improved performance than the traditional CSA and GWO with respect to accuracy, IoU, specificity, AUC, and sensitivity of the classification. The key contributions of our study are is three-fold. The remainder sections of the article are organized as: We explore related hybrid optimization methods in Sect. 2. We describe the fundamental notions of SVM, BGWO, and CSA in Sect. 3. A comprehensive depiction of the presented method is presented in Sect. 4. The experimental setup, evaluation metrics and numerical results are given in Sect. 5. Finally, the research is concluded at Sect. 6.

Related Work
Recently, several different hybrid metaheuristics approaches have been proposed and studied to resolve various complicated optimization issues [7]. Hybrid algorithms have become very popular due to their potential to manage the high-dimensional nonlinear problems of optimization [8]. Hitherto, extensive investigations have been carried out by several researchers on hybrid optimization approaches. Kao and Zahara proposed a novel optimization technique in which the Genetic Algorithm (GA) has been integrated with Particle Swarm Optimization (PSO) for hybridization [9]. To generate individuals, this hybrid technique not only used cross-over and mutation operators but also used searching operators of PSO. Tsai et al. introduced Taguchi-genetic algorithm-based hybridization to reach a global solution [10]. In this study, the Taguchi method is combined with GA to select higher quality genes to achieve better performance.
Jitkongchuen presented an enhanced mutation method using hybrid differential evolution approach. In this study, a differential evolution algorithm is combined with GWO for handling continuous global optimization problems [11]. Nabil proposed an improved version of the flower pollination algorithm (FPA). This work combines FPA with a clonal selection approach to achieve better accuracy than the traditional FPA [12]. The effectiveness of the modified FPA is evaluated via different benchmark datasets. Tawhid and Ali developed an innovative cohesive method based on GWO and GA. The authors applied this approach to reduce potential energy function [13]. Jayabarathi et al. developed a combined GWO approach by as simulating the original GWO with GA for getting a better efficiency for handling economic dispatch problems [14]. Singh and Singh hybridized the local search facility of PSO with the global search facility of GWO to improve efficiency and computational efforts [15]. Gaidhane and Nigam proposed a new hybrid approach to optimize the constraints for increasing the performance of complicated problems. This work integrates the bee's communication tactic from Artificial Bee Colony (ABC) algorithm, to realize effective exploration, with the leadership hierarchy of GWO retaining its exploitation ability [16].
Zhang et al. hybridized biogeography-based optimization (BBO) and GWO to make a suitable trade-off between local search and global search to achieve superior results than the original implementation of BBO and GWO. Hassanien et al. introduced the CSA to provide the global optimum result precisely. The proposed approach can effectively handle the roughness and impreciseness of the prevailing statistics about the global optimum and eventually increase the efficiency of the original CSA [17].The aforementioned research works proved that the hybrid algorithms have higher performance related to other optimization methods. In the domain of feature extraction, several combined methods are developed. In the year 2004, Oh et al. presented the first cohesive metaheuristic method for attribute selection which intended to hybrid local search methods and GA to control the searching process [18]. Talbi et al. proposed another hybrid approach for feature selection which intended to integrate PSO, GA, and SVM classification model and applied to microarray data classification [19].
Hybridization in ant colony optimization using the cuckoo search [20], and using GA [21,22] are also developed for effective feature selection process. Mafarja and Mirjalili developed a cohesive method to combine simulated annealing methods with the whale optimization approach for identifying the ideal set of features [23]. In recent times, a technique using evolutionary population dynamics [24] and grasshopper optimization [25] has been effectively implemented in [26] for handling optimal feature selection problems. For further information regarding feature selection and optimization approaches, interested readers can refer to [27] and [28].

Background Analysis
This section presents some basic concepts of SVM, BGWO, and CSA that have been utilized to construct and implement the proposed CS-BGWO-SVM classification model.

Support Vector machine
It is a discriminative classifier developed by Vapnik to find anomalies of biomedical signals owing to its excellence and strong ability to deal with the non-linear and high-dimensional dataset in the field of the healthcare industry [29].
The key notion behind this classifier is to differentiate the unknown validation set of data in to their appropriate classes based on the training set of some well-known data.
In binary classification, SVM constructs a hyperplane that optimally differentiates samples into two classes. Let is the training dataset. Here x denotes the input sample, and y ∈ {+1, −1} is the class label. Then the hyperplane is denoted as In above Eq. (1), a coefficient vector w is orthogonal to the hyperplane. The term b denotes the distance between the point in the dataset and origin. The key objective of the SVM is to determine b and w values. To get the ideal hyperplane ‖ ‖ w 2‖ ‖ should be minimized under the condition of y i w T .x i + b ≥ 1 as shown in Fig. 1. Hence, the optimization problem is formulated as In a linear problem, Lagrange multipliers are used to calculate w. The data points lay on the decision boundary are named as support vector. Therefore, the value of w can be calculated as To handle a nonlinear problem, SVM employs kernel trick. Then the decision function can be given by here Ŷ is the kernelized label for the unlabeled input x and the sgn function defines whether the anticipated classification comes out positive or negative. Usually, any positives definite kernel functions such as Gaussian function [30]. This subsection only gives a brief note on SVM. For further reading, readers can refer to [31] which give a comprehensive depiction of the SVM concepts.

Binary GWO
The conventional greywolf optimizer is a bio-inspired metaheuristic approach introduced by Mirjalili et al. [32]. GWO imitates the leadership level and hunting actions of grey wolves. Grey wolves are Canidae species and have a rigorous social dominant hierarchy. They are famous for their swarm intelligence-based hunting mechanism. Grey wolves are social predators or pack hunters which hunt their prey (target) by working together in a group of 5 to 12 wolves (i.e., pack). Within the pack, the leadership is divided into four levels like alpha (α), beta (β), delta (δ), and omega (ω). α(dominant) wolf is at the upper most hierarchy and acts as the leader of the cluster. Alpha is liable for making all the verdicts about hunting, sleeping place, wake uptime, maintaining discipline, etc. [33]. Its verdicts are dictated to and followed by the entire group. The wolves in the next level are beta which acts as an advisor to alpha and fortifies the orders of alpha all over the group. This means that β aids α in decision making or other tasks and dominate other lower-level wolves. Beta wolves being the subleaders in the group have the maximum possibility to become α wolves if one of α becomes very old or dead. The subsequent hierarchy is δ. They have to bow to α as well as β; in chorus, they instruct next lower-level wolves (i.e., ω wolves). Hunters, sentinels, scouts, caretakers, and elders belong to this level. Hunters support α and β when hunting targets and providing foodstuff for the group. Sentinels protect and assure the security of the group. Scouts inspect the borders of the search space and alert the group in case of any hazard. Caretakers are in charge of looking after the injured, sick, and weak wolves in the group. Finally, older wolves are experts who used to be α or β.
The wolves present at the lowermost level in the dominance structure are omega. They act as a scapegoat. They always subordinate to other higher-level wolves. The omega wolves are permitted to consume food last of all. It looks like ω is not an important wolf in the group, but it has been witnessed that the entire cluster realizes some complications and in-house fighting when dropping the ω wolf. This is owing to the expelling of frustration and violence of all wolves by ω. This helps in satisfying the whole group and preserving the hierarchy. In some scenarios, ω is also the baby-sitters in the group. Besides this dominant structure, pack hunting is an additional remark able characteristic of grey wolves. As stated by Munro et al. the basic steps of group hunting are given by [34]: (i) to track, chase and grasp the target (ii) to enfold and harass the target until it becomes stable, and (iii) to attack and kill the prey. Figure 2 illustrates these steps.
The following equations are used to define the encircling nature of the wolves [32].
where i specifies the current iteration. ⃗ P prey and ⃗ P SA are the position vectors of prey and search agent correspondingly. ��⃗ M and ⃗ K are coefficient vectors and are estimated as The classification process of SVM cant role in circumventing stagnation in local solutions. It appends some arbitrary weight to the target to make it a challenge for searching agents to denote the distance between itself and the target. If �� ⃗ K > 1, then �� ⃗ K is maximizing the impact of the target, and if �� ⃗ K < 1 , its influence will get minimize in a stochastic manner. In the entire course, ��⃗ M and �� ⃗ K are attuned cautiously and the exploitation or exploration processes get emphasized or deemphasized. Eventually, when a definite condition is met the algorithm will be halted and the optimum value of α will be returned.
Indeed, the original grey wolf optimizer is introduced to handle the continuous optimization problems. In the case of discrete optimization problems, a binary form of the optimizer is essential. This work implements a binary variant of GWO for the feature extraction process. Each searching agent in the proposed BGWO consists of a flag vector and its size is as same as attributes. As soon as the position of each searching agent was updated by Eq. (17), its discrete position is expressed as follows.
where, P i,j represents the j th place of the i th wolf. Figure 3 displays the algorithm of BGWO.

Crow Search Algorithm
CSA imitates the smart performances of crows in hiding and stealing foodstuffs. It is a metaheuristic approach presented by Askarzadeh. It is widely used in several scientific solicitations particularly in the field of optimization techniques [35]. Crows are now considered to be the world'sbest genius birds. Crows have the largest ratio of brain-to-bodyweight. They can recognize individuals and inform each other when a hostile one reaching. Likewise, they can exploit techniques, exchange information in matured ways, and recollect their feeding positions (places) across seasons [36]. Crows are famous for identifying their feeding positions, monitoring other crows, and stealing the food stuffs when the proprietor leaves. After a crow has executed larceny, it will provide additional fortifications such as shifting their feeding places to evade being an imminent target. They apply their intelligence of having been a thief to analyse the actions of other birds and can discover the optimum method to defend their foodstuffs from being taken off [37]. The rudimentary conventions of this CSA are (i) crows are reside in groups; (ii) crows reminisce their feeding positions; (iii) they monitor other crows to take foodstuff; and (iv) crows safeguard their feeding positions from being theft.
Consider a d-dimensional search area with a group of an N number of crows. The hiding place of the crow x at iteration i is denoted by P Here I max is the maximum iteration. All the crows have a memory in which their optimal place has been recalled. At iteration i , the place of search agent x is defined by P x,i . This is the optimal position that x has achieved hitherto. Furthermore, crows shift their place and search for enriched sources of food. Assume crow y desires to visit its feeding position,P y,i . At iteration i , x decides to monitor y to approach P y,i . Now two scenarios can arise.
(i) y is not cognizant of being monitored by x . Consequently, crow x will reach the feeding place of y . Now, the new place of x is measured by Eq. (20).
where r and is a randomly selected value with uniform distribution taking values in [0,1] and F x,i represents the flight length of x at iteration i. (ii) y is aware of being monitored by x . Hence, y will cheat crow x by flying to a new place in the search space to defend its foodstuff from being stolen.

Now, cases 1 and 2 becomes
where rand y represents an arbitrary value and PA y,i denotes the probability of awareness (PA) of y at iteration i . Usually, meta-heuristic algorithms should deliver an adequate trade-off between intensification and diversification [38]. In this algorithm, diversification and intensification are mostly regulated by PA . By decreasing PA , this approach attempts to perform the search on a local space. In chorus, for the lesser value of PA the intensification is augmented. On the other hand, by increasing PA , the possibility of searching the space for obtaining better solutions decreases, and this approach is apt to global search (randomization). Accordingly, the application of higher values of PA upturns the diversification. Figure 4 shows the pseudo-code of CSA.

Proposed CS-BGWO algorithm
For any population-based approach, the exploitation and exploration phases are mutually important to achieve better performance. In BGWO, the main problem is that all the locations of searching agents are rationalized according to the value of α, β, and δ in the entire process as given in Eq. (17). Mostly, this updating procedure causes early convergence since agents are not permitted to discover the search space effectively. Furthermore, the optimization process as given in Eq. (17) offers limited local search skills in the later part of the process which causes relaxed convergence. Therefore, to overcome these restrictions in BGWO is integrated with the crow search algorithm to realize an appropriate trade-off between local and global search. CSA includes a regulating factor F x,i in its position updating procedure as given in Eq. (21) which enables the agents to adapt the scope of the step movement in the direction of the new agents. This factor is significant in achieving the global optimum value as the higher value of F x,i result to global search where as a lesser value of F x,i leads to local search. As stated in the previous section, BGWO has better local search capability however deprived global search ability, hence in the CS-BGWO algorithm, a higher value of F x,i is used to applying with the outstanding exploration feature of CSA as given in Eq. (22). This denotes that CS-BGWO can effectively utilize the proficiencies of two algorithms and thus, it is universally applicable. In the CS-BGWO algorithm, an agent is permitted to appraise its location only by α and β rather than apprising from α, β, and δ as given in Eq. (22).
To preserve population diversity, instead of allowing all the search agents to update their position from α and β wolves, only the value of α is exploited in the CS-BGWO algorithm. This reduction approach facilitates the proposed algorithm to avoid local optimum effectively.
Although the CS-BGWO approach has the outstanding competencies of local and global search an appropriate trade-off between these two phases must be realized to obtain better results. In an ideal situation, an optimization technique can search a large space in the initial phases of the process to evade early convergence and use a small search space in the later phases of the iteration to effectively refine the results.
To assess the efficiency of the developed approach, more instinctively, the convergence curve of the developed approach is related to other approaches found in the literature. In Fig. 5, the X-axis designates the iteration counts, and the Y-axis indicates the best score. In an ideal situation, an optimization technique can search a large space in the initial phases of the process to evade early convergence and use  Figs. 5a, b, and c. The convergence speed of CS-BGWO is lesser than other approaches in the earlier phase, this indicates that the proposed approach has robust global search ability than other approaches in the initial phase, and the approach is difficult to drop into the local optimum result. Furthermore, the CS-BGWO algorithm still preserves a robust search aptitude in the later phase.
From the viewpoint of the complete iteration, the global and local search abilities of the CS-BGWO approach have been well balanced, which shows that the hybridization of CSA and BGWO can successfully enhance the performance of the classification process, and is more suitable for resolving diabetes diagnosis issues than other related algorithms. This indicates that to achieve the necessary global-local search ratio, a constant balance probability among Eq. (22) and Eq. (23) is not constructive. Thus the present study implements an adaptive balance probability ( ) that enables CS-BGWO to improve the accuracy and efficiency of searches. The value of can be calculated as follows where I max denotes the maximum iteration counts and i signifies the present one. This adaptive balance probability for various iterations is depicted in Fig. 6.
It is worth mentioning that ��⃗ M acts as a regulating element to balance local and global search. The value of ��⃗ M depends on �� ⃗ m which eventually regulates the direction of the searching. The higher value of �� ⃗ m enables global search and its smaller value enable the local search. This indicates that an apt selection of �� ⃗ m provides a good trade-off between local and global search which results in higher performance. In the BGWO approach, the value of �� ⃗ m is gradually reduced to 0 from 2. Hence, it can be witnessed that superior results can be realized if the value of �� ⃗ m is decreased non-linearly. By exploiting this notion, an enhanced approach, as given in Eq. (25), is developed to select the values of �� ⃗ m . This tactic enables the CS-BGWO approach to discover the search space efficiently as compared with the original GWO. Figure 7 shows the pseudo-code of the CS-BGWO algorithm.

CS-BGWO-SVM model
In this work, a novelCS-BGWO-SVM framework is developed to handle the feature selection problem in DDS. CS-BGWO-SVM model consists of two phases: (i) CSA is employed to define the initial locations of searching agents, and then BGWO is implemented to choose the optimum feature subset by exploring the search space adaptively. The optimum feature extraction process applies to the selected (25) m = 2 − cos(rand x ()) 1 I max Fig. 6 Adaptive balance probability Fig. 7 Algorithm of CS-BGWO features to increase the classification accuracy and reduce the number of designated attributes simultaneously. The mean accuracy of the classification over the tenfold crossvalidation (CV) scheme is employed as a fitness function for evaluating the selected feature subset; and (ii) the SVM classification model is simulated to achieve superior predictive accuracy according to the best feature subset. The dataset has been normalized in the range of[− 1, + 1] before classification. To achieve a more accurate estimate of a model's performance, the k-fold CV is employed [39]. The proposed work assumes k equals to10. This means that the entire dataset is partitioned into ten subsections. For each iteration, one subsection is considered as the testing dataset and the other subsections are integrated to construct the training dataset. Then the mean value of error across all 10 independentruns is calculated. The merit of this technique is that all testing data are independent and the consistency of the outcomes could be improved. It is noteworthy that only one replication of the tenfold CV will not produce adequate results for evaluation owing to the randomness in data splitting. Therefore, all the results are stated on an average of ten trials to gain precise estimation [40].

Experimental Setup
In our assessment, the experiments are carried out on 3.6 GHz with 16 GB RAM, Intel Core i7-4790 processor with Windows 10 operating system. All the algorithms are simulated using MATLAB R2009b. The proposed classification approach uses the Pima Indian diabetes patient records in the UCIand the diabetes type dataset in the Data World repository. Pima Indian diabetes patient records are collected from the National Institute of Diabetes and Digestive and Kidney Diseases. This database involves 768 medical data of Pima Indian heritage, inhabitants existing near Phoenix, Arizona, USA [48]. There are eight features related to this dataset (i.e., plasma glucose concentration, 2-h serum insulin, the number of times pregnant, diastolic blood pressure, function of diabetes nutrition, triceps skin fold thickness, body mass index (BMI), and age). Table 1 shows the statistical report of each instance in this dataset. The range of binary variables is limited to '0' or '1'. The target variable '1' represents a positive result for diabetes disease (i.e., diabetic), '0' is a negative result (i.e., non-diabetic). After preprocessing dataset, there are 392 cases with no missing values. The number of cases in label '1' is 130, and the number of cases in label '0' is 262.
The patient records in diabetes type dataset are collected from the Data World repository [49]. This database consist of 1099 records with 8 attributes including Plasma glucose test randomly taken at any time, plasma glucose test usually taken in the morning or 8 h after a meal, blood sugar while fasting, blood sugar 90 min after a meal, HbA1c (glycated hemoglobin), type, class, and age. Table 2 shows the statistical report of each instance in this dataset. After preprocessing dataset there are 1009 cases with no missing values. The number of cases in label '1' is 653, and the number of cases in label '0' is 356.
Selecting the penalty factor (C) and kernel width (γ) is extremely significant for implementing SVM models. The value of C affects the classification accuracy. If the value of C is very high then the accuracy is maximum in the training process, but it is minimum in the testing process. If C has a lesser value, then the accuracy is unacceptable, sorting the approach impractical. The value of γ has a greater impact on accuracy than C due to its influence on classification performance. An infinitesimal value of γ leads to under-fitting whereas a very high value results in over-fitting [50]. We assume both parameters are in the range of {2 −5 , 2 −4 ,…, 2 4 , 2 5 }. This study assumes C and γ as 32 and 0.125, respectively (i.e., C = 2 5 (32) and γ = 2 −3 (0.125)). The values of both parameters are selected through the trial and error technique. Table 3 shows the detail of the parameter setting.
Our proposed algorithm encounters challenging problems, such as outlier data and noise, which can reduce accuracy. The outlier detection method can be utilized in the pre-processing step to identify inconsistencies in data/outliers; thus, a good classifier can be generated for better decision-making. Eliminating the outliers and noise from the training dataset will enhance the classification accuracy. Past literature showed that by removing noise the quality of real datasets is enhanced [51]. In this study, DBSCAN is employed to identify the outlier data in a diabetes dataset [52]. The goal is to find objects that close to a given point to create dense regions. The points that are located outside dense regions are treated as outliers.
The resultant dataset after removing outlier is given in Table 4.

Evaluation MEASURES
To assess the efficacy of DDS using the CS-GWO-SVM model, we studied four important performance measures: classification accuracy, AUC [53], specificity, and sensitivity. These evaluation parameters are essential to be greater to increase the efficiency of the presented approach. The utility is calculated with respect to the accuracy as The intersection over-union (IoU) is a performance measure used to calculate the performance of any data classification technique. Given a data set, the IoU metric provides the similarity between the predicted data and the ground truth for a data existing in the dataset. The IoU measure can consider the class imbalance issue generally present in such a problem setting. It is defined as follows Sensitivity and specificity denote in what manner the classification approach differentiates negative and positive cases. Sensitivity defines the rate of disease prediction and specificity denotes the rate of false alarm.
In the above equations, true positive (TP) denote the number of persons properly categorized as diabetic patients; false negative (FN) denotes the number of diabetic patients wrongly identified as the non-diabetic. The term true negative (TN) denotes that the number of persons properly categorized as non-diabetic and false positive (FP) defines the number of non-diabetic person wrongly categorized as the diabetic one. The receiver operator characteristic (ROC) is a pictorial representation to define the classification accuracy. The curve shows the rate of TP and FP. The AUC represents the area under the ROC curve. This is an important measure for relating binary classification problems.
The Wilcoxon statistical test is carried out to decide whether the proposed CS-BGWO-SVM algorithm delivers     a noteworthy improvement related to other approaches or not [54]. The test was performed using the effects of the proposed CS-BGWO-SVM and related to each of the other approaches at a level of 5% statistical significance. Table 10 depicts the p values achieved through the test, where the p value < 0.05 indicate that the null hypothesis is precluded, i.e., there is a considerable variance at a level of 5% significance. In contrast, the p values > 0.05signify that there is no notable variance between the compared values. It can be evaluated from the results depicted in Table 10. From the results, it is observed that most of the p values are less than 0.05 which confirms that the improvement achieved by the CS-BGWO-SVM is statistically significant.

Experimental Results and Discussion
The proposed CS-GWO-SVM classification model is implemented. Table 5 shows the comprehensive results achieved by the proposed model. To facilitate the effectiveness of the proposed CS-GWO-SVM classification model, we relate the performance of the proposed approach to other prominent algorithms from the literature. The experimental results of various approaches are reported in Table 6. From Table 6, we can observe that the conventional SVM approach has realized nominal classification performance with 81.1% ACC, 86% IoU, 0.829 AUC, 81.4% sensitivity, and 85.4% specificity. With the intention of further increase in efficiency of the SVM approach, we combined optimization algorithms such as CSA, GA, PSO, GWO, EGWO, AGWO, and CS-BGWO with the original SVM. These models are implemented and examined rigorously on the same data. In the CSA-SVM classifier, PA is directly employed to regulate the algorithm diversity. Hence, this classifier has generated reasonable outcomes related to the original SVM. This algorithm has realized classification results with 87.1% ACC, 82.6% IoU, 0.869 AUC, 85.3% sensitivity, and 88.3% specificity. The CSA-SVM and GA-SVM approaches provide similar results; however, the CSA-SVM approach provides improved standard deviation as compared to GA-SVM. Since these classifiers hinge on the arbitrary generation of individuals and always there is a possibility to create a zero variable vector. Figure 8 illustrates the comprehensive outcomes acquired by CS-BGWO-SVM for various designated feature subset.
The benefits like fewer constraints and random walks in GWO enable an effective optimization approach for SVM classification. GWO-SVM realized classification results with 86.1% ACC, 87.6% IoU, 0.856 AUC, 83.5% sensitivity, and 87.7% specificity. EGWO-SVM with a better hunting mechanism has realized better performance with 89% ACC, 89.3% IoU, 0.861 AUC, 84.2% sensitivity, and 89.7% specificity by developing a superior trade-off between local and global search that provides improved outcomes. AGWO-SVM hit better results with 88% ACC, 89% IoU, 0.877 AUC, 85.7% sensitivity, and 91% specificity due to its exploitation ability. It is possible to conclude that the CS-BGWO-SVM model has achieved much better performance with 93.7% ACC, 95.2% IoU, 0.923 AUC, 90.4% sensitivity, and 94.9% specificity. As compared with the original SVM model, CS-GWO-SVM has boosted 12.6%, 9.4%, 9.0%, and 9.5% with respect to ACC, AUC, SEN, and SPE correspondingly. Also, it is stimulating to observe that the SD gained by the CS-GWO-SVM is lesser than that of all most all other classification methods which means that the CS-GWO-SVM can provide more robust diagnose results.
The mean value of performance measures and SD gained from Pima Indian dataset by each approach is illustrated in Figs. 9 and 10. It can be found that the CS-GWO-SVM outdoes all other approaches with respect to the performance metrics. The key reason behind the superior performance of CS-GWO-SVM is that the CSA-based initialization can increase the efficiency of the feature selection process to a certain level. From Fig. 10, it can be observed that the SD of the CS-BGWO-SVM was less than all other approaches in terms of the performance measures. Therefore, CS-BGWObased SVM model provides much more consistent results for diagnosing diabetes than the others. In other words, the CS-BGWO can not only enhance the enactment of SVM, but also enable better results for diagnosing diabetes disease. The comparative study demonstrates that CS-BGWO-SVM is a very competitive approach for diagnosing diabetes disease.   We can obtain similar results when we apply our proposed algorithm on diabetes type dataset. The experimental results obtained from diabetes type dataset using various approaches are reported in Table 7. The mean value of performance measures and SD gained from diabetes type dataset by each approach is illustrated in Figs. 11 and 12. From this table, we can observe that the conventional SVM approach has realized nominal classification performance with 81.3% ACC, 82% IoU, 0.839 AUC, 85.1% sensitivity, and 85.6% specificity. In the CSA-SVM classifier, PA is used to control the algorithm diversity. Therefore, it has produced reasonable outcomes related to the original SVM. This algorithm approach has realized classification results with 87.3% ACC, 78.8% IoU, 0.879 AUC, 89.2% sensitivity, and 88.5% specificity. The CSA-SVM approach provides improved standard deviation as compared to GA-SVM. Since these classifiers hinge on the arbitrary generation of individuals and always there is a possibility to create a zero variable vector.
As compared with CSA-SVM, GA-SVM, and PSO-SVM, IWOA has higher data classification accuracy better accuracy (87.6%), IoU (83.6%), AUC (0.879), sensitivity (90%), and specificity (89.4%), with minimum prediction SD. GWO-SVM realized classification results with 86.1% ACC, 87.6% IoU, 0.856 AUC, 83.5% sensitivity, and 87.7% specificity. EGWO-SVM with a better hunting mechanism has realized better performance with 89.2% ACC, 85.2%  .1% specificity. Also, it is observed that the SD gained by the CS-GWO-SVM is lesser than that of all most all other classification methods which means that the CS-GWO-SVM can provide more robust and reliable diagnose results. Figures 9, 10, 11, 12 also illustrate the p values of CS-BGWO-SVM in comparison with some meta-heuristic approaches achieved by the rank-sum test. This assessment is performed to decide whether the variation between the outcomes of CS-BGWO-SVM and the results of other approaches is noteworthy or not. Precisely, a p value is calculated, and if it is less than 0.05 then it designates that the results attained by CS-BGWO-SVM have considerable deviations from those of the other approaches. The p value > 0.05 denotes that there is no substantial variation between the outputs of CS-BGWO-SVM and the other approaches used for evaluation. In most cases, it can be easily found that the p values calculated by the Wilcoxon test are less than 0.05 which reveals the outstanding performance of CS-BGWO-SVM. Tables 8,9,10,11,12,13 illustrate the results of all the approaches for various folding on Pima Indian dataset. The mean and SD results gained by each approach are stated in each table and the optimum numerical values are emphasized in bold. It is found that the performance measures obtained by the CS-BGWO-SVM approach are superior to all other approaches. Hence, the results show that the integration of CSA and BGWO has revealed better results as compared to the original GWO, CSA, and all other approaches used in this research work. Moreover, it is worth noting that CS-BGWO-SVM exhibits better results as related to EGWO and AGWO in the majority of the scenarios. This demonstrates that the integration of CSA and BGWO has considerably increased the exploration aptitude and robustness.

Conclusion
Data classification in the medical diagnosis system supports physicians to learn hidden data samples by training a huge amount of real-world datasets. Diabetes mellitus is a severe global health problem. As stated in theInternational Diabetes Foundation report, now there are 300 million diabetic patients across the world, and an additional 300 million individuals are forestalling to be at greater possibility of diabetes in the year 2030. Therefore, it is indispensable for timely prediction, analysis, and management of hyperglycaemi a and its consequences. We develop an effective system to diagnose diabetes disease using a hybrid optimization-based SVM. The proposed hybrid optimization technique integrates a CSA and BGWO for exploiting the full potential of SVM in the diabetes diagnosis system. The hybrid optimizer exploits the strengths of both CSA and BGWO techniques successfully to provide promising candidate solutions and realize global optima effectively. In this work, a better initial seed is selected by the CSA and then BGWO is used to appraise the locations of the searching agent in the discrete searching area for obtaining the optimum parameters with better accuracy. The effectiveness of our proposed hybrid optimizationbased SVM (CS-BGWO-SVM) approach is rigorously evaluated on the technique, its performance is real-world database. To evaluate the CS-BGWO-SVM technique, its performance is related to many advanced metaheuristicbased SVM approaches. The empirical analysis illustrates that CS-BGWO-SVM can be considered as a promising approach with outstanding performance measures.
Funding Authors received no specific funding for this study.

Availability of Data and Material
The dataset for the experimentation is obtained from the National Institute of Diabetes and Digestive and Kidney Diseases. This dataset consists of 768 clinical records of Pima Indian heritage, a population living near Phoenix, Arizona, USA. There are eight attributes associated with this database (i.e., number of times pregnant, plasma glucose concentration estimated in 2 h after consuming a 75 g oral glucose, diastolic blood pressure, Triceps skinfold thickness, 2-h serum insulin, body mass index, diabetes nutrition function, and age). The Pima Indians Diabetes problem taken from Machine Learning Repository UCI: https:// github. com/ LamaH amadeh/ Pima-India ns-Diabe tes-DataS et-UCI/ blob/ master/ pima_ india ns_ diabe tes. txt

Conflict of interest
The authors declare that they have no conflict of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.