key: cord-0942190-7yzs8c08
authors: Shadkam, Elham
title: Parameter setting of meta-heuristic algorithms: a new hybrid method based on DEA and RSM
date: 2021-11-17
journal: Environ Sci Pollut Res Int
DOI: 10.1007/s11356-021-17364-y
sha: 5b106ec0e830ae97542f922318049d375a858ccb
doc_id: 942190
cord_uid: 7yzs8c08

The parameter setting of meta-heuristic algorithms is one of the most effective issues in the performance of meta-heuristic algorithms and is usually done experimentally which is very time-consuming. In this research, a new hybrid method for selecting the optimal parameters of meta-heuristic algorithms is presented. The proposed method is a combination of data envelopment analysis method and response surface methodology, called DSM. In addition to optimizing parameters, it also simultaneously maximizes efficiency. In this research, the hybrid DSM method has been used to set the parameters of the cuckoo optimization algorithm to optimize the standard and experimental functions of Ackley and Rastrigin. In addition to standard functions, in order to evaluate the performance of the proposed method in real problems, the parameter of reverse logistics problem for COVID-19 waste management has been adjusted using the DSM method, and the results show better performance of the DSM method in terms of solution time, number of iterations, efficiency, and accuracy of the objective function compared to other.

Many real-world problems can be modeled by operation research techniques and then optimized in a variety of ways. The goal of optimization is to find the best acceptable solution, given the constraints and assumptions of the problem. In each optimization problem, there are decision variables that are limited by a set of constraints or conditions. If the decision variables are discrete in nature, such problems are called combinatorial optimization problems. Because many potential solutions must be considered to solve combinatorial optimization problems, these problems are known as NP-hard problems. In many real-world problems, the solution time is increased non-linearly using exact optimization methods. To solve this problem, regardless of the accuracy of the solution, it is possible to create an acceptable solution in a much shorter and more logical time by using approximate optimization methods. Therefore, optimization methods and algorithms are divided into two categories: exact and approximate methods. Exact methods are able to find the optimal solution accurately but are not efficient enough in NP-hard optimization problems. Thus, heuristic algorithms were first developed, which are a set of approximate algorithms, but these algorithms are generally designed for a specific problem and often cannot be used for other problems. Therefore, meta-heuristic algorithms were proposed that are independent of a specific problem. The performance of these algorithms is highly dependent on their parameters. In order for these algorithms to work optimally, it is necessary to find specific values for each parameter. Therefore, the problem of parameter setting of meta-heuristic algorithms has been raised. Adjusting the parameters of meta-heuristic algorithms is an effective factor in the performance of these algorithms, so this problem has been considered in many previous studies and various methods have been proposed for it (Lobo et al., 2007; Yazdani, et al., 2017a Yazdani, et al., , 2017b . The importance of the parameters of any optimization algorithm, especially meta-heuristic algorithms, is inevitable . The optimal values of these parameters, which generally depend on the characteristics of the parameter and the problem, have a significant impact on the performance of these algorithms and better search of the solution space (Najafi & Behnoud, 2015; Yazdani, et al., 2017a Yazdani, et al., , 2017b . One of the recently introduced meta-heuristic algorithms is the cuckoo optimization algorithm (COA). The use of the COA has many applications in various engineering sciences. This algorithm, like other meta-heuristic algorithms, is sensitive to parameter setting. In using the COA, precise adjustments must be made in the execution of the algorithm, determination of the appropriate combination of parameters, final execution of the algorithm, and validation of the solutions. Previous research has used methods such as trial and error, experimental design, and Taguchi method to adjust the parameters of this algorithm (Azadeh et al., 2013) . The problem of setting the parameters of meta-heuristic algorithms is so important and noteworthy that an article entitled Why Tuning the Control Parameters of Meta-heuristic Algorithms Is So Important for Fair Comparison? has been proposed in this regard (Kazikova et al., 2020) . In a study, Cheng and Chang (2007) proposed a genetic algorithm to optimize the scheduling problem and proposed a hybrid method for optimizing the genetic algorithm parameters using the design of Taguchi experiments. They tested and optimized seven algorithm parameters including initial solutions, selection method, intersection approach, mutation rate, initial population size, algorithm velocity, and mutation method. Naderi et al. (2009) used the Taguchi method to adjust the parameters of the simulated annealing (SA) algorithm. They considered parameters such as initial solutions, coding pattern, and local search structure as controllable factors of SA and different levels for each factor. They then considered and optimized these factors using the Taguchi method. In a study, Bartz-Beielstein and Markon (2004) proposed a method for parameters setting of complex real-world optimization problems. This method is a combination of statistical design of experiment methods, regression analysis, design and analysis of computer experiments, and regression tree that has been used to set the parameters of the SA and evolutionary algorithm. In a paper, Ridge and Kudenko (2007) set the parameter of the ant colony optimization algorithm by designing experiments to solve the traveling salesman problem, which used the response surface methodology to evaluate the quality and time of the solutions. In their research, Fallahi et al. (2014) examined the parameters of the ant colony algorithm for traveling salesman problem, which were performed in five stages: parameter profile screening, clustering, parameter screening, response procedure methodology, and finally optimization. In a paper, solving the problem of vehicle routing using the Memetic algorithm, Saremi et al. (2007) considered setting the desired parameters using experimental design method and analysis of variance to investigate the interactions. Xu et al. (1998) examined the parameter setting in the Tabu search algorithm using statistical tests and experimental design. In this study, the decision tree minimization problem is considered by explaining three types of neighbor structure selection methods. This method has been used by other researchers to adjust the parameters of the Tabu search algorithm to solve the problem of vehicle routing. Beielstein et al. (2002) examined the parameters of the particle swarm optimization algorithm using experimental design techniques and analyzed their important parameters and interactions. In an article, Ridge (2007) adjusts the parameters of one type of the ant colony algorithms to solve the traveling salesman problem, and the results show the better performance of the mentioned algorithm with standardized parameters. In a research, Šilc et al. (2015) proposed a data mining approach to set the performance of meta-heuristic algorithms. In this paper, the idea of learning, a model for algorithm behavior using data mining is presented as the results obtained in parameter setting. A paper presented by Veček et al. (2016) includes a method called the chess ranking system. This approach has been used to compare and score evolutionary algorithms. In that method, the desirability and performance of a combination of parameters are obtained by comparing with other structures based on their score, deviation from the score, and the distance of the score. During the performance and in each iteration, the bad structures are significantly eliminated and new structures are obtained through mutations and combinations of existing structures. The results of this method show that there is not always a significant difference in the assumption of zero statistical tests, but the method mentioned in the paper shows better performance when using similar parametric structures. Also, parameter setting with this method has been faster than other methods such as the competitive method, and there is no limit in setting absolute parameters. The researches done in this field are very diverse and extensive (Table 1) , and the studies done in this field are reviewed.

One of the methods used for parameter setting of metaheuristic algorithms is response level methodology (RSM). RSM is a combination of mathematical and statistical techniques that is widely used for modeling and data analysis. In this method, response variables (outputs) are influenced by several independent variables (inputs), and the goal is to obtain the best values of inputs in order to optimize outputs (Hinkelmann & Kempthorne, 1994) . One of the first researches to determine the parameters of meta-heuristic algorithms using RSM is Wang and Wu's (1999) research, in which a six-step RSM-based process for identifying and optimizing the SA algorithm with computational time constraints has been performed. In a paper, Najafi and Behnoud (2015) presented a combined algorithm of RSM and data envelopment analysis (DEA) to adjust the parameters of meta-heuristic algorithms. Using the RSM, they set the parameters of a genetic algorithm to optimize the machine scheduling problem assuming sequence-dependent preparation times. In their method, in order to optimize the parameters of the algorithm, only the optimization of the objective function has been considered and the solution time of the algorithm has not been addressed. The use of the DEA method in optimization problems is very common and can be referred to Shadkam (2021a Shadkam ( , 2021b Shadkam ( , 2021c ) that in a twophase algorithm using DEA to optimize the objective function and maximize the efficiency of the problem simultaneously. Also, Shadkam and Bijari (2015) in a research have used the combination of RSM and DEA in order to optimize multi-objective problems for bank branches. For the first time in a review paper, Phan et al. (2020) have reviewed and classified the various methods of parameter setting of evolutionary algorithms. These methods are divided into two groups: parameter tuning and parameter control methods. Parameter tuning techniques are used to select offline parameters before the optimization algorithm is executed, but online parameter control techniques are used in dynamic optimization problems during an algorithm execution or after each algorithm execution. According to this classification, one of the most widely used and popular methods in the field is the experimental design of experiments method (DOE). In this method, the parameters are adjusted to create the best solution for the objective function. The Taguchi method is mostly used in this field. In fact, because only one solution can be examined, the value of the objective function is usually considered, and if we want to consider other response such as the number of iterations and the execution time of the algorithm, this process must be repeated Joshi and Bansal (2020) Ant colony optimization algorithm Parameter specification screening, clustering, parameter screening, RSM, and optimization 7 Saremi et al. (2007) Memetic algorithm Design of experiments and variance analysis 8 Šilc et al. (2015) Meta-heuristic algorithm Data mining 9

Veček et al. (2016) Evolutionary algorithms Comparison and scoring 10 Najafi and Behnoud (2015) Genetic algorithm RSM and DEA 11 Saeheaw (2020) HCSCROCFO-3Opt algorithm A random manner 12 Odili and Fatokun (2020) African buffalo optimization algorithm Mathematical modeling 13 Gomes and de Almeida (2020) Sunflower optimization Statistical method of mixture design 14 Joshi and Bansal (2020) Gravitational search algorithm Topological characteristics of the given optimization problem 15 Do Ngoc Luu and Anh (2021) Deep belief networks to predict time series data Harmony search algorithm 16 Özakın and Kaya (2020) Air-based PVT system Taguchi method and ANOVA 17 Walker and Craven (2020) Evolutionary multi-and many-objective optimization Multi-objective evolutionary algorithms A novel MOCRS-tuning method 21 Cheng et al. (2021) Genetic algorithm and particle swarm optimization Auto-tuning symbiotic organisms search algorithm 22 Alavi et al. (2021) A variable neighborhood search meta-heuristic method Conventional sequential optimization method 23 Devarapalli and Bhattacharyya (2021) Power system stabilizer Sine-cosine algorithm 24 Mergos and Yang (2021) Flower pollination algorithm A simple non-iterative, single-stage sampling tuning method 25 Tien Bui et al. (2021) Neural computing Whale optimization algorithm separately for each response. Finally, the parameters created will not perform well at all responses because each response is optimized separately. In our study, by combining DEA and RSM methods, a hybrid method called DSM is presented. The name of the proposed hybrid method is taken from the methods used in it. This method has various applications that will be discussed in future research by the authors of this paper. In principle, this method is presented to optimize multi-objective problems and another application is to parameterize the meta-heuristic algorithm, which has been used in this paper to optimize the parameters of the COA. The COA is one of the newest and most powerful meta-heuristic algorithms that has the ability to solve discrete and continuous problems. The parameter setting process of this algorithm is performed by the proposed DSM method, with the help of data obtained from the design of experiments and the implementation of the algorithm for both standard and real-time functions. In the hybrid DSM method, first, a number of decision-making units (DMUs) are selected, and the values of their inputs and outputs are determined. Then, the efficiency values of the DMUs are calculated using the DEA method. Then, with the RSM method, the response surfaces are made, which includes an efficiency surface and the response surfaces to the number of output variables. The efficiency surface is optimized, and the values of the optimal independent variables are obtained, then according to the optimal independent variables, the optimal response variables are obtained. The obtained independent variables and the optimal response variables are considered as a suitable configuration for the COA. Similarly, two tools, DEA and RSM, have been used for optimizing multi-surface problems in power plants of IRAN (Shadkam, 2021a (Shadkam, , 2021b (Shadkam, , 2021c . The proposed method of the paper overcome the shortcoming of DOE method and Najafi and Behnoud (2015) algorithm. In this method, in order to set the parameters, in addition to the value of the objective function, other responses (execution time and number of iterations) are also examined. Also, all three responses are optimized simultaneously, which ultimately leads to a better configuration for the algorithm.

In order to implement the proposed method in the real world, the medical waste management problem is examined. Medical wastes are one of the most important hazardous wastes and carry a large number of pathogenic bacteria (Sangkham, 2020) . This issue has special importance and place in research, and many articles on the design of medical wastes logistics networks have been presented (Kargar et al., 2020b) . The problem of waste disposal in hospitals is a very important problem that has been the subject of much research, including Wang et al. (2020) , Fletcher et al. (2021) , Adu et al. (2020) , and Homayouni and Pishvaee (2020). Due to the global outbreak of the coronavirus, the volume of medical wastes and the dangers posed by it have increased dramatically (Ozkaya et al., 2020) . In this regard, Shadkam (2021a Shadkam ( , 2021b Shadkam ( , 2021c has presented a mixed integer model for the management of COVID-19 wastes, especially its vaccine. The purpose of this model is to design a reverse supply chain model for safe and optimal disposal of medical wastes (focusing on corona vaccine wastes). Due to the containers and equipment used in COVID-19 vaccine packaging for easy transportation and maintenance, this model can be a great help in returning these containers to the factory and recycling them, saving you money and time. Therefore, the proposed model directs non-recyclable wastes to disposal and returns recyclable wastes such as COVID-19 vaccine containers to the plant after disinfection. Many researches have been done in the field of medical waste management, an example of which can be seen in Table 2 .

In this paper, the parameter setting of inverse logistics problem of corona waste management has been investigated. Today, reverse logistics is one of the main activities of supply chain management that covers all physical activities related to returned products (such as collection, recovery, recycling, and destruction) (Julianelli et al., 2020) . In this regard, proper design and implementation of reverse logistics, in addition to increasing customer satisfaction, reduce inventory and shipping costs. In this paper, the parameter setting of a mixed integer linear programming model for the design of an integrated direct and reverse logistics network is investigated. Due to the prevalence of coronavirus (COVID-19) and the increase in medical waste, the need for a reverse logistics system for waste management is strongly felt. For this purpose, the inverse logistics model in the field of coronavirus waste management, especially vaccine waste, has been investigated. This model was introduced by Shadkam (2021a Shadkam ( , 2021b Shadkam ( , 2021c . This model is solved using the cuckoo optimization algorithm, and before solving using the proposed DSM method, the cuckoo algorithm parameter will be adjusted for the inverse logistics problem. In the following research, in the next section, the COA is described, then the proposed DSM method is presented, and how to implement it to set the parameters of meta-heuristic algorithms will be examined. In the next section, the parameters of the COA are adjusted with the help of the Ackley and Rastrigin functions. The verification of the proposed method is examined through four approaches, and the results will be reviewed and interpreted. These approaches include the Taguchi method, RSM, and Najafi and Behnoud (2015) method. In the next section, the cuckoo algorithm parameter is adjusted to solve the inverse logistics problem of COVID-19 waste and is compared with the experimental approach of parameter setting.

This algorithm was proposed by Rajabioun (2011) and is inspired by the life of a cuckoo bird. In various papers, the performance of this algorithm has been compared with other similar algorithms, which can be referred to (Kahramanli, 2012) and (Shadkam & Bijari, 2014) for more information in this regard. This algorithm has been used in various fields, such as production planning (Akbarzadeh & Shadkam, 2015) , portfolio selection , supplier selection (Shadkam & Bijari, 2017) , and solving multi-objective optimization problems (Gorjestani et al., 2015; Shadkam & Jahani, 2015; and Borhanifar & Shadkam, 2016) , supply chain problem , and resource leveling in control project management (Shadkam, 2021a (Shadkam, , 2021b (Shadkam, , 2021c . Due to the good performance of this algorithm, improved versions of this algorithm have also been presented (Shadkam & Bijari, 2020) . The flowchart of this algorithm is presented in Fig. 1 . In the COA algorithm, similar to other evolutionary algorithms, it starts with a random population of initial cuckoos. Each cuckoo has positions called a habitat (relation 1), which indicates the array of the problem-solving variables.

Each cuckoo lays a number of eggs according to its assigned egg-laying radius (relation 2). Some of these eggs are destroyed, and some of the eggs remain and become adult cuckoos. Adult cuckoos look for places where they have a better chance of survival or where the average value of the objective function is optimal in those areas. After identifying the optimal location, each cuckoo migrates with some deviation from the optimal location and starts laying eggs. Due to dietary restrictions, only a few cuckoos can survive. Therefore, cuckoos with inappropriate values of the objective function are removed, and only the more efficient cuckoos will be transferred to the next iteration of the algorithm. This process is repeated in the same order as mentioned to establish the stop condition.

(1) (Shadkam, 2021a (Shadkam, , 2021b (Shadkam, , 2021c No Authors Multi-objective Multi-period Multi-product Uncertainty Case study Approach Software 

In most real-world problems, we have more than one response or output and optimizing multiple responses at the same time is one of the major problems in problems. Through design of experiments of data, the response surfaces is obtained, which is an equation between a response variable and independent (input) variables. In such cases, it is difficult to optimize several response surfaces simultaneously in order to obtain the best independent variables. Various methods have been proposed to solve such problems, which include constraining the problem using other objectives, the utility function, and the least-cost response surface. Since the proposed method includes two methods DEA and RSM, it is called DSM method. This proposed algorithm combines the RSM as one of the optimization methods. In order to make the optimization process more efficient, in addition to the RSM, the DEA method has been used and a hybrid model of both is presented. The main advantage of the DSM algorithm is to create an efficiency surface instead of generating multiple surfaces for each of the output or objective functions. In DSM algorithm, first, a suitable experimental design is selected to generate the experiments, then the simulation process is performed for the designed experiments and the inputs and outputs generated from the simulation are analyzed using DEA and RSM tools. Finally, using the surface obtained for efficiency, the input parameters of the system can be determined in such a way that the system under review has the highest value of efficiency, and according to these input values, the corresponding outputs can be output from the response surfaces obtained for each output. The main advantage of the proposed model is its applicability to an unlimited number of objective functions and simultaneous attention to efficiency maximization along with optimization of problem outputs. The DSM method steps Also, the pseudo-code of this method is as follows:

1. Design of experiment or data collection: Control factors and response variables of problem are determined. 2. Data normalization: To reduce the effect of different scales on the data, the collected data have been normalized. 3. Determination of efficiency for each experiment by DEA: Using the input and output values of step 2, the efficiency value of each experiment is calculated. 4. Making efficiency surface by RSM: Using the input values of step 2 and the efficiency values obtained from step 3, an efficient surface is creased using the RSM. 5. Calculate optimal inputs: The optimal values of input variables or control factors are obtained. 6. Calculate optimal outputs: The optimal values of outputs variables or responses are obtained.

Input and output variables are determined by previous research and expert opinion. This data is collected for a number of specific decision-making units. Normalizations are performed to reduce the effect of different scales on the data. In this study, the Euclidean norm mentioned in Eq. 3 is used.

Measuring the efficiency in the field of performance evaluation in companies and organizations is very important. DEA is a non-parametric method in performance evaluation. There are several models in the field of DEA. The CCR method in model 4 is a basic model that is used in this study to evaluate the efficiency of decision-making units.

The purpose of using the RSM (Eq. 5) is to determine the best variables (inputs) that can create the best values for the objective function (system output). The main advantage of the proposed method is in this section. In fact, instead of generating multiple response surfaces for the outputs separately, an efficiency surface according to Eq. 5 is presented to calculate the optimal level of parameters.

In fact, to construct the efficiency surface, we use normalized inputs as independent variables, as well as the functions calculated as response variables. After setting the efficiency surface by the Minitab software, the surface is optimized by the Lingo software, and we introduce the resulting solutions as optimal system inputs. It then sets a separate surface for each output, except that here each output is placed as a response variable. Then, the optimal independent variables are placed in the response surfaces, and the optimal response variables are obtained which are considered as the optimal outputs of the system.

In this section, in order to investigate the proposed DSM method, the optimal parameters of the COA are set. The details of the implementation of the proposed approach for setting the algorithm parameter are shown in Fig. 2 .

The Ackley and Rastrigin functions are common non-convex functions used as a performance test problem for optimization algorithms. These functions are typical examples of non-linear multimodal functions. The Ackley function was proposed by David Ackley in his 1987 PhD Dissertation (Ackley 1987). Also, the Rastrigin function was presented in 1974 by Rastrigin as a 2-dimensional function (Rastrigin, 1974) and has been expended by Rudolph (Rudolph, 1990) . The Ackley and Rastrigin functions are expressed as relations 6 and 7, respectively. Also, the diagrams of the Ackley and Rastrigin functions are shown in Fig. 3 . As can be seen, these functions have many local minimum, and therefore, it is difficult to find the global minimum. However, the cuckoo algorithm can find as well the optimal solution for these two functions. Also, n indicates the number of variables or dimensions of the problem, which is equal to 2 in Fig. 3 , but we consider 100 to implement the DSM method and set the parameter.

Step 1: Determine the problem variables

In this step, the independent variables (inputs) of the COA, which are the parameters of the algorithm, are determined, which are as follows:

1. The number of initial cuckoos ( The response variables (outputs) are as follows:

1. The value of the problem objective function 2. Execution time to solve the problem Step 2: Design of experiments By specifying the upper and lower bounds for each of the input variables (these bounds are written in front of the input parameters in parentheses), the design of experiments is performed using the Minitab software, the output of which is 32 test designs. In fact, each experiment is a configuration for the parameters of the meta-heuristic algorithm and consider as a DMU. Then, by placing each of these configurations of parameters in the COA, the optimization of Ackley and Rastrigin functions is performed with a number of variables equal to 100, and the values of the output variables determined in the previous step are measured, and the results for Ackley and Rastrigin functions are shown in Table 3 .

According to previous statements, the data are normalized with Euclidean norm. It should be noted that the type of the outputs values are cost, and Eq. 6 is used to normalize them (Eq. 8).

Step 4: Calculate efficiency by DEA Each of the experiments or configurations of the previous step is considered as a DMU, and the efficiency values for each algorithm configuration in both functions are calculated using the DEA method and Lingo software according to last columns of Table 3 . The calculation of efficiency values is done with the CCR approach.

(8)

Step 5: Efficiency surface and response surfaces by RSM model In this step, using the information obtained from steps 4 and 5, the surfaces for efficiency and output variables are constructed. Surfaces are created using the Minitab software, inputs are considered as independent variables, and efficiency calculated in step 5 is considered as a response or dependent variable. The efficiency surface for the Ackley and Rastrigin functions is in accordance with Eqs. 9 and 10.

Also, for the three mentioned outputs, separate procedures are considered in which the inputs are considered as independent variables and the outputs as the response variable. The Ackley function response surfaces are in the form of relations 11, 12, and 13. Also, the response functions of the Rastrigin function are in the form of 14, 15, and 16 relations.

The objective function surface for the Ackley function (y 1 ) is as follows: The time execution surface for the Ackley function (y 2 ):

The number of iteration surface for the Ackley function (y 3 ):

The objective function surface for the Rastrigin function (y 1 ):

The time execution surface for the Rastrigin function (y 2 ):

(12) y 1 = 3.19427 − 9.64494E − 04x 1 − 0.248918x 2 − 0.260315x 3 + 0.00982920x 4 − 0.00220977x 5 + 0.000213333x 1 x 2 − 2.21111E − 04x 1 x 3 + 0.000149136x 1 x 4 + 1.35556E − 06x 1 x 5 + 0.0283750x 2 x 3 − 0.00336111x 2 x 4 + 0.000143906x 2 x 5 + 0.00816667x 3 x 4 + 0.000201250x 3 x 5 + 1.07639E − 05x 4 x 5 (13) y 3 = 11.5064 + 0.0201019x 1 + 0.512674x 2 + 0.559201x 3 − 0.211343x 4

The number of iteration surface for the Rastrigin function (y 3 ):

Step 6: Optimize efficiency surfaces and response surfaces To find the optimal parameters of the COA, the surfaces obtained for efficiency (Eqs. 7 and 128 are optimized with Lingo software and the values of the independent variables are determined in Table 4 (Input) .

Then, to calculate the values of the corresponding response variables, the obtained values of the independent (15) Y 2 = 4.26218 + 0.000353820x 1 − 0.439723x 2 − 0.426627x 3 + 0.00826078x 4 − 0.00226608x 5 − 2.55500E − 04x 1 x 2 + 0.000104389x 1 x 3 + 8.09630E − 05x 1 x 4 + 4.83472E − 07x 1 x 5 + 0.0583781x 2 x 3 − 0.00410972x 2 x 4 + 0.000151867x 2 x 5 + 0.000165278x 3 x 4 + 0.000198133x 3 x 5 + 1.38924E − 05x 4 x 5 (16) 

In this section, the proposed DSM method is compared with similar methods through different approaches including numerical results, Najafi and Behnoud (2015) method, Taguchi method and general comparison of all methods.

In this step, using the optimal configuration obtained from the proposed algorithm of the paper, the COA is executed 10 times with the help of the MATLAB software. This process is performed separately for both the Ackley and Rastrigin functions. It should be noted that because the independent variables (inputs) must be integers, the previously obtained optimal variables are rounded out according to Table 4 . After executing the COA, the average of each response variable is calculated in Table 5 . This process is to validate the results of the proposed DSM method with the results of the actual implementation of the algorithm. In fact, the results of the analytical approach are compared with numerical results. To facilitate the comparison of the results of both methods, the absolute value of the difference of the response variables in both approaches is calculated according to last column of Table 5 . As can be seen, the difference between the optimal values obtained from the DSM method and the values obtained from the implementation of the algorithm is very small, and this indicates the optimality of the proposed DSM method.

In this section, the results of the proposed method are compared with the method of Najafi and Behnoud (2015) that were observed in the literature review. In this method, first, a central composite design was created, and then using the data related to this design, the complete second-order model was fitted. Then the optimal values of the parameters are determined by solving the non-linear model. Then, the efficiency of the problems solved by the algorithm is measured by considering the obtained parameters as inputs of the DEA method and the most efficient configuration of the genetic algorithm is determined. In this section, in order to validate the proposed DSM method, Najafi and Behnoud's (2015) method are implemented for the Ackley and Rastrigin functions. As can be seen from the process of this method, the algorithm parameter setting dose only pays attention to the response value of the objective function, and other important outputs such as execution time and number of iterations are not considered. The process of the Najafi and Behnoud (2015) methods are as follows.

Step 1: Determine the problem variables Input variables (parameters affecting the performance of the algorithm) are identified, which are the same as input variables of the DSM method. Also, the output is the value of the objective function of the problem. Unlike the proposed DSM method, only one response is considered that is the value of the objective function of the problem, and the execution time and number of iterations are ignored. Step 2: Calculate the values of the variables

The values of the inputs and the corresponding value of the objective function are determined for the COA similar to the DSM method (Table 3) .

Step

In order to implement the RSM method for the step 2 data, the response surface is created according to the values of the inputs for the output of the value of the objective function. Similarly, the response surface is the value of the objective function for the Ackley and Rastrigin functions according to Eqs. 9 and 12.

Step 4: Optimization of objective function surface

In order to determine the input parameters, the response surface of the value of the objective function (Eqs. 9 and 12) is considered, and the constraints related to the boundaries of the parameters of meta-heuristic are considered. Finally, the mathematical model of parameter setting related to the Ackley and Rastrigin functions is 17 and 18, respectively. These models are optimized with the Lingo software, and the optimal values obtained for both functions are shown in Table 5 (input).

y 1 = 2.87342 + 0.000230282x 1 + 0.00919271x 2 + 0.00825863x 3 + 0.00982920x 4

5 ≤ x 1 ≤ 20 2 ≤ x 2 ≤ 4 3 ≤ x 3 ≤ 5 2 ≤ x 4 ≤ 5 10 ≤ x 5 ≤ 50 S.t.

The optimal values obtained for the COA parameters are placed in the response surfaces (Eqs. 11, 12, 13, 14, 15, and 16) , and the values are calculated according to Table 5 (output).

In this step, using the input variables obtained in step 4, the COA for both the Ackley and Rastrigin functions, with the help of the MATLAB software, is executed 10 times, and the average of the obtained response variables (numerical) to compare with the result obtained in step 4 (analytical). For ease of comparison, the absolute value of the difference in the results is calculated in Table 5 . According to the results, it can be seen that Najafi and Behnoud's (2015) methods have a big difference in terms of numerical and analytical, while in the proposed method (Table 5) , this difference was small.

As mentioned in the literature review, one of the most common tools used to parameter setting of meta-heuristic algorithms is the Taguchi method. In the following, the Taguchi and RSM methods are implemented for both of the Ackley and Rastrigin functions, and for each function, responses are considered separately, and the results are shown in Fig. 4 . As mentioned earlier, this method only optimizes one response. Figure 4 (a) shows that considering the response value of the objective function, the best values for the input parameters are 20, 4, 3, 2, and 50. But if the execution time is examined (b), these optimal values change to 5, 2, 3, 2, and 10. Also, considering the results of response, the number of iterations (c) in Taguchi method are 20, 2, 3, 2, and 10. According to the different results of Taguchi method, the DSM method is implemented and the efficiency response is considered (d).

To determine the most important factor, RSM method was implemented, the results of which can be seen in Fig. 4 . When the value of the objective function is considered as the y 1 = 38.5052 + 0.0296916x 1 + 8.45343x 2 + 3.47633x 3 + 3.01223x 4 − 0.0229248x 5 − 0.0188956x 1 x 2 + 0.00478616x 1 x 3 + 0.00731104x 1 x 4

response, the most important parameter is mini number of egg (e), and considering the execution time and the number of iterations as the response (f and g), the most important input variables will be maximum number of cuckoo and number of clusters, respectively. Here, too, the efficiency response is used because of the different results, and as can be seen (h), all inputs are of almost equal importance. This process is performed similarly for the Rastrigin function, and the results are shown in Fig. 4 (i-p) . Figure 5 is based on the results obtained from the previous sections. As can be seen, the proposed DSM method performs better than other methods in all three responses of the value of the objective function, the execution time and the number of iterations of the algorithm. In the Ackley function, it was much better in terms of execution time, and in Rastrigin function, it was much better in terms of number of iteration. It should be noted that according to the different results of Taguchi method, the average of these values has been considered.

In order to compare the entire configuration created by the proposed DSM method with the Najafi and Behnoud and Taguchi methods in terms of both input and output values, the configuration efficiency obtained from the DSM method is compared with the configuration efficiency of other methods. The values for the DSM method, Najafi and Behnoud (2015) method, and Taghuchi on the Ackley function are, 45.3779, 3.64499, and 2.78653, respectively. Also for the Rastrigin function are 48.8649, 2.13645, and 2.03467. In order to better compare the three methods, we allow the DEA model to have an efficiency value of more than one. It can be seen that the configuration created by the DSM method in both Ackley and Rastrigin functions is more efficient than other methods, and this indicates the superiority of the DSM method. Due to the importance of parameter setting, the proposed method can be used in various problems such as repairs' scheduling (Al-Refaie & Al-Hawadi, 2021), supply chain network design (Hasani et al., 2021; Hosnavi et al., 2019; Khalili et al., 2017) , and vehicle routing Latorre-Biel et al., 2021) .

In this section, the reverse logistics model of COVID-19 waste presented by Shadkam (2021a Shadkam ( , 2021b Shadkam ( , 2021c is reviewed. After describing the model, it will be optimized using the proposed method of DSM and cuckoo optimization algorithm. 

The proposed model is a direct and reverse single-cycle and single-purpose logistics network that includes three parts: production, consumption, and recycling. These three sections are shown in Fig. 6 . The production part includes suppliers, factories, warehouses, and distribution centers. These departments are responsible for the production and production of coronavirus (vaccine) drugs. The raw materials are sent to the factory, and from the factory, coronavirus (vaccine) drugs are transported to distribution centers, while the rest are stored. Coronavirus drugs (vaccines) are delivered from distribution centers to the consumption department. Consumption includes hospitals, clinics, laboratories, and residential areas that use corona (vaccine) drugs. The recycling part includes a collection/ disinfection center, a recycling center, and a landfill. In the reverse flow, the wastes of coronavirus drugs (vaccine) are transferred from the consumption section to the collection/disinfection center. In this center, medical wastes are collected and disinfected, and from there, depending on the type of waste, it is transferred to distribution centers, factories, landfills, or recycling sites. It is assumed that corona vaccine wastes after disinfection are divided into two categories: recyclable and non-recyclable. Recyclable wastes are also divided into three categories.

A batch must be completely recycled and reused as a raw material. This group is recyclable and not dangerous and is sent to the recycling center for reuse in the factory. The second group is sent directly to the factory for further use in the production of corona medicine, and the third group is sent to the distribution center for use in the process of transporting drugs and vaccines. Non-recyclable Shadkam (2021a Shadkam ( , 2021b Shadkam ( , 2021c wastes are sent directly to the landfill for safe disposal and disposal. In this way, the proposed network will be a closed-loop logistics network.

Index of potential supplier centers i ∈ I J Index of fixed factory centers j ∈ J K Index of potential distribution centers k ∈ K L Index of fixed consumer centers (including hospital, clinic, laboratory, residential area)l ∈ L M Index of potential collection/disinfection centers m ∈ M P Index of potential recycle centers p ∈ P N Index of potential landfill centers n ∈ N If the collection/disinfection center m is constructed 1 otherwise 0, if the distribution center k is constructed 1 otherwise 0, if the recycle center p is constructed 1 otherwise 0, if the landfill center n is constructed 1 otherwise 0 X ij , X jk , X kl , X lm The amount of medical products flow from the supplier's center i to the factory's center j, the factory's center j to the distribution center k, the distribution center k to the consumption center l, the consumption center l to the collection/ disinfection center m Q jj

The amount of medical products flow from the factory's center j to its warehouse Q jk

The amount of medical products flow from the warehouse j to the distribution center k X mo , o {p, n, j, k}

The amount of return flow of medical products from the collection/disinfection center m to the recycle center p, landfill centers n , factory's center j, distribution center k

If the collection/disinfection center m is constructed 1 otherwise 0, if the distribution center k is constructed 1 otherwise 0, if the recycle center p is constructed 1 otherwise 0, if the landfill center n is constructed 1 otherwise 0 U j The amount of medical products in warehouse j as inventory

The mixed integer model for integrated reverse logistics network with the aim of reducing costs is presented as follows: In Eq. 19, the objective functions f 1 to f 5 represent fixed cost, direct transportation costs, warehouse transportation costs, reverse transportation costs, and inventory costs. These objective functions indicate the minimization of the cost of building facilities, the cost of transporting medical products directly and inversely, and inventory in accordance with Eqs. 38, 39, 40, 41, and 42. Constraints (20) and (21) ensure that all medical needs are met in direct flow and that all returned waste is collected from consumers in reverse flow. Constraints (22) to (25) are related to the flow level in the centers. Constraint (26) ensures that the output flow of the warehouse is less than the total revenue stream of the warehouse. Constraints (27) to (33) ensure that the flow is only between the centers where the facility is constructed and that the total flow in each center does not exceed its capacity. Constraint (34) ensures that at least one potential center is active. Constraints (35) ensure that the sum of the return waste coefficients is 1. Constraints (36) and (37) determine the types of model decision variables.

In this section, the parameters of the cuckoo optimization algorithm for use in the mentioned logistics network are examined. These parameters are examined and compared using two approaches, including the experimental approach and the DSM method. The optimal parameters of the cuckoo algorithm using the experimental approach are shown in Table 6 . Also, after implementing the DSM method in a similar way in the previous sections, the results are as shown in Table 6 . In order to compare the performance of both methods, the criteria of the value of the objective function, execution time, and number of iterations are calculated using the configuration of both methods (Table 6) . Also, the convergence diagram of the value of the objective function and the number of iterations for both approaches is shown in Fig. 7 . As can be seen, using the parameters created from the proposed method of the paper, a better objective function can be achieved in the same iterations with less execution time than the experimental approach, which shows the optimal performance of the proposed method. Also, in order to compare the efficiency of both methods, the values of Table 6 are considered as inputs, and the values of Table 6 are considered as outputs, and after implementing the CCR model of data envelopment analysis, the efficiency values are calculated. To better differentiate between the two methods, the CCR model is allowed to obtain performance values greater than one. The results show an efficiency value of 1.23456 for the proposed method and efficiency value of 0.98753 for the experimental method. Therefore, it can be seen that the efficiency of the proposed method is better than the experimental method.

In this paper, the hybrid DSM method was presented by combining two methods of RSM and DEA. One of the advantages of the DSM method is the reduction in the number of response surfaces to one efficiency surface, and instead of optimizing multiple response surfaces for the problem outputs, only efficiency surface optimization is considered. With this proposed method, in addition to optimizing the objective function of the problem, it is possible to maximize its efficiency simultaneously. In this paper, the hybrid DSM method was used to set the optimal parameters of the COA. In order to evaluate the performance of the DSM method, two types of problems were investigated, which included standard functions and COVID-19 waste management problem. Standard functions include Ackley and Rastrigin functions. In order to validate the proposed method, the Taguchi and the Najafi and Behnoud methods were implemented and the results of the proposed method were compared with these methods. The numerical results indicate the higher efficiency of the proposed DSM hybrid algorithm. Then, the cuckoo algorithm parameter was adjusted using the proposed method for the inverse logistics . 7 The implementation of the cuckoo algorithm on the proposed inverse logistics network based on the a experimental method and b DSM method problem of COVID-19 waste, and the results were compared with the experimental approach which shows the superiority of the proposed method over the experimental method according to the objective function, execution time and efficiency value. Finally, it can be said that by setting the parameters of meta-heuristic algorithms scientifically, errors related to trial and error can be avoided. As future research, the DSM method can be used to parameterize other meta-heuristic algorithms, and by assigning appropriate values to the effective parameters in the algorithms, a suitable and acceptable solution can be achieved in the shortest time and with the least error and high reliability. Also, the process in this paper can also be applied similarly to other real-world problems.

Author contribution The whole process of the preparing the paper, including coding, modeling, and writing, etc., has been done by Elham Shadkam.

Code availability All paper codes are written by the author with the Lingo and Mititab software and can be presented if needed.

The author agrees to participate and publish.

The author declares no competing interests.

Medical waste-sorting and management practices in five hospitals in Ghana

The study of cuckoo optimization algorithm for production planning problem

Optimal fuzzy repairs' scheduling and sequencing of failure types over multiple periods

Simultaneous optimization approach for combined control-structural design versus the conventional sequential optimization method

A cost optimization model for hazardous medical waste management in Portugal

Stochastic modeling and optimization of medical waste collection in Northern Jordan

An integrated multi-criteria Taguchi computer simulation-DEA approach for optimum maintenance policy and planning by incorporating learning effects

Tuning search algorithms for real-world applications: a regression tree based approach

Tuning PSO parameters through sensitivity analysis

The new hybrid COAW method for solving multi-objective problems

Reverse logistics optimisation for waste collection and disposal in health institutions: the case of Turkey

A study on flowshop scheduling problem combining Taguchi experimental design and genetic algorithm

Auto-tuning SOS algorithm for two-dimensional orthogonal cutting optimization

Tuning multiobjective evolutionary algorithms on different sized problem sets

Optimal controller parameter tuning of PSS using sine-cosine algorithm

Tuning parameters in deep belief networks for time series prediction through harmony search

A parameter tuning methodology for metaheuristics based on design of experiments

A framework for assessing the circularity and technological maturity of plastic waste management strategies in hospitals

Clustering approach using artificial bee colony algorithm for healthcare waste disposal facility location problem

Tuning metaheuristic algorithms using mixture design: application of sunflower optimization for structural damage identification

A hybrid COA-DEA method for solving multi-objective problems

A multi-objective optimization approach for green and resilient supply chain network design: a real-life case study

Wiley Online Library Homayouni Z, Pishvaee MS (2020) A bi-objective robust optimization model for hazardous hospital waste collection and disposal network design problem

A resilient supply portfolio considering political and disruption risks

Parameter tuning for meta-heuristics

2020) Interplay between reverse logistics and circular economy: critical success factors-based taxonomy and framework

A modified cuckoo optimization algorithm for engineering optimization

A reverse supply chain for medical waste: a case study in Babol healthcare sector

Reverse logistics network design for medical waste management in the epidemic outbreak of the novel coronavirus (COVID-19)

Why tuning the control parameters of metaheuristic algorithms is so important for fair comparison?

A bi-objective model for relief supply location in postdisaster management

Integrated production-distribution planning in two-echelon systems: a resilience view

Combining simheuristics with Petri nets for solving the stochastic vehicle routing problem with correlated demands

Parameter setting in evolutionary algorithms

Flower pollination algorithm parameters tuning

Scheduling hybrid flowshops with sequence dependent setup times to minimize makespan and maximum tardiness

Consolidated technique of response surface methodology and data envelopment analysis for setting the parameters of meta-heuristic algorithms-case study: production scheduling problem

A stochastic inventory routing problem for infectious medical waste collection

The mathematical model, implementation and the parameter-tuning of the African buffalo optimization algorithm

A discrete and improved bat algorithm for solving a medical goods distribution problem with pharmacological waste collection

Experimental thermodynamic analysis of air-based PVT system using fins in different materials: optimization of control parameters by Taguchi method and ANOVA

Significant density change in medical waste during Covid19 pandemic and effects on medical waste management in Istanbul. 5th Eurasia Waste Management Symposium

Design of experiments for the tuning of optimisation algorithms

Analyzing heuristic performance with response surface models: prediction, optimization and robustness

Globale Optimierung mit parallelen Evolutionsstrategien Diplomarbeit

Parameter tuning of the HCSCROCFO-3Opt algorithm for solving the capacitated vehicle routing problem

Face mask and medical waste disposal during the novel COVID-19 pandemic in Asia

Tuning the parameters of a memetic algorithm to solve vehicle routing problem with backhauls using design of experiments

Cuckoo optimization algorithm in reverse logistics: a network design for COVID-19 waste management

Evaluation the efficiency of cuckoo optimization algorithm

The optimization of bank branches efficiency by means of response surface method and data envelopment analysis: a case of Iran

Multi-objective simulation optimization for selection and determination of order quantity in supplier selection problem under uncertainty and quality criteria

A novel improved cuckoo optimisation algorithm for engineering optimisation

Portfolio selection by the means of cuckoo optimization algorithm

A hybrid COAε-constraint method for solving multi-objective problems

The new TOPCO hybrid algorithm to solve multi-objective optimisation problems: the integrated stochastic problem of production-distribution planning in the supply chain

The problem of resource leveling in multi-project mode by cuckoo optimization algorithm

A novel two-phase algorithm for a centralized production planning problem by symmetric weighted DEA approach: a case study in energy efficiency

Network model and optimization of medical waste reverse logistics by improved genetic algorithm

Multicriteria optimization for infectious medical waste collection system planning. Pract Period Hazard Toxic Radioact

Data mining-assisted parameter tuning of a search algorithm

Experimental investigation and multi response optimization of turning process parameters for Inconel 718 using TOPSIS approach

Fine-tuning of neural computing using whale optimization algorithm for predicting compressive strength of concrete

Parameter tuning with chess rating system (CRS-Tuning) for meta-heuristic algorithms

Animal communication: when I'm calling you, will you answer too?

Identifying good algorithm parameters in evolutionary multi-and many-objective optimisation: a visualisation approach

Disinfection technology of hospital wastes and wastewater: suggestions for disinfection strategy during coronavirus Disease 2019 (COVID-19) pandemic in China

A parameter set design procedure for the simulated annealing algorithm under the computational time constraint

Fine-tuning a tabu search algorithm with statistical tests

Optimizing the sum of maximum earliness and tardiness of the job shop scheduling problem

A singlemachine scheduling problem with multiple unavailability constraints: a mathematical model and an enhanced variable neighborhood search approach

A parallel machine scheduling problem with two-agent and tool change activities: an efficient hybrid metaheuristic algorithm