Interactive Reference Point Procedure Based on the Conic Scalarizing Function

In multiobjective optimization methods, multiple conflicting objectives are typically converted into a single objective optimization problem with the help of scalarizing functions. The conic scalarizing function is a general characterization of Benson proper efficient solutions of non-convex multiobjective problems in terms of saddle points of scalar Lagrangian functions. This approach preserves convexity. The conic scalarizing function, as a part of a posteriori or a priori methods, has successfully been applied to several real-life problems. In this paper, we propose a conic scalarizing function based interactive reference point procedure where the decision maker actively takes part in the solution process and directs the search according to her or his preferences. An algorithmic framework for the interactive solution of multiple objective optimization problems is presented and is utilized for solving some illustrative examples.


Introduction
A multiobjective optimization (MOP) problem can be stated as min ∈ ( 1 ( ) , . . . , ( )) , where is the feasible region in decision space; is the feasible region in the objective space, where ∈ if and only if there exists an ∈ such that = ( 1 ( ), . . . , ( )).
The vector of values, , is usually referred to as the nadir point while the vector of values, , is often labeled the ideal point.
A reference point consists of desirable values for each objective function. For decision makers (DMs), reference points are a natural way of expressing desires in solutions because DMs do not have to learn to use new, artificial concepts. Instead, objective function values are used that as such are meaningful and understandable for DMs.
The use of reference points appeared in the early development of multiple objective programming as part of the work of Charnes and Cooper [1] on goal programming. 2 The Scientific World Journal Wierzbicki [2] produced seminal research on reference point methods, including an investigation of the characteristics of various achievement functions for allowing the search for attractive efficient solutions to be controlled by reference points. These achievement functions were designed to have a significant advantage over goal programming by producing only efficient, or Pareto-optimal, solutions. Reference point methodology provides the foundation for many interactive search procedures in multiple objective programming. Over the last decades, various interactive methods and decision support systems have been developed to deal with MOPs. Among interactive approaches, methods using reference points (for the idea see, e.g., [3,4]) have been popular (for some comparative studies, see, e.g., [5,6]) because of their straightforward nature. Examples of methods utilizing reference points include the reference point method [4], the visual interactive approach [7], STOM [8], GUESS [9], and light beam search [10]. In addition, methods based on classification are closely related to reference point methods because a reference point can be formed once a classification has been made [11,12]. Interactive methods enable the reduction of the computational effort and aid the DM in the decision process. In interactive methods, the set of efficient solutions is explored by a progressive articulation of the DM's preferences. This is a shared feature of all interactive methods, but there are different paradigms.
Interactive procedures are closely related to scalarization functions. Scalarization means combining different objectives into a single one such that the obtained single objective optimization problem allows finding (all) efficient (or properly efficient) solutions of the initial multiobjective problem. A variety of scalarization methods for finding efficient solutions of multiple objective programs MOPs have been developed. Some of the methods were designed specifically for linear problems and others work well only on problems with concave objective functions and a convex feasible region. However, most of the mathematical programming models of real life problems have discrete variables. Since the set of efficient solutions for problems with discrete variables is not convex, weighted sums of the objective functions do not provide a way of reaching every efficient solution. Besides the supported, there exist unsupported efficient pointspoints that are dominated by convex combinations of other efficient points. Conic scalarization based programs have the advantage over weighted-sums programs of being able to reach not only supported, but also unsupported efficient points [13][14][15]. The conic scalarization function falls into the class of achievement scalarizing functions [16,17]. Depending on the properties of these functions, they can be used to identify (weakly, properly) efficient solutions. Most of them differ from the conic scalarizing functions by using a norm term of higher order (e.g., the Euclidean or maximum norm) in addition to the weighted-sum term. The conic scalarizing functions presented in [13,14] guarantee a complete characterization of all properly efficient solutions.
A general characterization for the set of Benson proper efficient point was first proposed by Gasimov [13,14]. Gasimov introduced an explicit class of increasing convex functions which serve for combining different objectives into a single one without any restrictions on objectives and constraints of the problem under consideration, such as convexity and/or boundedness. These functions are used to characterize the Benson proper efficient points of nonconvex multiobjective problems in terms of saddle points of scalar Lagrangian functions introduced in the paper [13,14]. The conic scalarization method has been successfully applied to several multiobjective problems. Firstly, the conic scalarization method was applied by Ozdemir and Gasimov [18] to the nonconvex multiobjective faculty course-assignment problem with nine objective functions. This integrated approach consists of an analytic hierarchy process (AHP), the conic scalarization, and the modified subgradient algorithm. In the study, the weights of the objectives were calculated by using AHP. The calculated weights were used for combining the objectives by using the conic scalarization method. Ismayilova et al. [19] proposed a similar approach to solve the multiobjective faculty course-time-slot assignment problem. AHP and analytic network process (ANP) were used to determine the weights of conflicting objectives. Efficient points corresponding to the obtained weights by using AHP and ANP have been calculated and the results compared. The conic scalarization method has also been applied to the 1.5dimensional assortment problem with two objective functions by Gasimov et al. [20]. In the study, the authors focused on the possible nonsupported or "hidden" efficient points. Since the problem under consideration involves integer variables, it becomes nonconvex. In these situations a DM may not be sure whether all the efficient points corresponding to own preferences (weights) are found by using the weightedsum scalarization. To overcome this drawback, Gasimov et al. [20] proposed a way to obtain nonsupported efficient points based on the conic scalarization function.
Another application area of the conic scalarization method is the multiobjective portfolio optimization [17,21]. Ehrgott et al. [17] compared multiattribute utility theory (MAUT) with the conic scalarization method. The calculated weights by using UTADIS method were used for combining the objectives in the conic scalarization method. The nonsupported or "hidden" efficient points can be obtained by using the conic scalarization method. Ustun and Kasimbeyli [21] calculated different efficient points of a multiobjective portfolio optimization problem based on combined forecasts with 11 objectives.
These approaches based on the conic scalarizing function take place either in the class of a priori methods or the class of a posteriori methods. The conic scalarization method is used as a posteriori method for the multiobjective portfolio optimization problem and the 1.5-dimensional assortment problem. A posteriori methods can also be called methods for generating efficient points. After the set of efficient points (or a part of it) has been generated, it is presented to the DM, who selects the most preferred among the alternatives. The drawbacks here are that the generation is usually computationally expensive and sometimes at least partly difficult or impossible. On the other hand, it is hard for the DM to select from among a large set of alternatives. If there are only two objective functions in a mathematical model, a set of efficient points can be generated parametrically (see, e.g., [20]). When The Scientific World Journal 3 there are more than two objectives, the problem becomes more complicated [17]. The conic scalarization method is also used as an a priori method for the nonconvex multiobjective faculty course-assignment problem and the multiobjective faculty course-time-slot assignment problem. In the a priori methods, the DM must specify the preferences, hopes, and opinions before the solution process. The main difficulty is that the DM does not necessarily know beforehand what is possible to attain in the problem and how realistic such expectations are. The existing conic scalarization approaches in the literature do not allow changes of the weights and/or reference points by the DM during the solution process according to learning. Additionally, these approaches do not support the DM in the search for more satisfactory solutions. If the DM wants to trade off among some objective function values for a current solution, she/he should tune the weights of the objective functions and/or the reference point. Unfortunately, the tuning process for obtaining more satisfactory solutions is usually computationally expensive or sometimes impossible using the existing conic scalarization approaches in the literature.
The main aim of this paper is to proposean interactive reference point procedure based on a conic scalarizing function to cope with these drawbacks to the use of the conic scalarizing function in the literature. We use the main idea of interactive classification based multiobjective optimization methods [3,22]. The main idea of interactive classification based multiobjective optimization methods is that the DM examines the values of the objective functions calculated at a current efficient point and classifies the objective functions into different classes. This means that the DM is asked to indicate (by means of a classification) what kind of a solution would be more satisfactory than the current one [3,11,22]. Many scalarizing functions produce only weakly efficient solutions but such solutions are not necessarily interesting to a DM because some objective function may be improved without impairing any of the others [6]. The conic scalarization method guarantees obtaining a Benson proper efficient solution according to the multiobjective problem without solving an additional subproblem. The proposed procedure allows the DM being asked to say which objective functions have acceptable values, which should be improved, and which could be impaired. In this way, the DM has a flexible decision-making procedure to move around in the set of efficient points towards more satisfactory solutions. By using classification, we can avoid the wellknown drawbacks related to expressing preference information in the form of weighting coefficients or generating the set of efficient points in the existing conic scalarization approaches.
The rest of this paper is organized as follows. We present some concepts, basic theorems, related results, and methods in Section 2. In Section 3, we introduce classification based conic scalarizing functions and the related single objective subproblems with some theoretical results for them. The conic scalarization based interactive reference point procedure is given in Section 4. Then, we apply the proposed procedure to some illustrative or comparative examples in Section 5. The paper is concluded in Section 6.
Definition 2. Let be a nonempty subset of .   ∈ is called an efficient solution of (1) if ( ) ∈ min( ( )); ∈ is called a proper efficient point of (1) (in the sense of Benson), if ( ) ∈ − min( ( )). Let Theorem 4 (see [13]). Suppose that for some ( , ) ∈ an element ∈ is an optimal solution to the following scalar minimization problem: Then, the solution ∈ is a proper efficient solution of (1).
Theorem 5 (see [13]). Suppose that ∈ is a proper efficient solution of (1). Then there exists an element ( , ) ∈ , such that ∈ is an optimal solution to the following scalar minimization problem: Theorem 4 asserts that any solution of the (scalar) problem (4) is an efficient solution of the problem (1). On the other hand, Theorem 5 claims that every efficient solution of (1) can be calculated by solving a scalar problem of the form (5) for some pair ( , ) ∈ . Thus, these theorems assert that the problem (1) can be scalarized in the form (4) and/or (5), and all efficient solutions of (1) can be calculated by solving scalar problems of these forms. Regarding the solutions of problem (1), note that in the case where the signs of the objective functions remain unchanged in the set of feasible solutions, the absolute values in (4) are not essential in the whole expression, and it reduces to an expression 4 The Scientific World Journal representing the weighted-sum scalarization. This situation can be overcome as follows. It is well known and easy to show that the following two problems have the same set of efficient solutions: where 1 , . . . , are arbitrary fixed numbers. By choosing these numbers in the interior of ranges of 1 , . . . , , respectively, we can replace the objective functions by new ones for which the conic scalarization approach will work more effectively. When the multiobjective problem under consideration is nonconvex, by taking different values for the parameter and the numbers 1 , . . . , , different efficient solutions corresponding to the same set of weights 1 , . . . , can be obtained by solving the scalarized version of the problem (7) as follows: Additionally, the analyst must call attention to properties of the conic scalarizing function (8).
Note that the conic scalarizing function (5) for = 0 or ≤ reduces to the weighted-sum scalarization. Therefore, nonsupported efficient solutions can only be found with > 0 and > . We This is equivalent to the weighted-sum scalarizing function.

A Priori Method Based on the Conic Scalarizing Function.
A priori methods based on the conic scalarizing function were proposed by Ozdemir and Gasimov [18] and Ismayilova et al. [19] to solve the multiobjective faculty courseassignment problems. The main steps of these methods can be given as follows.
Step 1. Construct a multiobjective model for a problem under consideration.
Step 2. Calculate the weights of the objectives using the AHP and/or the ANP proposed by Saaty [23,24]. Choose a reference point.
Step 3. Synthesize the multiobjective problem by using the weights and the reference point determined at Step 2 for the conic scalarizing function to obtain a single objective problem.
Step 4. Solve the synthesized problem for the different values of . Show the obtained efficient points to DM to obtain a final solution.

A Posteriori Method Based on the Conic Scalarizing Function.
A posteriori methods were used to solve multiobjective portfolio optimization problems by Ehrgott et al. [17] and Ustun and Kasimbeyli [21] and the 1.5-dimensional assortment problem with two objective functions by Gasimov et al. [20]. The main steps of these approaches can be given as follows.
Step 1. Generate a group of well-dispersed weight vectors and choose a reference point from the interior of ranges of the objective functions. Scalarize the given problem by using the conic scalarizing function.
Step 2. Solve the scalarized problem for each weight vector. Show the obtained efficient points to DM for a final solution.
In the existing a posteriori methods, authors tried to reach the different nonsupported or "hidden" efficient points that cannot be calculated by the weighted-sum scalarization. To reach possible nonsupported efficient points, different weight vectors and reference points were firstly generated. The scalarized problem is solved by using these parameters. Then obtained hidden efficient points are diagrammatically illustrated in a figure. The most systematical approach among these studies is given by Gasimov et al. [20]. Gasimov et al. [20] propose the following way to obtain a group of welldispersed weight vectors and to select a reference point in Step 1 for a biobjective programming problem.
Authors utilized formula (8) for all possible integer triples ( , 1 , 2 ) ∈ with 1 , 2 = 1, 2, . . . , 50, 1 + 2 = 50, and = 1, 2, . . ., min{ 1 , 2 }. The reference points initially can be arbitrarily chosen between the "relatively distant" from each other. After calculating additional efficient points, the intervals for choosing the nonzero values for reference points in the conic scalarizing function can be sequentially tightened and DM can find other nonsupported efficient solutions (if they exist) in this way, by which the authors obtained different linearly nonsupported efficient points.

Classification Based Conic Scalarizing Function and Related Subproblems
In this section, we introduce some background for the classification based conic scalarizing function to be used in this study. We use the general idea of classification in multiobjective optimization given by Miettinen [3,22]. The idea of interactive classification based multiobjective optimization methods is that the DM examines the values of the objective functions calculated at a current efficient solution and classifies the objective functions into different classes. This means The Scientific World Journal 5 that the DM is asked to indicate (by means of classification) what kind of a solution would be more satisfactory than the current one. In this paper, we use up to four different classes and they are for functions whose values (2) should be decreased down to some aspiration level The DM is asked to specify the aspiration levels for ∈ ≤ satisfying ≤ ( ) and the upper bounds for ∈ ≥ such that ≥ ( ). The difference between the classes < and ≤ is that functions in < are to be minimized as far as possible but functions in ≤ only to the aspiration level. The idea is that with classification, the DM directs the solution process in order to find the most preferred efficient solution. That is why a classification is feasible only if at least some of the objective functions should improve and some are allowed to increase from their current levels, that is, four different classes. We can also assume that in a feasible classification, an objective function that has reached its ideal value cannot be assigned to either of the classes < ∪ ≤ .
We can update the conic scalarization function by using the classification information expressed by DM. We update the conic scalarizing function by using the original objective functions and the preference information about the four different classes specified. We propose the following single objective subproblem: ∈ .
The constraints (12) provide lower bounds for ( ) for all ∈ ≤ , because these objective functions should be decreased down to . There is another way to provide that the functions should be decreased down to for all ∈ ≤ . Weights of ( ) for all ∈ ≤ can be taken equal to the value of parameter in (8). In this way, the constraints (12) are eliminated. This result can be shown as follows: If is taken equal to in (14), we obtain Thus it provides that the functions should be decreased down to for all ∈ ≤ . But the taking of = for any contradicts the theoretical results of the conic scalarization such that ( , ) ∈ . To eliminate this contradiction, we take the related weighting values to be = + > 0, where is a small scalar for all ∈ ≥ ∪ ≤ ∪ = . On the other hand, the values ( ) for ∈ ≥ should be allowed to increase up to . The DM can specify the aspiration levels = for ∈ ≥ for this requirement. Additionally, the values ( ) for ∈ = for a current efficient solution are required to remain at its current value, ( ).
We have also used a denominator to scale each term in the objective function of (8) to be of a similar magnitude. This aims to capture the preferences of the DM better. Note that by the definitions of ideal and nadir point, we have − > 0 for all = 1, . . . , .
The subproblem (16)- (19) combines the conic scalarization method with the -constraint method. It is a hybrid method which considers both of DM's aspiration levels and reservation levels by the classification point of view.
is an optimal solution of the subproblem (16)- (19) if and only if * ∈ is a proper efficient solution of (1).

Conic Scalarizing Function Based Interactive Reference Point Procedure
After summarizing the necessary mathematical background, we can introduce the conic scalarization based interactive reference point procedure.
Before the solution process starts, some information is given to the DM about the problem. If possible, the ideal point and the nadir point are presented to illustrate the ranges of the set of efficient points . Another possibility is to minimize and maximize the objective functions individually in the feasible region. Thus if the feasible region is bounded, then the DM is informed about the best and worst values of the objective functions. If calculating the ideal or nadir points is difficult in practice and/or any element of these points is infinite, the denominators − in the objective function (16) can be removed for all . The theoretical result given in Theorem 7 is valid in this situation.
The main steps of the conic scalarization based interactive reference point procedure are the following.
Step 1 (initialization). Set = 0. Convert the maximization objectives into minimization ones by multiplying the objective functions by −1. Present information about the problem to the DM. If possible, calculate the ideal and nadir points. Present them to the DM. Ask the DM to specify a finite reference point ∈ .
Step 2 (weights determination). Ask the DM to specify weights ∈ of the objective functions for ∈ < . If = 0, then all objective functions are in the class < and set = . The analytic hierarchy process (AHP) is proposed to the DM for determining the weights of the objective functions ∈ [23] in this step. Ask the DM to construct a pairwise comparison matrix of the relative contributions of each objective function by using Saaty's 9-point scale as follows: where is the relative importance of objective function compared to objective function .
Compute the consistency ratio (CR) of the pairwise comparison matrix. Accept the matrix if CR is less than 0.10. Otherwise, check the judgments on the relative importance of the objective functions to eliminate the inconsistency until CR is less than 0.10.
Compute the priorities for each objective function. Use these priorities as the weights of objective functions = ( 1 , . . . , ) after normalization.
Step 3 (most preferred solution). If = 0, then all objective functions are in the class < and use ( ) = , where is nadir point (or the worst values of the objective functions in the feasible region) and ∈ , by satisfying < ( ).
Present the solutions of the subproblems according to the values of to the DM. Ask the DM to determine the most preferred solution among the calculated solutions for different values of . Set = , where the th efficient point is the most preferred by the DM. Let be an efficient solution related to .
Step 5 (decision). Present the alternatives to the DM. If the DM finds one of the + 1 efficient points satisfactory, then the DM chooses it as the final solution. Otherwise, set = + 1. Ask the DM to classify the objective functions into four different classes such as = , < , ≤ , and ≥ . Let be the number of objectives in the class < . Also ask the DM to specify aspiration levels for ∈ < ∪ ≤ satisfying < ( −1 ) and upper bounds such that ≥ ( −1 ) for ∈ ≥ . Go to Step 2. It is clear that if increases, then the DM can select a more sensible efficient solution in Step 3. But large values of The Scientific World Journal 7 require more computational time. In this study the parameters of the conic scalarization function are interactively determined by the DM. The weights of the objectives are calculated by using AHP at each iteration of the procedure. The value of the parameter is chosen by the DM after solving subproblems. The proposed procedure must solve + subproblems at each iteration. Thus at each iteration, the DM evaluates + efficient solutions (some of these can be the same according to the shape of the efficient frontier) and learns whether the goals are realistic or not. The advantage of perturbing the reference point in Step 4 is that the DM gets a better conception about possible solutions. If the reference point is far from the set of efficient points, the DM gets a wider description of the set of efficient points, and if the reference point is near the set of efficient points, then a finer description of the set of efficient points is given.

Illustrative Examples
We now examine the conic scalarization based interactive reference point procedure in the context of some nonconvex MOP problems. Utility functions were used to simulate DM preferences with different weights of objective functions. Thus the solution quality of a reference point and an efficient point can be evaluated and compared by using DM's utility function. The mathematical models in this study are solved by using the Lingo 11.0 solver. There is no limitation for the number of variables, integer variables, and constraints. Lingo 11.0 can solve linear, integer, mixed integer, and nonlinear programming models [25]. We also compute the priority for each objective function and consistency ratios of pairwise comparison matrices by using the trial version of Expert Choice 11.5. The trial version of Expert Choice 11.5 can be downloaded at the URL (http://expertchoice.com/). An Intel Core 2 Duo 2.27 GHz Notebook PC with WINDOWS 7 operating system on one processor is used in the computations.
The objective functions 1 ( , ), 2 ( ), and 3 ( ) are, respectively, the sum of the material cost and the order cost, the expected number of defective products, and thetotal value of purchasing. We assume that the DM wants to join the interactive decision-making process and that DM's implicit utility function is UT( ( , )) = (− 1 ( , ) − 2 ( ) − 3 ( ))/3. Let be the feasible set determined by the constraints (30)-(36). The main steps of the conic scalarization based interactive reference point procedure are simulated as follows.
Step 2 (weights determination). Because = 0, all objective functions are in the class < and set = 3. The DM constructs a pairwise comparison matrix of the relative contribution of each objective function by using Saaty's 9-point scale as follows: Although CR = 0 is less than 0.10, the pairwise comparison matrix is consistent. The calculated priorities of the objective functions are 0.33, 0.33, and 0.33, respectively. Use these priorities as the weights of the objective functions 0 = ( 0 1 , 0 2 , 0 3 ) = (0.33, 0.33, 0.33).
The efficient points according to the perturbed reference points are given in Table 3.
Step 5 (decision). The efficient points obtained in Table 3 are the same with 0 . The DM wants to continue the process for a more attractive solution. Set = + 1. The DM classifies the objective functions into two different classes. There are two objective functions 1 ( , ) and 3 ( ) in the class < . The objective function 2 ( ) is in the class ≥ . The DM specifies aspiration levels 1 1 = 118000 and 1 3 = −2600 for ∈ < satisfying 1 < ( 0 ) and the upper bounds 2 = 38 ≥ ( 0 ) = 33.5 for ∈ ≥ . Go to Step 2.

Second Iteration
Step 2 (weights determination). The DM constructs a pairwise comparison matrix of the relative contribution on each objective function by using Saaty's 9-point scale as follows: Although CR = 0 is less than 0.10, the pairwise comparison matrix is consistent. The calculated priorities of objective functions are 0.5 and 0.5, respectively. Use these priorities as the weights of objective functions 1 = ( 1 1 , 1 3 ) = (0.5, 0.5).
Ask the DM to determine the most preferred solution among the calculated solutions for different values of . Set Step 4 (alternative solutions). The analyst calculates three other efficient points by solving the subproblem in Step 3 for = 9 = 0.45 with perturbed reference points 1 ( ) = 1 + 1 , where 1 = ‖ 1 − 1 ‖ = 414.8 and is the th unit vector for = 1, 2, 3.
The obtained efficient points according to the perturbed reference points are given in Table 5.
Example 2. Let us consider a mixed integer nonlinear multiobjective problem with 5 objectives, 40 binary variables, 40 positive variables corresponding to assets, and 3 constraints [17,27]. Let be the number of available assets, and the fraction of the available capital invested in asset and the expected return of asset for = 1, . . . , , and the covariance of the returns of assets and . For asset The Scientific World Journal 9  ∈ {1, . . . , }, let 12 be the 12-month performance, 36 the 36-month (long-term) performance (expected return), ≥ 0 the relative annual dividend, and ∈ {1, . . . , } the number of stars assigned to asset (one star indicates relatively poor performance of assets and five stars indicate very good performance). We have used a dataset of = 40 investment funds from Standard and Poor's 1999 database. Consider the following five objective functions.
is the 12-month performance. This objective function is a measure for the short-term expected return. 02 ( ) = ∑ =1 36 , the 3-year performance, is a measure for the long-term expected return.
represents the relative annual dividend of a portfolio.
is the average star ranking of portfolio risk. Standard and Poor's Fund Service evaluates the performance of most investment funds contained in their data based on an annual basis which results in a performance ranking.
is the usual variance measure of portfolio risk.
Let , = 1, . . . , , denote binary variables with = 1 if and only if asset is contained in the portfolio. It is assumed that exactly assets must be held in the portfolio. Thus, we consider the cardinality-constrained ( = 10) portfolio optimization problem of size = 40 in Example 2: subject to This problem is an extended version of the Markowitz meanvariance model. We assume that DM's implicit utility function is UT( ( )) = −0.163 1 ( ) − 0.003 2 ( ) − 0.032 3 ( ) − 0.535 4 ( ) − 0.268 5 ( ). Let be the feasible set determined by the constraints (46)-(50). The main steps of the conic scalarization based interactive reference point procedure are simulated for Example 2 as follows.

First Iteration
Step 1 (initialization). Set = 0. The analyst or computer program calculates the ideal point and the    Table 6.
The solutions of the subproblems according to the values of are given in Table 7.
Ask the DM to determine the most preferred solution among the calculated solutions for different values of .  Table 7 is the most preferred by the DM.
The efficient points obtained are given according to the perturbed reference points in Table 8.
Step 5 (decision). The efficient points obtained are the same as with 0 given at
We compare the proposed interactive method with the existing approaches based on the conic scalarization by using the example. The existing approaches and the proposed method are employed to maximize DM's implicit utility function.
Case 1. Firstly, we applied the a priori method given in Section 2.1 to Example 3. The objective functions are transformed to the minimization form by multiplying −1.
Step 1. We consider the multiobjective problem given in Example 3.
Step 2. The weights of the objectives are calculated by using the AHP. We assume that DM constructs a pairwise comparison matrix of the relative contribution on each objective function by using Saaty's 9-point scale as follows: Step 3. We synthesized the multiobjective problem by using the weights and the reference point determined at Step 2. The obtained single objective problem due to the value of can be given as follows: Step 4. The synthesized problems were solved for the different values of such as 0.1, 0.2, 0.3, and 0.4. The obtained efficient points are given in Table 9. It is clear that DM chooses the efficient solution (0.7995, 0.6007) with the highest utility value 0.6441 among the efficient solutions given in Table 9.
Case 2. We applied the a posteriori method given in Section 2.2 to Example 3.
Step 2. The scalarized problem is solved for each weight vector. The highest utility value among the 650 number of  Case 3 (the proposed interactive procedure). The objective functions are transformed to the minimization form by multiplying −1.
Step 2 (weights determination). Because = 0, all objective functions are in the class < and set = 2. The DM constructs a pairwise comparison matrix of the relative contribution of each objective function as follows: Although CR = 0 is less than 0.10, the pairwise comparison matrix is consistent. The calculated priorities of the objective functions are 0.571 and 0.429, respectively. Use these priorities as the weights of the objective functions 0 = ( 0 1 , 0 2 ) = (0.571, 0.429).
The solutions of the subproblems according to the values of are given in Table 10.
Ask the DM to determine the most preferred solution among the calculated solutions for the different values of . Set = 1 = 0, where the first efficient point 0 = (0.7995, 0.6007) is the most preferred by the DM.
The efficient points according to the perturbed reference points are given in Table 11.
Step 5 (decision). The efficient points obtained in Table 11 are the same with 0 . The DM wants to continue the process for a more attractive solution. Set = + 1. The DM classifies the objective functions into two different classes. There is the second objective function 2 ( ) in the class < . The objective function 1 ( ) is in the class ≥ . The DM specifies aspiration levels 1 2 = −0.62 for ∈ < satisfying 1 < ( 0 ) and the upper bounds 1 = −0.75 ≥ 2 ( 0 ) = −0.7995 for ∈ ≥ . Go to Step 2.

Second Iteration
Step 2 (weights determination). The DM constructs a pairwise comparison matrix of the relative contribution on each objective function by using Saaty's 9-point scale as follows:  Although CR = 0 is less than 0.10, the pairwise comparison matrix is consistent. The calculated priority of the objective function is 1. Use this priority as the weight of the objective function 1 = ( 1 2 ) = (1).
Ask the DM to determine the most preferred solution among the calculated solutions for different values of . Set = 0 = 0, where the first efficient point 1 = (0.75, 0.6614) is selected as the most preferred efficient point by the DM.
The obtained efficient points according to the perturbed reference points are given in Table 13.
Step 5 (decision). The DM finds one of the two efficient points as satisfactory, and then DM chooses 1 (2) = (0.75, 0.6614) as the final efficient point.
The results obtained by using the a priori method, the a posteriori method, and the proposed method are summarized in Table 14.
As it can be seen in Table 14, the highest utility value among the methods is obtained by using the proposed  method. The proposed method allowed DM to reach the more satisfactory solution. It is clear that the existing approaches are either computationally expensive or not flexible for DM's learning process. In the a posteriori method, 650 efficient points should be compared by DM for a simple biobjective problem. On the other hand, although the more satisfactory efficient solution exists in the feasible region for Example 3, the a priori method causes DM to have to select the less satisfactory efficient solution.

Conclusions
In this study, a novel interactive approach called the conicscalarizing function based interactive reference point procedure is proposed. The proposed procedure aims to eliminate the drawbacks in current practice of the use of the conic scalarization based methods. This procedure allows a flexible multiobjective decision-making process to select a more preferred efficient solution for the DM. The basic technique of the proposed method is the subproblem (16)- (19) which combines the conic scalarization method with the -constraint method. It is a hybrid method which considers both of DM's aspiration levels and reservation levels by the classification point of view. DMs direct the proposed procedure by using four different classes to obtain more satisfactory solutions. In this way, the proposed procedure allows tradeoffs among the objectives' values by classifying. The theoretical background of the proposed method is given in detail. The proposed method was applied to illustrative examples. It is also compared with the existing methods according to the multiplicative-additive utility function by using the comparative example. The highest utility value among the methods is obtained by using the proposed method. Further directions of the study can be given as follows: from the theoretical point of view, the proposed algorithm can be combined with different optimization techniques for effectively solving nonconvex subproblems. Also, different sampling techniques can be adapted to obtain different efficient solution sets. The main difficulty about the use of the conic scalarization technique is the tuning of the value of the parameter according to the preferences of the DM. In this study, we proposed a sampling method related to values. Different search methods can be used for this purpose. For example, a kind of interval reduction technique can be used. From the practical point of view, the proposed algorithm can be applied to different test problems and compared with other interactive procedures. The weights of the objectives can be calculated by using the analytic network process at each iteration of the procedure. Fuzzy numbers could be used to represent uncertain parameters of the procedure. The proposed procedure could also be used to develop software with decision support systems methodology.