Multiple Criteria Decision-Making in Heterogeneous Groups of Management Experts

In commercial organizations operations, frequently some dynamic events occur which involve operational, managerial, and valuable information aspects. Then, in order to make a sound decision, the business professional could be supported by a Multi Criteria Decision-Making (MCDM) system for taking an external course of action, as, for instance, forecasting a new market or product, up to an inner decision concerning for instance, the volume of manufacture. Thus, managers need, in a collective manner, to analyze the actual problems, to evaluate various options according to diverse criteria, and finally choose the best solution from a set of various alternatives. Throughout these processes, uncertainty and hesitancy easily arise, when it comes to define and judge criteria or alternatives. Several approaches have been introduced to allow Decision Makers (DMs) to deal with. The Interval Multiplicative Preference Relations (IMPRs) approach is a useful technique and the basis of our proposed methodology to provide reliable consistent and in consensus IMPRs. In this manner, DMs’ choices are implicitly including their uncertainty while maintaining both an acceptable individual consistency, as well as group consensus levels. The present method is based on some recent results and an optimization algorithm to derive reliable consistent and in consensus IMPRs. In order to illustrate our results and compare them with other methodologies, a few examples are addressed and solved.


Introduction
Among the vast world of commercial operations organizations and despite their differences, they have common business operations or activities such as acquiring inventory, hiring employees and cashing from customers.Nowadays, inside each modern organization, we can frequently find an information system working in synchrony with these business operations.Furthermore, several important information systems are fed by operating departments (work centers of the organization), and, as a result, these systems' outcomes can be used to manage these operations.As a consequence, managers analyze their corresponding information system in light of the work that the organization performs.For instance, when a marketing manager is required to advise management and to have some reports for management decision-making, she/he must understand the organization's product cycles.
There are various events which occur while organizations engage their business operations as, for instance, diverse trends in purchases and sales.This dynamic data coming from these events are frequently recorded and kept up in a database to mirror and supervise business operations.These records include operational, managerial, and valuable information details.Thus, in order to design and use a group decision-making system, the business professional must previously consider what kind of event data is needed and the necessary process to extract the useful information.
In essence, the type of decision under consideration rules what kind, amount and quality of the information must be used to make a sound decision.Furthermore, information is more valuable when it recognizes and doesn't disregard the personal management styles, the main choices of the the decision maker, the weighting of each decision maker with respect to the entire decision-making process, and managers' uncertainty and hesitancy when a heterogenous group tries to achieve an integral decision.
Multi criteria decision-making then could support human-centered management in taking a course of action for a sound decision-making to elucidate several management problems.For instance, what products to sell and the suitable targeted market to sell them, or the structure of the organization better suited for this process or even the direction and motivation of employees based on some known standards.Another type of management problem can be an inner decision as for instance, human resources information, volume of manufacture and the available delivering chains, which is useful to develop alternative methods for manufacturing and delivering a new product.
On the other hand, the essence of the information provided to managers must be according to the management level.For example, strategic level managers need information allowing them to assess the environment to forecast future events and conditions, where s/he may not be as concerned with the timeliness or accuracy of the information as her/his interest is in the trends.In this case, most of the information is exogenous to the organization, or tactical management needs information mainly coming from pertinent operational units.Some exogenous information is needed, as well as more detailed and accurate information than the information required at the strategic level.Finally, operational management unit usually requires exogenous narrower information but more detailed and accurate information.It comes largely from within the organization where frequent decisions are made, with shorter lead times to respond in a timely manner to current variations as, for instance, in sales patterns.
As we noted in this section, managers support the entire firms' making decision to achieve diverse kind of goals.They have basically different kinds of perspectives and local goals to reach a high efficiency level at their respective departments.However, the main goal of the firm needs, at certain moments, to combine local goals and perspectives of the work centers, with particular global targets.
In this scenario of heterogeneous managers group, an MCDM system for supporting the management for group decision-making becomes paramount.
Fortunately, there are a vast variety of MCDM techniques ranging from: analytic hierarchy process (AHP) [1], multi attribute utility theory (MAUT) [2], simple multi attribute rating technique (SMART) [3], fuzzy set theory (FST) [4], data envelopment analysis (DEA) [5], case-based reasoning (CBR) [6], simple additive weighting (SAW) [7], elimination et choice translating reality (ELECTRE) [8], technique for order of preference by similarity to ideal solution (TOPSIS) [9], preference ranking and organization method for enrichment evaluation (PROMETHEE) [10], and goal programming (GP) [11].Another interesting method for addressing uncertainty is the Hesitant Fuzzy Sets (HFS) [12] where some recent contributions for heterogeneous information are found [13][14][15][16][17][18].These methods are appropriate for uncertainty problems, since fuzzy logic aims to represent human preferences based on individual opinions expressed through a linguistic setting.Its major drawback is that membership functions (which can be seen as intervals) are fixed and also, until now, there does not exist detailed work related to an appropriate analytical tool (norms, aggregation operators, etc.) for a particular study case.Thus, it results in a greater uncertainty and vagueness to the problem.
Group Decision Making (GDM) is a main MCDM issue, where multiple DMs (managers in our case) act collaboratively and collectively, analyze decision-making problems, evaluate goals according to a set of criteria, and finally choose the best solution from a set of alternatives [36].As noted in [37], when organizations gather specialized groups or larger groups, and the number of alternatives increase, unanimity may be difficult to attain.For this reason, flexible or milder benchmarks (definitions) of Group Consensus and Individual Consistency have been employed.
Group Consensus (GC) has to do with group cooperation and agreement since the alternative, option, or goal to be achieved is the best course of action for the whole organization.On the other hand, Individual Consistency (IC) concerns each DM to have her/his information, and, consequently, her/his judgments, free of contradictions.
Derived of the blend of heterogenous DM group, a common problem in big organizations is that managers often can accurately state definitions and assessments on the priority rate of the set of criteria or alternatives, when it concerns their own operating departments.However, when it comes to define and judge any other set of criteria or alternatives, they face uncertainty and hesitation problems.
In this paper, we address and solve these problems by allowing the DM to utilize blended assessments.For example, DMs can use crisp values for her/his ratio judgments of criteria or alternatives when s/he is confident; and then they can use intervals where these are used to express her/his uncertainty and hesitant assessments.
The aim of this paper is to synthesize a novel approach in order to provide a reliable measure of both the IC and the GC of a set of these blended assessments, which we called Interval-Multiplicative Preference Relations (I-MPRs).Then, in the next step, the approach is verified by a constrained optimization algorithm where the improved I-MPRs will finally fulfill the Individual Consistency and Group Consensus Indices since both restrictions are involved in the constraints.
In order to verify the requirements of acceptable IC and GC levels, we use the Hadamard's operator.As soon as IC and then GC are validated, an Interval Priority Vector can be obtained from this set of ordered judgment I-MPRs, in order to rank the alternatives as a final result of the DM analysis.Various techniques are revised and the prioritization method which indicates the entire order of the intervals and the preference degrees, is addressed according to our results.
The main advantages of our approach follow: • It is provided through a couple of algorithms and a nonlinear optimization approach concurrently applied.

•
Through the Hadamard's operator and some easy algebraic manipulations, objective functionals are synthesized (as it will be detailed further on), to be used in the optimization algorithm.

•
When the I-MPRs improved by the methodology are reduced into an MPR (defined in the I-MPR), our approach can still give reliable results.For example, for this MPR, we can verify the results of IC or GC, with an alternative method.

•
The IC or the GC accepted indices (threshold values) have been previously investigated and fixed.Nevertheless, the project designer could assign a different value depending on the project requirements.

•
Obtained results are independent of the method of prioritization utilized in the consensus operation.
This paper is organized as follows: In Section 2, some preliminaries (Definitions, Theorems and Lemmata) are given to support a basis of the main methodologies and techniques previously described above and the approach introduced here.In Section 3, an extension of analysis and results derived in the former section is used in the I-MPR framework for obtaining our main results.Then, in Section 4, a slight modification of a prioritization method is given in light of the results.In Section 5, some numerical examples are solved and compared with other methodologies.Finally, in Section 6, concluding remarks are given about the main advantages, drawbacks and future research of our methodology.

Preliminaries
Let us consider a Group Decision-Making problem and let D = {d 1 , d 2 , • • • , d m } be the set of DMs, and C = {c 1 , c 2 , • • • , c n } be a finite set of criteria (or alternatives), where c i denotes the ith criteria.
In the AHP framework [1], a pairwise comparison matrix or an MPR is given by a DM where s/he provides judgments through a ratio (c i /c j ) for every pair of criteria (or alternatives) (c i and c j ) to represent the preference degree of the first criteria (or alternative) over the second.A Saaty's Scale is frequently used to pick up a value for this preference ratio where SS = [1/9 9].
Thus, an MPR for instance A = (a ij ) n×n , is a positive reciprocal n × n matrix, a ij > 0, such that Let λ = {λ 1 , • • • , λ m } be the weight vector of the m − th DM, where λ s > 0, s ∈ M, ∑ m s=1 λ s = 1, which can be derived with several techniques (see for instance [38] and the references cited therein).
An MPR n × n matrix is called a completely consistent MPR (cf.[1]) if Thus, their corresponding completely consistent MPR K = (k ij ) n×n can be constructed from each MPR n × n matrix as follows: Furthermore, let A c be the group MPR (cf.[39]) which represents the group opinion utilizing the geometric average operator: Let us denote (cf.[40]), an I-MPR A t given by the t − th expert as where ji and (5)

Measuring the Dissimilarity between Matrices
A useful operator to measure the degree of dissimilarity between two MPRs is the Hadamard Product (HP).The HP of A = (a ij ) n×n and B = (b ij ) n×n is defined by Consequently, the degree of dissimilarity between A and B is defined as d(A, B) = 1 n 2 e T A • B T e or: 1 where e = (1, In order to have an assessment of Individual Consistency (CI), one can measure the compatibility of A l with respect to (w.r.t.) its own completely consistent matrix K given by Equation (2).Thus, where In a similar manner, the group consensus index of each MPR, i.e., GCI A c (A l ), l = 1, 2, • • • , m, is based on the compatibility of A l w.r.t. the group opinion given by A c in Equation (3).Thus, the assessment of group consensus for each MPR A l given by GCI A c (A l ) read as: where the index is usually set at GCI = 1.1 and A l is an individual MPR.
From Equations ( 7) and ( 8), for CI K , it follows Respectively from Equations ( 7) and ( 9), for GCI A c evaluated for an A l follows: An MPR A t is completely consistent if and only if CI K (A t ) = 1.Thus, a threshold useful to measure the similarity of two MPRs up to an acceptable level of consistency was suggested by [41,42] as CI = 1.1.Similarly, an MPR A t is completely in consensus if and only if GCI A c (A t ) = 1 and an acceptable level of consensus is GCI = 1.1.
In the following up to the end of the section, we utilize Definitions and Theorems recently introduced to measure the Individual Consistency Index and Group Consensus Index of a set of I-MPRs (cf.[43,44]).For details and proofs, please refer to these works. .
Definition 2. The Individual Consistency Index of the I-MPRs (A t ) n×n given by Equation ( 4) when one has generated the set of MPRs { o A t } n×n given by Definition 1, is defined by Definition 3. The smallest Individual Consistency Index of the I-MPRs (A t ) n×n given by Equation ( 4) when one has generated the set of MPR { o At}n×n given by Definition 1, is defined by From Definition 2, the next theorem is set forth.given by the I-MPRs (A t ) n×n , then In order to illustrate how useful are these Definitions and Theorem, we apply them on two I-MPRs given by two experts to assessing three criteria.

Example 1.
Let us consider two experts evaluating a set of three criteria through: where for A 1 it implies w 1 w 2 w 3 , and for A 2 one gets w 1 w 2 w 3 .Note: Experts coincide in their rankings of the three criteria.Nevertheless, the assessment ratios are different.
In the following, the Definitions 1-3 are figured out for these two I-MPRs.
From Definition 1 applied to A 1 : In a similar manner, from Definition 1 applied to A 2 : where, from now on, symbol * corresponds to the inverted entries, respectively.
From Definition 2 applied to A 1 and A 2 , it read as: The first expert has the higher Individual Consistency Index which means that s/he has a weaker consistency in her/his judgments.It coincides with our logical expectations by analyzing and comparing both I-MPRs.
By defining for A 1 and A 2 an MPR within their own intervals, for instance and by applying Theorem 1 to both MPR, one has: In order to reach a Group Consensus solution, the experts are supposed to participate in rounds of discussion to meet a common goal.Then, definitions are next introduced to compute the dissimilarity amongst DM' assessments for finally obtaining the GCI level.
In a similar manner as Definition 1, one can compute the whole combination of the values given by the set of I-MPRs through: Definition 4. From the m I-MPRs given as in Equation ( 4), (A t ) n×n , t = 1, 2, • • • , m, their set of MPRs { Ǎ} n×n which is the whole combination of the interval values of the set of (A t ) n×n I-MPRs with , is given by Furthermore, a set of { Ǎc } n×n which represents each group opinion utilizing the geometric average operator corresponding to each element of the set { Ǎ} n×n , can be defined.Definition 5. Based on Definition 4, one can compute for the set of I-MPRs (A t ) n×n and { Ǎ} n×n , t = 1, 2, • • • , m, its set of MPR { Ǎc } n×n .Then, the corresponding set of { Ǎc } n×n which is the Group I-MPRs Opinion utilizing the geometric operator is given by where { ǎc ij } n×n is the combination of the whole set of values ( ).
Definition 6.The Group Consensus Index of a set of I-MPRs (A t ) n×n , t = 1, 2, • • • , m, given by Equation ( 4) Ǎc }, given by Definition 5, is defined by and Ȟt is associated with the corresponding p Ǎc matrix and p ∈ {1, where ν = 2  16) and ( 17), it follows: From Definition 5 applied to Equation ( 25), one has: where, from Equation (3), it implies: From Definition 6, it results: These results indicate that the first expert A 1 , which has the lowest weight (λ 1 = 2/7), is not in consensus.
In addition, from Definition 7, it follows: As soon as we obtain ICI and GCI of a set of I-MPRs, the accepted threshold values (ICI and GCI respectively) give how acceptable in Consistency and Consensus the assessments of the DMs are.Whenever some I-MPRs are over their index values, we can apply the following strategy to improve their Individual Consistency and/or their Group Consensus Indices.

Reliable Intervals for Individual Consistency and Group Consensus
Once the ICI has been calculated, it may happen that one or more I-MPRs are not consistent.For example, r I-MPRs were not consistent.In that case, we need to improve the Individual Consistency of those I-MPRs and then let us define this set of MPRs as: Naturally, each one of the A o I , I = 1, 2, • • • , r doesn't verify for a special combination of its values (given by Equation ( 12)), the inequality given in Equation (10).
In a similar manner, once the GCI has been calculated, it may happen that one or more I-MPRs are not in consensus.For example, for instance y I-MPRs were not in consensus.In that case, we need to improve the Group Consensus of those I-MPRs and then let us define this set of MPRs as Once again, each one of the A P J , J = 1, 2, • • • , y doesn't verify, for a special combination of values (given by Equation ( 23)), the inequality given in Equation (11).

Sequential Quadratic Programming Methodology
When one has a constrained nonlinear optimization problem (NLP) and wants a numerical solution, the Sequential Quadratic Programming (SQP) [45] is a useful approach.Let us consider the implementation of the sequential quadratic programming methodology to NLP of the form: where f : R n −→ R is the objective functional, the functions h : R n −→ R m and g : R n −→ R p describe the equality and inequality constraints.The sequential quadratic programming is an iterative method for modelling the NLP for a given iterate x k , k ∈ N 0 , by a Quadratic Programming (QP) subproblem.
As soon as the method solves that QP subproblem, the solution is used to build a new iterate x k+1 .This is figured out in such a way that the sequence (x k ) k∈N 0 converges to a local minimum x * of the NLP Equation ( 32) as k → ∞.
Then, as we can observe, the problem of minimizing an objective functional (given in our case by Equation ( 7)), where one has some inequalities constraints from Equation ( 10), (related to g(x)) with the additional observation of constraints imposed by Equation ( 30) (related to h(x)), defines an NLP to improve the Individual Consistency Index of an I-MPR.
In a similar manner, the problem of minimizing an objective functional (given by Equation ( 7)) providing some inequalities constraints from Equation (11), (related to g(x)) with the additional observation of constraints imposed by Equation (31) (related to h(x)), defines an NLP to improve the Group Consensus Index of a set of I-MPRs.
An additional observation is that since we have interval judgments in I-MPRs, we should apply the SQP to minimize the objective functional to found the minimal values verifying the required constraints.After that, we should then apply the SQP to maximize the objective functional to found the maximal values.In this manner, we will finally find reliable intervals where are fulfilled inequalities and equalities constraints and conditions of Equations ( 10) and (11), to obtain acceptable Individual Consistency and Group Consensus, respectively.

Matching the Problem with the SQP for Improving I-MPRs
From Equation (30) (Equation (31) resp.), the following inequalities can be stated in terms of optimization variables x i , as follows: where δ stands for o in the Individual Consistency assessments and for p in the Group Consensus assessments.For the Individual Consistency, the inequalities given from Equation (10) are 1 which impose some inequality constraints to be fulfilled.For Group Consensus improvement, they are given from Equation (11) as 1.0 ≤ d( ǍJ , ȞJ ) ≤ 1.1.
The relationship between and the set of I-MPRs (A t ) n×n , t = 1, 2, • • • , m under analysis is given by: x where ij stands for the corresponding crisp or interval value of the I-MPR.From the definition of constrained nonlinear optimization problem, we note that when the I-MPR is in fact just an MPR (when the respective expert has a high confidence in her/his assessments), the SQP could not give any different solution.This is because the optimization variables will remain unchanged since the variables x i , i = 1, 2, 3, • • • , (n 2 − n)/2 will be defined as a crisp values (cf.Equation ( 33)).
Thus, in order to provide a general benchmark to address and solve any possible case of assessments given through I-MPRs, we will modify slightly the SQP.We introduce a design parameter which will be used as an additive or subtractive element, for high and low bounds, respectively, of the I-MPRs to be improved.
For example, we modify Equation (33) as follows: where > 0. For the sake of compactness, let us integrate in the same terms the parameter .
For example, a In that manner, when the SQP does not provide a feasible solution, we will again iterate with an incremental -value, until a solution is found.
It is worth mentioning that must be initialized, with a small value, to keep the maximum of the information provided by the expert.

Individual Consistency Objective Functional
Thus, when an I-MPR is inconsistent, we need to synthesize the objective functional.For example, from Equations ( 7) and (10), one obtains: After some algebraic manipulations, it implies:

Group Consensus Objective Functional
In a similar manner, when an I-MPR (A o J ) n×n is not in consensus, we need to synthesize the objective functional.For example, from Equation (11), one obtains: After some algebraic manipulations, it implies for (A J ) n×n : where α 1 equals to the term over the first brace, α 2 equals to the term over the second brace, and so on.Furthermore, note that the term a excepts the term a , the term a (1)(λ 1 ) 13 • a (2)(λ 2 ) 13 excepts the term a , and so on.
Finally, the initialization point x 0 used in the NLP can be set at any point within the corresponding interval given by Equation (33).
In the following, an algorithm based on an NLP is used to obtain reliable intervals for the assessment of based decision models such as those given by the set of I-MPRs for both indices (CI and GCI).Two similar methods are blended to obtain both reliable I-MPR Consistency and Consensus levels.For example, for Individual Consistency improvement, we can find an I-MPR verifying Individual Consistency.For the other case, for Group Consensus improvement, a reliable Consensus Index for I-MPRs is obtained.In order to use an NLP, an SQP algorithm can be found in [45].In the following, our algorithm is described in detail.

Improving the Individual Consistency of an I-MPR
A scheme of the Algorithm 1 implementation is depicted in Figure 1.Once the set of I-MPRs are all acceptably consistent, we can improve the Group Consensus level through the following algorithm.

Algorithm 1: Algorithm IC-I-MPR
ij ) n×n : the initial interval I-MPR; x 0 : the initialization value for the nonlinear optimization, which is to be defined within the corresponding interval or as the corresponding crisp value; the threshold value of CI for the Individual Consistency assessment; c : Design parameter value allowing the enlargement pace of the searching space of the algorithm.Output: A o I : the consistency interval matrix computed and verifying interval conditions given by Equation (10).
Step 1: Get the function for the assessment of Individual Consistency given by Equation (37).
Step 2: Define for the nonlinear optimization algorithm, the allowed intervals: Thus, assign the former linear inequality constraints as follows: Step 3: Obtain the acceptable index of Individual Consistency, 1.0 ≤ d( * AI, * K) ≤ CI.Thus, based on matrix K given by Equation ( 2), the nonlinear inequality constraints imposed are given by Equation (10).
Step 4: Solve the former nonlinear optimization problem using the SQP algorithm to minimize it.
Step 6: Solve again the same nonlinear optimization problem but this time in order to maximize it.In order to do so, assign the objective functional as − f (x). where , stands for the interval or crisp value obtained.

Start Decision Makers Board
Decision Targets and Framework Information Sources I-MPR:

Improving the Group Consensus of a Set of I-MPRs
A scheme of the Algorithm 2 implementation is depicted in Figure 2.

Algorithm 2: Algorithm GC-I-MPR
ij ) n×n : the initial interval I-MPRs; x 0 : the initialization value for the nonlinear optimization, which is to be defined within the corresponding interval or as the corresponding crisp value; GCI for the Group consensus assessment; g : design parameter allowing the enlargement of the searching space of the algorithm.
Output: A P J : the I-MPRs computed and verifying interval conditions given by Equation (11).
Step 1: Get the function for the assessment of Group Consensus given by Equation (39).
Step 2: Define for the nonlinear optimization algorithm, the allowed intervals: Thus, assign the former linear inequality constraints as follows: Step 3: Obtain the acceptable Group Consensus Index d( ǍJ , ȞJ ).Thus, based on the set of { Ǎc } n×n given by Equation ( 22), the nonlinear inequality constraints imposed is given by 1.0 ≤ d( ǍJ , ȞJ ) ≤ 1.1.
Step 4: Otherwise, continue with the next step.
Step 5: Solve the former nonlinear optimization problem (NLP) using the SQP algorithm to minimize it.
Step Step 9: Compose the J − th Group Consensus Interval Matrix A P J as follows: where , stands for the interval or crisp value obtained.
Section Remarks: -In the case that an expert has provided a crisp value(s) in her/his judgement(s), this value(s) drives the process of nonlinear optimization since they will slightly change with the pace of .It is very useful since precisely in that value(s), the expert has shown her/his highest confidence level.-At the end of both algorithms, one gets reliable I-MPRs, i.e., where the consistency and consensus constraints are fulfilled.

Prioritization Method and Methodology Application
The process of deriving a priority weight vector of the alternatives, w i = (w 1 , w 2 , • • • , w n ) T from an I-MPR, for instance an A t , is called a prioritization method, where w l ≥ 0 and ∑ n l=1 w l = 1.Then, when one has the evaluation for Ǎ which is associated to the maximum value (cf.Equation ( 23)) of the Group Consensus Index of I-MPRs A t , it follows from Equation ( 9): where w i ( Ǎ) and w i ( Ȟt ) are the weights of Ǎ and Ȟt , respectively.In a similar manner, for A which is associated to the minimum value (cf.Equation ( 13)) of the Group Consensus Index of I-MPRs A t , one has: where w i (A) and w i (H t ) are the weights of A and H t , respectively.Consequently: As soon as the set of I-MPRs is acceptably consistent and in consensus, the interval priority vector to rank alternatives is obtained.To do so, an interval ranking is used.

Interval Priority Vector Synthesis
An interval priority vector should reflect different expert's risk preferences for her/his interval judgments.There are mainly two prioritization methods: (among others cf.[46] and the references therein) Eigenvalue-based Methods (EM), (cf.[1,47]) and the Row Geometric Mean Method (RGMM), (cf.[48]) that are utilized to derive a priority weight vector from an ordered judgment matrix which is the method here utilized.
Based on the results given in [49], a slight modification of their method is addressed below.
where if P(a > b) > 0.5, then a > b; if P(a > b) = 0.5, then a = b and if P(a > b) < 0.5, then a < b; and in this manner Equation (49) indicates the total order of intervals and their preference degrees [50].
In order to provide access to the public as on-line security and resource conditions allow, these complete methodology (definitions, applied theorems, algorithms and interval priority vector synthesis) will be soon available in a site.Then, this benchmark will be useful to test a submitted set of I-MPRs for assessment on their Individual Consistency (ICI) and the Group Consensus Indices (GCI).This site will test different data on various based decision models.
In the following, our main methodology is applied to numerical examples under diverse considerations in order to test different situations when the system is working in a real scenario.

Illustration of the Methodology through Numerical Examples
In the following example, let us suppose that three managers from different operating departments are participating in a group decision about the evaluation of a new market to have the best return of investment of a set of products.There are four decision criteria which involve the details of judgments, which are: development, legal restrictions for economic activities, society and infrastructure denoted as C 1 , C 2 , C 3 , and C 4 , respectively.The three mangers state their preferences over the four criteria and, since they are confident on the analysis carried out by their operating department, they have a crisp value for criteria judgments related to their department.Nevertheless, they provide an interval for others criteria assessments.Their associated weight vector is λ = [1/3, 1/3, 1/3].

Example 3. The three managers provide their I-MPRs through
where the initial ranking given by expert 1 is w and by expert 3 is: w 3 0.8179 > w 4 .Note that the first manager has more information and expertise on the third criteria.On the other hand, the second manager is focused on the first criteria and finally the third manager is focused on the second criteria.
Let us now apply Definition 1 to A 1 , A 2 and A 3 as follows: where Ar}, r = 1, 2, 3 are obtained in a similar manner as in Example 1.Then, CI * K ( * Ai), i = 1, 2, 3 is calculated through Equation ( 2) where, for the first expert, one obtains: and its corresponding consistent MPR (calculated by Equation ( 2) Similar calculations are carried out for A 2 and A 3 .Then, from Definition 3, CI K ( Ãi ), i = 1, 2, 3 is calculated through Equation ( 13) where for the first expert one obtains: and its corresponding consistent MPR (calculated by Equation ( 2)) is Similar calculations are carried out for A 2 and A 3 .
From Equation ( 10 From Equation (30) Let us provide the solution for A 1 which has the highest CI.Then, from Equation ( 33) and the first expert, the following inequalities will be utilized in the optimization process: where c = 0.001, but for the first iteration c = 0.0.
From Equation ( 37) and variable's definition given by Equation (34), the objective functional is obtained by replacing a 12 = x 1 , a 13 = x 2 , a 14 = x 3 , a 23 = x 4 , a 24 = x 5 and a 34 = x 6 .Thus, one obtains: where β 1 equals the term over the first brace, β 2 equals the term over the second brace, and so on.In addition, sign = 1 is defined to obtain the minimization and sign = −1 for the maximization.The same procedure is used for the second expert (A 2 ).After applying the SQP, the results follow: Now that we have the set of I-MPRs in acceptable Individual Consistency, we can proceed with solving the Group Consensus for A 1 , A 2 and A 3 .
Let us now apply Definition 4 to these I-MPRs.The set { 1 Ǎ} is calculated as in Equation (25).For example, (59) From Definition 5, let us apply it to A 1 , A 2 and A 3 where, by Equation (3), one obtains: where (a Since we have detailed the calculations of the Algorithm 1, in the next example, we provide details of the Algorithm 2.
In the next example, we address and solve a problem of partner selection that is determined for the formation of a virtual enterprise already considered in [51,52].
A virtual enterprise is a dynamic association of enterprises, working together to benefit from a market opportunity by giving a solution that could not be delivered individually.Naturally, it comes with a market opportunity and then it is stopped after that market particular event.Thus, a virtual enterprise is, by definition, a non-permanent alliance of diverse, autonomous, and in some cases geographically distributed organizations sharing resources and skills having common objectives and profiting from a benefit window in the market opportunities [51].
Example 4. In this example, the main enterprise needs to select a partner to grasp a new market opportunity where four candidates (alternatives) must be analyzed by four DMs.The CEO's management staff (DMs) is involved in the partner selection process which should select the most suitable partner.Then, each DM compares each pair of alternatives c i and c j , and gives her/his preference assessment through the next I-MPRs: where the initial individual order of intervals for each DM is given by x 4 ; DM3:x 2 1 x 1 1 x 4 1 x 3 & DM4:x 3 .6382 x 1 .9904 x 2 1 x 4 .In [52], it is found that the third expert is more consistent and then it is assigned a higher weight.They obtained the final ranking of the alternatives as: x 2 x 1 x 4 x 3 .On the other hand, when all experts are considered equally weighted, they obtained the final ranking of the alternatives as x 1 x 2 x 3 x 4 .
Let us apply our method to check for consistency and consensus indices.By applying the first part of our Definitions, the Individual Consistency of each DM (by using Equation ( 12)) is given by: where it is noted that the third DM (A 3 ) has provided a slightly inconsistent I-MPR.Since the third expert's judgments A 3 are not consistent, let us apply the Algorithm IC-I-MPR.
From Equation (31) A p 1 = A 4 and from Equation (33) the following inequalities will be utilized in the optimization process: where g = 0.001, but for the first iteration g = 0.0.From Equation ( 39) and variable's definition given by Equation ( 34), the objective functional is obtained by replacing in it a 13 , a 14 , a 23 , a 24 and a (4) 34 .Then, one obtains: where α 1 equals the term over the first brace, α 2 equals the term over the second brace, and so on.In addition, sign = 1 is defined to obtain the minimization and sign = −1 for the maximization.
Let us then apply the SQP algorithm for A ij ) with GCI = 1.1 and = 0.001.As soon as the optimization algorithm converges, the new Group Consensus Indices from Equation ( 23) is read as: GCI Ǎ1 = 1.0132,GCI Ǎ2 = 1.0921,GCI Ǎ3 = 1.0508,GCI Ǎ4 = 1.0697. (67) In addition, from Equation (24): The priority vector is w Second case: For an equal experts' weighting (For example, λ = [1/4, 1/4, 1/4, 1/4]), one has that the Group Consensus Indices for each decision maker (by using Equation ( 23)) are: where the third and fourth DM have provided their I-MPRs not in consensus.Then, apply again Algorithm GC-I-MPR for A As soon as the optimization algorithm converges, the new Group Consensus Indices from Equation ( 23) is read as: In addition, from Equation (24): where the set of I-MPRs are Individually Consistent and in an acceptable Group Consensus.They are read as: We note that, through Definitions, Theorem and Algorithms introduced here, we obtain the same results of [51,52] when the experts' weight is the same.However, when the weighting of the experts is different, some small differences are found.The above is due to the variation in the weighting provided here and in those articles.

Case Study Discussions and Managerial Implications
In the first example, the second and third DMs have provided inconsistent I-MPRs, and it was sufficient to improve their Individual Consistency to have the three I-MPRs in Group Consensus.Note that the ranking order of each DM was initially slightly different among them; nevertheless, at the end of the process, the Group Consensus points to an acceptable and analytic decision.
As we have have seen along the first example, when a group of DMs needs to state a point of agreement, they can define their assessments through I-MPRs which do not necessarily need to be exactly in the same judgment direction.Sometimes, they only need to be confident in the quality of their decisions based on an MCDM system which analyzes them through a well-known decision model.
For the second example, only two DMs evaluate in the same direction the set of criteria (alternatives).The other two experts are even in clear contradiction, since the first DM states that the third criteria is lesser in importance and the second DM states that this criteria is the most important.Furthermore, by assigning to the third expert a highest evaluation confidence (weight), the methodology produced a very interesting result, by selecting a global ranking order that none of them had chosen.
On the other hand, when every DM has the same weight, the methodology produced also a different global ranking order that none of them had chosen.The DMs A 1 , A 2 and A 3 , have preserved their individual ranking order at the end of the process, and only the fourth DM has had his individual ranking order changed.This result can be used in his/her operating department as an internal feedback to reconsider their position with respect to the other departments and the organization objectives.
As soon as a DM obtains the results, s/he can use them so that, with this information, he can state new ways of organizing her/his operating department.For example, this method can also be carried out within the operating departments since each one of their business operations have criteria and alternatives that can be better emphasized.
In summary, the MCDM methodology can provide results which could reinforce the position of a set of DMs or point out to a new direction.

Concluding Remarks and Future Work
In this paper, we provide a methodology based on a couple of algorithms and a nonlinear optimization approach to be used when a heterogeneous managers group needs to solve an MCDM problem.
Our approach can use Interval Multiplicative Preference Relations or Multiplicative Preference Relations, and demonstrates the utilization of the methodology to synthesize reliable intervals where consistency and consensus constraints hold.Once decision makers have proposed their I-MPRs, our method can solve for these I-MPRs from well-known decision support models.One advantage of our algorithm is that DM can re-express their preferences within an interval where, usually, have to observe some constraints based on decision targets, framework rules and advice.When the DM is confident on the pair of criteria (o alternatives) under evaluation, s/he can utilize a crisp value.On other other hand, when s/he is hesitant or uncertain about the assessment, s/he could use an interval.Our algorithm can solve independently of the used approach.
In this work, reliable I-MPRs provide a distinct advantage in interpretation of hesitancy and uncertainty about the final consistency and consensus.
Main advantages of the present approach: • It is provided through a couple of algorithms and a nonlinear optimization approach (Sequential Quadratic Programming) concurrently applied.

•
Through the Hadamard's operator and some easy algebraic manipulations, objective functionals were synthesized to be used in the optimization algorithm.

•
When the I-MPRs improved by the methodology are reduced into an MPR (defined in the I-MPR), our approach can still give reliable results.For example, for this MPR, we can verify the results of IC or GC with an alternative method.

•
The IC or the GC accepted indices (threshold values) have been previously investigated and fixed.Nevertheless, the project designer could assign a different value depending on the project requirements.

•
Obtained results are independent of the method of prioritization utilized in the consensus operation.
Main drawbacks of the present approach: • The computational cost increases as the I-MPRs dimension and the number of DMs involved in the evaluation process are increased.

•
For a real project where a high number of criteria and experts participate, it can be necessary to program this method through an exhaustive parallel computation system.

•
For a real project where a high number of criteria and experts participate, the notation can be cumbersome.
Future works aim to make an implementation on: • The application of our approach to various study cases where heterogenous groups of DMs with different weights participate in a collaborative manner.

•
The integration of the complete methodology in a benchmark to compare the results of a diverse set of MCDM tools.

•
The definition or employment of this methodology on different frameworks, v.gr.fuzzy or hesitant MCDM.

Theorem 1 .
Let (A t ) n×n = (a (t) ij ) n×n be an I-MPR given by Equation (4) and generate the set of µ MPRs { o At}n×n by using the Definition 1.If one has an MPR, (A x ) n×n = (a (x) ij ) n×n within the intervals − At + At

Step 7 :
Compose the Consistency Interval Matrix A o I as follows:

Figure 1 .
Figure 1.Process flowchart for improving consistency of an I-MPR.

7 :ij ) n×n . Step 8 :
Obtain ij ) n×n .Solve again the same nonlinear optimization problem but this time in order to maximize it.Goto to Step 2.

Figure 2 .
Figure 2. Process flowchart for improving consensus of I-MPRs.

Let us consider for instance two intervals a =b
> 0 and a, b are positioned on the x and y axis, respectively.A uniform probability distribution is assumed on the constrained area composed of a and b.For upper left points on y = x, the y values are larger than x values and viceversa for the lower right points (cf.Figure3where one possible case is shown).Then, the preference degree of P(a > b) is equal to S 1 d(a)d(b) , where d(a) ranking interval method can be stated.Definition 8.

1 Figure 3 .
Figure 3.One case for the interval analysis.