An enhanced version of Harris Hawks Optimization by dimension learning-based hunting for Breast Cancer Detection

Swarm intelligence techniques have a vast range of real world applications.Some applications are in the domain of medical data mining where, main attention is on structure models for the classification and expectation of numerous diseases. These biomedical applications have grabbed the interest of numerous researchers because these are most serious and prevalent causes of death among the human whole world out of which breast cancer is the most serious issue. Mammography is the initial screening assessment of breast cancer. In this study, an enhanced version of Harris Hawks Optimization (HHO) approach has been developed for biomedical databases, known as DLHO. This approach has been introduced by integrating the merits of dimension learning-based hunting (DLH) search strategy with HHO. The main objective of this study is to alleviate the lack of crowd diversity, premature convergence of the HHO and the imbalance amid the exploration and exploitation. DLH search strategy utilizes a dissimilar method to paradigm a neighborhood for each search member in which the neighboring information can be shared amid search agents. This strategy helps in maintaining the diversity and the balance amid global and local search. To evaluate the DLHO lot of experiments have been taken such as (i) the performance of optimizers have analysed by using 29-CEC -2017 test suites, (ii) to demonstrate the effectiveness of the DLHO it has been tested on different biomedical databases out of which we have used two different databases for Breast i.e. MIAS and second database has been taken from the University of California at Irvine (UCI) Machine Learning Repository.Also to test the robustness of the proposed method its been tested on two other databases of such as Balloon and Heart taken from the UCI Machine Learning Repository. All the results are in the favour of the proposed technique.

Breast Cancer is the leading cause of deaths of women in all over the world and it happens to more than 8% of women in their lifetime 1 . The combination of medical science and technology contributes a lot to human health and their quality of life.It is common among women while rare among men. Breast cancer commonly affects women more than 40 years of age however younger women can also be affected especially with genetic predisposition (a genetic characteristic that influences the development of an individual organism under the influence of environmental conditions). It arises from the breast tissues mostly from the ductal carcinoma (the inner lining of milk ducts) or less frequently from the lobular carcinoma (the lobules that supply milk to the ducts). The risk factors for breast cancer are age, genetics, obesity, family history or late pregnancy 2 .Due to the factors related to cost and professional experience, in the last two decades computer systems to support detection and diagnosis have been developed in order to assist experts in early detection of abnormalities in their initial stages. Despite the large number of researches on computer-aided systems, there is still a need for improved computerized methods 3 . Breast cancer in men does not occur very often, less than 1% of all breast cancers occur in men. For men, the risk of being diagnosed with breast cancer during lifetime is about 1 in 1,000. Men and women all have breast tissues. The various hormones in women's bodies exhilarate the breast tissue to grow into full breasts while men's bodies normally don't form much of the breast-stimulating hormones. As a result, their breast tissue usually stays flat and small.Mammography is the most widely used modality for detecting and characterizing breast cancer. It is a medical imaging that uses low-dose X-ray system to see inside the breasts. A mammography exam, called a mammogram, helps in the early detection and diagnosis of breast diseases. It has high sensitivity and specificity Radiologists visually search mammograms for specific abnormalities. Some of the important signs of breast abnormalities that radiologists look for are: • Calcifications that are tiny mineral deposits (calcium) scattered throughout the mammary gland, or occur in clusters. They appear as small bright spots on the mammogram. Concentration of microcalcifications in one place is also in favor of malignancy while the scattered calcifications are usually benign. • Masses defined as a space occupying lesion. These are areas that look abnormal and they can be many things, including cysts (non-cancerous, fluid-filled sacs) and non-cancerous solid tumors (such as fibroadenomas) , but may sometimes may be a sign of cancer. • Architectural distortions are when the normal architecture is distorted with no definite mass visible. This includes spiculations radiating from a point, and focal retraction or distortion of the edge of the parenchyma. It appears as a distortion in which surrounding breast tissues appear to be "pulled inward" into a focal point.
In the last few decades, many methods have been adopted for making such systems which can help in assisting radiologists in detecting any type of abnormalities in a mammogram for initial investigation. But due to large amount of data and their complexities no single algorithm is sufficient for that so we needed more enhanced versions for tackling these issues. Some works already done by the researchers have been described in this phase. Tariq Sadad et al. 5 proposed a CAD system in which segmentation technique i.e. FCMRG (Fuzzy C-mean Region Growing) is applied to obtain the mass in the image. Different features are extracted such as LBP-GLCM and LPQ (local Phase Quantization) and the feature selection is done on the basis of mRMR (minimum-redundancy maximum-relevancy) algorithm. At the end classification is carried out using different classifiers such as Decision Tree (DT), Support Vector Machine (SVM), Logistic Regression (LR), Linear Discriminant Analysis (LDA), K-nearest Neighbors (KNN) and Ensemble classifier to differentiate the benign tumors from the malignant ones achieving an accuracy of 98.2%. Elmoufidi et al. 6 proposed a method to segment and detect the boundary of different breast tissue regions in mammograms by using dynamic K-means clustering algorithm and Seed Based Region Growing (SBRG) techniques. Mammographic Image Analysis Society (MIAS) database is used for evaluation. Pratiwi et al 7 proposed comparison of two classification methods: Radial Basis propagation Neural based on Gray-level Co-occurrence Matrix (GLCM) texture based features and BPNN. The computational experiments are performed on MIAS database showing that RBFNN is better than BPNN in breast cancer classification. Jen et al. 8 proposed an efficient abnormality detection classifier (ADC) in mammogram images. Firstly, preprocessing is performed that included global equalization transformation, image denoising, binarization, breast  www.nature.com/scientificreports/ object extraction, determination of breast orientation and the pectoral muscle suppression. On the obtained segmented images, gray level quantization is performed and further five features are extracted from the ROI and PCA is applied for determining feature weights. Azevedo et al. 9 worked on IRMA database of mammograms that contains four types of tissues: fat, fibroid, dense and extremely dense. In this work, Morphological Extreme Learning machine is proposed with a hidden layer kernel based on morphological operators: dilation and erosion for classifying the masses as benign or malignant. Pereira et al. 3 presents abnormality detection method in CC and MLO view of mammograms. Preprocessing included artifact removal algorithm followed by an image denoising and enhancement based on wavelet transform and Wiener filter. For segmentation of masses three techniques are used: multiple thresholding, wavelet transforms and genetic algorithm.DDSM database is used for experimentation Area overlap metric (AOM) achieved by the proposed method is 79%. Chen et al. 10 proposed a clustering approach based on combination of fuzzy C-mean clustering (FCM) with PSO. As fuzzy C-mean clustering has some drawbacks such as the number of clusters needs to be specified in advance and also we need to have knowledge of the ground truth. The data points in overlapping areas cannot be correctly categorized so to overcome all these drawbacks PSO is used in combination with FCM. The result of this algorithm shows that it can automatically find the optimal number of clusters. Shanmugavadivu et al. 11 proposed an intuitive segmentation technique to separate the microcalcification regions from the mammogram. This mechanism first enhances the input image using an image-dependent threshold value and binarizes it to obtain an enhanced image. Then the pixels constituting the edges of microcalcification regions are grown in the enhanced image, with respect to the neighbourhood pixels. Lastly, the edge intensities of the enhanced image are remapped into the original image, using which the regions of interest are segmented. The results of the present work are compared with the ground realities of the sample images obtained from MIAS database. Rao et al. 12 proposed a new optimization method known as Teaching Learning-based optimization (TLBO) which is based on the influence of a teacher on the learners. It is a population-based method that uses a population of solution to reach to a global solution. Rangayyan et al. 13 presents an overview of various digital image processing and pattern analysis techniques for problem solution in several areas of breast cancer diagnosis. This includes: contrast enhancement, detection and analysis of calcifications, masses and tumors, asymmetric shapes analysis and detection of architectural distortion. Mirjalili 14 has utilized the GWO algorithm for training MLP for evaluating the accuracy and classification rate of three database issues. On the basis of experimental outcomes, proved that the GWO strategy is competent in giving the superior quality of outcomes in terms of enhanced local optima avoidance. After that the various authors have utilized the many recent enhanced and hybrid techniques on these database for effective solutions 15,16 .For locating the cancerous region in the mammogram image a comprehensive algorithm has been utilized by Sha 17 .This strategy has been applied on image noise reduction, optimal image segmentation for feature selection and extraction, thereby reducing the evaluating cost and enhancing precision etc. In addition, the proposed algorithm was utilized by the MIAS and DDSM databases. For verification, the performance of the presented approach has been validated with ten recent optimizers. Through tabulated outcomes it has been proved that the presented approach is competent to give the 96% Sensitivity, 93% Specificity, 85% PPV, 97% NPV, 92% accuracy, and better efficiency as compared to others. Shaikh et al. 18 has introduced a new hybrid version by merging the features of harmony search (HS) and simulated annealing (SA) for precise and accurate breast malignancy. In this work, has been utilized 02 breast databases such as (i) benchmark BCDR-F03 database and (ii) local mammographic database. On the basis of tabulated outcomes have been proven the robustness of the proposed approach for these databases. Houssein et al. 19 gave a method to overcome the problem of abundant data analysis in cheminformatics.This paper proposed a hybrid method named CHHO-CS in which harris hawk optimizer is combined with two other operators i.e. cuckoo search and chaotic maps. Controlling of position vectors of HHO algorithm is done by cuckoo search. For prohibiting the control energy parameters to fall in local optima, chaotic maps are used. SVM is used as objective function alongwith CHHO-CS for identifying best features and removing repetitive data. This proposed method is tested on various databases and gave very promisisng results over traditional methods like HHO, moth-flame optimization, grey wolf optimizer and others. Houssein et al. 20 proposed 3 modified versions named HHOCM i.e HHO with crossover and mutation CM, OBLHHOCM i.e opposite-based learning HHOCM and ROBLHHOCM i.e random opposition based learning HHOCM. These methods are tested on two datasets and results ensured best performance in finding subsets of molecular descriptors over other algorithms such as salp swarm optimization, original HHO, dragonfly algorithm, whale optimization algorithm in literature. Houssein et al. 21 used Harris hawk optimization algorithm in large scale sensor network to find the sink node which helps the whole network in processing and analysis of information. As finding a sink node is a challenging task, HHO performed very well as compared to other optimization algorithms such as particle swarm optimization, sine cosine algorithm, flower pollination algorithms and others in literature. HHO gave good results in terms of increasing network lifetime and energy consumption for different size of networks as well as for multiple sink nodes. Recently a newly hybrid SCHHO approach has been developed by Hussain et al. 22 for optimization problems and feature selection, which incorporates the features of sine-cosine algorithm (SCA) in Harris hawks optimization (HHO) algorithms for improving the exploration and exploitation performance. This approach is assessed through using recent test suites and sixteen datasets and compared with recent metaheuristics.
Nowadays the various hybrid and enhanced algorithms are developing by the researcher such as OHGBPPSO 23  www.nature.com/scientificreports/ the proposed algorithm with the most recent optimizer strategies. To sum up, the main contributions of this new work are: • A new enhanced version namely DLHO that includes features from DLH search strategy is proposed.
• To segment the cancerous region from mammogram, based on a fresh structure of the multilayer perceptron (MLP) neural network using DLHO algorithm. • Utilizing DLHO algorithm for disease detection in other biomedical databases also. • Statistical and qualitative experimental analyses to assess the robustness and effectiveness of the DLHO algorithm as compared to recent algorithms.

Harris Hawks optimizer algorithm (HHO)
Harris Hawks optimizer approach is recently developed by Heidari et al. 31 . It's a population based method and inspired through the intelligence of the crowd in which main inspiration is the co-operative behavior and chasing style of Harris' hawks in nature called surprise pounce. During the hunting of the target each search agent or hawks jointly pounce prey by different positions. This methodology is implemented for finding the best target in the complex search space.
Exploration stage. On the basis of the following equations, each search member randomly visits each location and waits to find a target; where X(t + 1) and t, X rabbit (t) are shows the the position vector of search agent and rabbit in the next generation, X(t) is show the recent location of search agents, r 1 , r 2 , r 3 , r 4 , and q are random numbers inside (0,1), which are updated in each iteration, LB and UB show the upper and lower bounds of variables, X rand (t) is a randomly selected search agent from the recent crowd, and X m is the average position of the recent crowd of search agents. The mean location of the search member is attained by the following Eq. (2): where X i (t) , t and N illustrates the position of each search agent, iteration and total number of search agents respectively.
Transition from exploration to exploitation. During this phase, energy of the prey is measured through the following Eq. 3; where E,T and E 0 illustrates the escaping energy of the target, total generations and initial state of its energy.
Exploitation phase. Soft besiege The behavior of the search agents is attained through the following conditions; where J = 2(1 − r 5 ) , �X(t) , r 5 and t are presents the random jump strength of the rabbit throughout the escaping procedure,the difference between the location vector of the rabbit, random number and recent position in generation.The J value changes randomly in all iterations to simulate the nature of prey or rabbit motions. Hard besiege In this phase, the present positions of all agents are updated by Eq.(6):

Soft besiege with progressive rapid dives
The next move of the search member has been evaluated through the following (Eq. 7): Each member are drive by the following rule; where D, 1 × D are presents the dimension of function, random vector and LF is the levy flight function, Which is calculated through the Eq. (9): (1) where Y and Z are obtained using Eqs. (7) and (8).

Hard besiege with progressive rapid dives
The following rule is performed in hard besiege condition: where Y and Z are obtained using new rules in Eqs. (12) and (13).
where X m (t) is obtained using Eq. (2).

Enhanced DLHO Algorithm
The complex optimization applications are challenges for the optimization of meta-heuristics. According to literature, each optimization method is not able to show the best solution for all types of complex problems. All algorithms may face some drawbacks, so due to these weaknesses these could fail to find the solution of complex functions.
Cancer related issues are a big challenge for the bio-medical field researchers. Due to their complexity each optimization method is not competent to tackle these issues. Therefore, we always need the most robust optimizer method for the future demand. After this inspiration, to address these issues, a technique named DLHO has been developed by the merits of the HHO and dimension learning-based hunting (DLH) search strategy.In DLHO,the exploitation and exploration phases of HHO has been enhanced by DLH method to demonstrate the merits of DLH search strategy.DLH search strategy utilizes a dissimilar method to paradigm a neighborhood for each search member in which the neighboring information can be shared amid search agents.This strategy helps in maintaining the diversity and the balance amid global and local search. The exclusive motivation for overdue mixing modifications in HHO is to advantage the process to evade immature convergence and to steer the search in the way of possible exploration or search area in a faster direction.
By this mechanism, a hawk outcome is competent to escape from the local optima outcome. Also the accuracy or quality of the outcome is extended with faster convergence speed. With that the search members have explored the large search areas for trapping the superior outcome or goal. This procedure is repeated again and again until a new outcome or position has not fulfilled the termination conditions.
Mathematical phases of DLHO algorithm are as follows: • Parameters The following parameter settings (see in Table 1) have been taken during the implementation of the algorithms. • Crowd initialization Firstly, the crowd in the exploration area has been initialized randomly. In this, optimizers allocate a random n d for the i th hawk; H i (i = 1, 2, . . . , n) . In the exploration area each slime has been allocated as: . . iter  500  500  500  500  500  500  500   Crowd size  30  30  30  30  30  30  30   lb  10  10  10  10  10  10  10   ub  100  100  100  100  100  100  100 No. of (runs) 30  www.nature.com/scientificreports/ Where n and d illustrate the hawks and size etc. The fitness of each slime has been amended by the equation (15); where fh i illustrates the fitness outcome of the ith hawk during the searching procedure in the exploration space. In addition the following two distinct matrix can be formulated for target as in Eqs. ( (16)- (17)); where ft i denotes the fitness outcome for the final target. • Fitness Evaluation The fitness value of the hawk is calculated through the following equations (18)-(19); where F best and F worst are illustrates the best and worst fitness value for functions.
• Drive stage During this phase, the DLH helps in specifically chasing the prey or goal of the search member in the search domain. DLH is used to construct a local area for each search member in which the best area or neighbors information can be shared between search members. Additionally, next phase represents how DLH and canonical HHO phases make two different search agents. • Canonical HHO phase In HHO method phase, the position of the all member in the search domain is updated by equations ((1)- (13)).Finally the first search member for the new location or position of HHO X i (t) named X HHO (t + 1) is updated by equation (11). • DLH phase In DLH, the position of the search agents is modified by an equation in which this different search agent is learned by its different neighbors or local optima and a randomly chosen search member from the population.After that, besides X i−HHO (t + 1) , the DLH phase creates other candidate for the new position of hawk X i (t) named X i−DLH (t + 1) . To do this, a radius R t (n) is evaluated applying Euclidean distance among present location of hawk X i (t) and the candidate position X i−HHO by equation (20).
Then, the near hawks of X i (t) indicated by L t (n) are calculated by equation (21) with respect to equation (20).Here M i illustrates the Euclidean distance amid X i and X i .
After that, the multi-neighbors learning is executed by the following mathematical formulation of equation (22); where d is illustrate the dimension. • Position update stage In this phase, the best hawk is nominated by comparing the fitness outputs of two hawks X i−HHO (t + 1) and X i−DLH (t + 1) by equation; www.nature.com/scientificreports/ So, as per above equation if the fitness output of the chosen member is < X i (t) , then X i (t) is modified by the chosen member. Otherwise, X i (t) remains unchanged in the crowd. Pseudocode of DLHO. The pseudocode of DLHO algorithm is reported in Algorithm 1.

Analysis and discussions
Analysis. Test suites and constants. All methods have been coded in R-Matlab-2018a software and tested on 8GB Ram with 64 for bit operating system and Core i3, 8th Gen system for evaluating the robustness performance of the algorithms. Under this implementation of these methods the various constant settings have been fixed such as search agents ( n = 30 ), size of problems (10-100), upper and lower size taken amind -100 to 100 respectively. The performance of the proposed method has been tested on 29 CEC functions, these shown in Table 2 39 . The three dimension graphs of CEC are illustrated by Fig. 2. Normally, these suites could be divided into four phases as uni-modal, simple multi-modal, hybrid and composition etc.
Testing and evaluation. To confirm the effectiveness, robustness and accuracy of the proposed optimizer it has been run on 29-CEC test suite. Numerical and statistical results of the proposed optimizer proves that the proposed method is able to give the highly effective and accurate solutions in terms of min and max objective scores, mean and standard scores etc. Results of the optimizers are denoted by Tables 3-6  www.nature.com/scientificreports/ Further, in the following subsections, the brief details of the results analysis and discussions have been described.
Discussions on the experimental results. Generally the test suites can be split into four phases and each phase is utilized to prove the robustness of the optimizer on different levels such as; • Uni-modal: These are utilized to present the exploitation capability of the optimizers.
• Multi-modal: These functions involves various local optima that are applied to examine the capability of optimizers for many optima. • Hybrid and Composite: These functions are used to assess the exploration capability of the optimizers.
Exploitation competence. Uni-modal functions involve a single global optima normally used for evaluating the exploitation performance of the optimizers. Results in Tables (3)(4)(5), illustrated that the DLHO is competent to give the best exploitation performance in the search space as compared to other optimizers. Experiments have proven that the DLHO method can handle the uni-modal functions easily and is able to give the best and accurate score on these functions than others. Here, it could be concluded that the enhancement of HHO method allowed the functions to reach the best global optima. So, the DLHO method could tackle the high domain and complex functions easily.
As specified previously, the CEC test suites are more suitable functions for evaluating or testing the robustness of the optimizers. All Simulation proves that the DLHO is extremely functional. Results give strong evidence that the DLHO is competent to give the most effective and accurate optima solutions for the complex domain functions.
Competence valuation. Multi-modal functions are more suitable for testing the suitability of the optimizer since these involve many local optima, and the number of variables exceed exponentially against the size of the Balance Competence valuation. In hybrid and composite functions various local and global optima are involved normally which are utilized for evaluating the exploration robustness of the optimizers. In addition, these suites are utilized for evaluating and verifying the balance amid exploration and exploitation phases of the optimizers. Results showed that the DLHO is competent to make a strong balance amid exploration and exploitation phases. So, by this enhancement we can find the most effective solutions for the complex functions which illustrate the www.nature.com/scientificreports/ strong exploitation and exploration behavior of the DLHO method. Furthermore, this modification helps to make the new position or other best fitness location for each search agent which helps to amend their present locations. Experimental measurements illustrate the effective detection behavior of the DLHO for trapping the best scores against the least number of generations. Since hybrid and composite suites involve more complex space which shows the robustness of the new method. Hence these experiments showed that the proposed method leads to best optima by the most effective scan behavior. Additionally, the DLHO involves robust local and global optima prevention which helps in trapping the best goal quickly.
Accuracy. In this phase the accuracy behavior of the proposed optimizer will be discussed. In general, the least average score shows the accuracy of the optimizer against the best score. The best average scores of the optimizer are illustrated in a Table 7 into two different phases such as best (B) and worst (W) respectively. The best and worst average scores have been assigned in Table 7 by the outcomes of the Table 5. Results of the table show that the DLHO is competent to trap the best score against the least average score than others. Hence it can be said that the new enhanced version is competent to trap the best and accurate outcomes for the complex issues as comparison to others.
Stability. In general, if the standard scores lie near to '0' then these illustrates the stability of the optimizers for the global optima in the search space. The standard scores of the optimizers have been plotted by the Fig. 3 against the Table 6. In Fig. 3, it can be easily seen that the DLHO method is finding the best global optima for each function against the least standard score than others. These results show better stability performance of the new version than others on all functions. Additionally the least sd score also illustrates the convergence speed of the optimizers. On the basis of the all outcomes of Table 6 and Fig. 3 it can be concluded that the DLHO is competent to finding the best outcomes for complex issues much more fast without losing their track. www.nature.com/scientificreports/ Convergence graphs analysis and discussion. The performance graphs of the methods have been plotted between maximum number of generations and best outcomes so far, these shown in Figs. 4-7. These figures show that the best outcomes are obtained beside each generation which are shown by the best outcomes. According to Berg et al. 40 , assures that the optimizers ultimately trace to a goal and catch the best outcomes in search areas. So, the DLHO expands the fitness outcome for every member and guarantees best goals for issues as generation increases. Now, it can be said that it happens by the enhancement of HHO method.Each member travels from highest to lowest optima and with this strategy the overall members are amends their locations with every generation.
With this strategy, the best outcomes are stored for searching the next best location for each member of the crowd in the search space. All these graphs show that the DLHO is competent in trapping the best outcomes for each function against the least number of generations which illustrates the fast convergence speed of the DLHO algorithm.
Assessment of robustness. The accuracy and robustness of the presented approach has been confirmed by utilizing the Wilcoxon signed ranks method for a superior assessment 41 .It is a non-parametric strategy that is applied on different samples for judging the significance difference amid them.The significance level shows the best sample amid them. However, this strategy generally helps in locating the significant difference of the behaviors of two algorithms.
The results of the Table 8 have been evaluated by the results of the Table 7. In the Table 8,if p < 1 , then it shows a rejection of the null hypothesis ( H 0 ), whereas p > 1 , shows a failure to reject the ( H 0 ).So, p < 0.05 , it shows that the proposed strategy is significantly better than the other algorithms.Other hand, if p > 0.05 , it presents that the attained enhancements are not statistically significant.
Here, in this phase for assessment the robustness of the present strategy the wilcoxon method has been utilized against the average values of the algorithms so that it could be concluded that the significant difference is amid the algorithms or not.In Table 8,it can be easily seen that the proposed strategy has better characteristics such as strength of the global optima goal and superiority of the optimal solution. Also, significant importance may be placed in local exploitation and global exploration.The Wilcoxon method results illustrated that the proposed www.nature.com/scientificreports/ strategy is superior among-others-in-comparison.Hence, the proposed strategy is-statistically superior and this has not happened by likelihood/or chance.

Proposed technique(DLHO) for breast cancer detection
In this work, firstly the mammography pictures are fed into an MLP neural network (NN) for pre-processing to categorize the nonlinear distinguishable structures of the pictures. Secondly, preprocessed output has been passed to DLHO.The structure of the proposed methodology has been illustrated in following subsections: Picture pre-processing. Occasionally, Mammograms have a slight sound due to fluctuation and accidental alterations in dignified motions. Sound is a serious issue for picture-processing developments, mostly when boundaries inside the pictures have to be recognized, which needs separation. Separation rises the consequence of extraordinary frequency pixels, which comprises sound. In addition the median filter is utilized for the preprocessing procedure earlier picture division. it recalls boundaries although eliminating sound. It is a lowpass filter that needs a lengthier dispensation time than other filters. The median filter exchanges the particular pixel by the median of the nearby pixels. It is evaluate by mathematical is as follow; Where M illustrate the median and δ shows the neighbourhood mid nearby the position (m, n) in an picture.

MLP neural networks (NNs).
The NNs are unique of the extreme developments in the domain of artificial and computational intelligence. They mimic the neurons of human brain to frequently resolve the big databases related to bio-medical science. Every link amid the neurons is allocated a unusual weight, which specifies the amount of influence of the production on the contribution to the subsequent neuron. Usually, a neuron also has its individual weight, mentioned to as a bias, which defines the influence of the neuron on itself 42 .The basic structure of the MLP Neural Networks (NNs) has been illustrated by Fig. 8.  www.nature.com/scientificreports/ The outcome of every node can be calculated by two phases (1) firstly, the sum of the weight is calculated by the following equation; Where x i , w ij and h j represents the input, weight and hidden neuron j.Secondly, the activation function (AF) is utilized to evaluate the outcome of the neurons. During this study, the sigmoid function has been utilized in AF and it can be evaluate by the following equation; At end, the outcome of all neuron can be calculated by the following formula; DLHO based MLP trainer. The multilayer perceptron (MLP) plays an important role in the field of biomedical science. Firstly, the optimizer method was utilized on these databases by. The weight and biases are the most important decision variables in these databases. Normally these variables are utilized for obtaining the superior prediction accuracy, approximation and classification respectively.Here the DLHO method accepts these decision variables in the form of a vector and it could be evaluate is as follows; . . , W n,n , h, θ 1 , θ 2 , . . . , θ h www.nature.com/scientificreports/ where n, θ j and W i,j illustrates the number of input nodes, connection weight and bias of the j th node. Further, the mean square error has been obtained by utilizing the following equation;  where s and m illustrates the number of training trails and number of outcomes respectively.At end, the fitness function has been obtained by the above equations is as follows; DLHO Algorithm and MSE. In this stage of implementation the proposed DLHO method has been applied to detect the breast cancer.Breast cancer detection is one of the major challenging well-known issue in the biomedical field. From last few decades, scientists of different fields are trying to solve this issue by various robust optimizers. Due to the complexities of this issue, a robust technique is needed so that we could tackle this problem. By merging the merits of two powerful techniques such as HHO and DHL 43 search strategy, we have tried to present a new method for tackling this issue so that the best solution can be generated.
Metric MSE(mean square error) has been utilized to evaluate the difference between desired and actual outcomes for each sample or object. MSE can be calculated by the following equation; where y j (k) and y * j (k) are denotes the actual and desire outcomes.The main objective of this study, the new hybrid version is utilization as a tool for breast cancer identification.
Post-processing. During this work, the kapur's strategy has been utilized for thresholding. Let L show the levels of a picture with m pixels. Then the average occurrence of a orange area i has been evaluated is as follows; where k(i) is illustrate the number of orange areas i in the picture.During the segmentation of the picture into various classes (m) we need m − 1 outcomes. So limit of the orange levels for all class concerning optimal thresholds is defined as; where z 0 = 0 and z m − 1 = L − 1.On the basis of above equation orange levels w k has been calculated by the following equation;  www.nature.com/scientificreports/ On the basis of above fitness function the optimum threshold has been obtained.
Database report. To verify the applicability of the proposed method for breast cancer detection, experiments have been performed on following data sets: • First database First MIAS database has been taken from the Pilot European Image Processing Archive (PEIPA) at the University of Essex 44 .It includes 322 digitized mammography pictures or images (Amid it contain 202 normal and 120 abnormal pictures) and size of each image is 1024 × 1024. The sample of these images are illustrated by Fig. 9. Figure 7. Performance graphs of Optimizers on composition test suites. www.nature.com/scientificreports/  Table 13.
Results and discussion. Statistical and numerical computations are done on the first database that included 322 digitized mammographic images out of which 202 are normal and 120 are abnormal pictures and the size of this database is taken as 1024 × 1024. This work has been implement on Matlab R2018a with a 64-bit operating system, Core i3, 8th GEN and 8GB Ram respectively. For verifying the robustness of the proposed method the obtained outcomes has been compared with recent algorithms such as BP 46 , GA 47 , PSO 48 , GWO 14 and WOA 49 respectively in terms of correct detection rate (CDR),  where D r ,C p and T p are illustrates the correct detection rate, number of pixels correctly classified and total number of pixels in the database respectively.
where F a ,H p and T p denote the false acceptance rate, number of healthy pixels classified as breast cancer and total number of pixels in the database.
where F r ,C p and T p denote the false rejection rate, number of cancer pixels classified as healthy and total number of pixels in the database. Experimental outcomes of the algorithms on the MIAS database have been illustrated through Table 9 and Fig. 10.The graph is plotted (Fig. (10)) to compare the cancer detection rate (%) of various algorithms.
During the implementation of the proposed algorithm on the second database for breast cancer,various input parameters are set such as number of attributes (9), number of training objects (599), number of test objects (100) and number of classes (2) respectively. The algorithm is run iteratively to find the best optima.Tabulated outcomes of the methods are illustrated in the Table 11. The classification rate of the methods against this database are drawn by Fig. 13. Figure (12) depicts the detection of left and right breast cancer detection in eight mammogram pictures using the DLHO algorithm.
Again the values in the Tables 9-11, show the superiority of the proposed method.The statistical score shows that DLHO method has a high ability to ignore local optima and approximate the superior global optima outcomes for biases and weights.This issue has the utmost complexity compared to the earlier discussed database Here noticeably the cancer detection rate (%) of the maximum optimizer on this issue is very poor due the complexity of this data set.However, all the methods have been tested on the same parameter settings so that we can correctly find out the robustness of the algorithms.Tabulated values of the methods strongly shows that the proposed method is competent to solve this problem with best optima, accuracy and high rate of cancer detection as compared to others.The computational time of the algorithms for cancer detection have been reported in the Tables 10-12. All the values in the tables have proved that the proposed method is competent to detect the cancer with highest detection rate as well as have least computational time as compared to others. Figures 11-14, are shows the mean square error of the algorithm on the first and second database. As stated, the metaheuristics were run various times to confirm the robustness. These graphs show that the proposed method successfully solves these databases with least mean square errors.
By experiments, it could be determined that a huge drawback of the recent algorithms is that numerous weaknesses such as getting stuck at a local optima, slow convergence, weak exploration and exploitation balance and premature convergence cause them to fail to resolve complex issues in certain cases (see .

Robustness of DLHO for other biomedical databases
In this part of research, two different bio-medical issues such as balloon and heart are taken from the University of California at Irvine (UCI) Machine Learning Repository 45 . On these databases various population based optimizers are applied 14,50-52 for evaluating the high rate of classification and accuracy of the outcomes.The weight and biases range has been fixed -10 to 10 for each database. The crowd size is set 50 for Balloon, 200 for the rest databases and 250 maximum number of generations have taken during the implementation of the codes.
The input values for these databases have been illustrated by Table13. Each method has been run various times to generate the best outcomes for each databases. The outcomes of the optimizers have been presented in the terms of average, standard and classification rate (%) respectively. Noticeably, the least mean and standard score of mean square error in the end of the generation shows superior performance and accuracy.
The normalization used in this work is called min-max Normalization which is the most important stage of MLP when we solve these databases with attributes in different sizes or ranges. During this work, the stage is known as min and max normalization and is evaluate is as follows; where the above mathematical equation maps x in the range of [l, m]- [n, o]. Figure 11. MSE graphs of algorithms on first database for cancer detection. On the basis of experimental outcomes, firstly we can seen that the classification rate of each algorithm is same, however the average and standard scores of the algorithms for this database are unique. The least average and standard score shows the accuracy, robustness and stability of the algorithms for the functions. So on the www.nature.com/scientificreports/   www.nature.com/scientificreports/   www.nature.com/scientificreports/ basis of these outcomes we can concluded that the DLHO algorithm is competent to provide the high accurate, robustness and stability performance of this issues than others. The statistical outcomes proved the robustness of the proposed method.
Heart. In lastly, the heart database has been solved by the optimizer methods which have number of attributes (22), number of training objects (80), number of test objects (187) and number of classes (2). Multi-layer perceptron's against the construction of 22-45-1 has been trained through the metaheuristics. Obtained outcomes of the metaheuristics have been reported by table 15 and the classification rate of the algorithms on the heart database are plotted by Fig. 15. Experimental outcomes revealed that the DLHO is competent at giving the superior classification rate with the best statistical outcomes. Here, it could be concluded that the proposed method most effective in issue approximation heart database. The least mean and standard deviation outcomes reveals the better local optima avoidance of the proposed method.All simulation also illustrates that the superior error fits the proposed method. This validates the accuracy and best performance of the DLHO method as well.

Conclusion
This work points out an enhanced version with the integration of merits of Harris Hawks Optimization (HHO) algorithm and dimension learning based hunting (DLH) search strategy for biomedical issues, it is called DLHO algorithm. In this modification the DLH search strategy has been utilized for enhancing the exploration and exploitation stages of the HHO algorithm, so that the above weakness could be removed. For verification of effectiveness the DLHO method has utilized 29 standard suites and three database. However, the robustness of the classifiers has been verified on all features and particular attributes distinctly to get and compare the attained accuracy. All simulations reveal that the presented algorithm is able i) to detect the high quality of global solutions and ii) to define the effective optima outcomes for biases and weight in terms of local avoidance and detection or classification rate or produces high accuracy for the biomedical issues as compared to others.
For future work, we shall present a more enhanced version for biomedical and engineering application. At end, we believe this research will inspire every young scientist, who is recently working on meta-heuristics and engineering applications.