International Journal of Computational Intelligence Systems

Volume 9, Issue 1, January 2016, Pages 153 - 167

ELECTRE I Method Using Hesitant Linguistic Term Sets: An Application to Supplier Selection

Authors
Ali Fahmi1, fahmi@itu.edu.tr, Cengiz Kahraman2, kahramanc@itu.edu.tr, Ümran Bilen3, umranbilen@gmail.com
Received 7 September 2015, Accepted 18 December 2015, Available Online 1 January 2016.
DOI
10.1080/18756891.2016.1146532How to use a DOI?
Keywords
Multiple criteria decision analysis (MCDA); outranking methods; ELECTRE method; hesitant fuzzy linguistic terms set (HFLTS); ordered weighted averaging (OWA); computing with words (CWW)
Abstract

Decision making is a common process in human activities. Every person or organization needs to make decisions besides dealing with uncertainty and vagueness associated with human cognition. The theory of fuzzy logic provides a mathematical base to model the uncertainities. Hesitant fuzzy linguistic term set (HFLTS) creates an appropriate method to deal with uncertainty in decision making. Managerial decision making generally implies that decision making process conducts multiple and conflicting criteria. Multi criteria decision analysis (MCDA) is a widely applied decision making method. Outranking methods are one type of MCDA methods which facilitate the decision making process through comparing binary relations in order to rank the alternatives. Elimination et Choix Traduisant la Réalité (ELECTRE), means elimination and choice that translates reality, is an outranking method. In this paper, an extended version of ELECTRE I method using HFLTS is proposed. Finally, a real case problem is provided to illustrate the HFLTS-ELECTRE I method.

Copyright
© 2016. the authors. Co-published by Atlantis Press and Taylor & Francis
Open Access
This is an open access article under the CC BY-NC licence (http://creativecommons.org/licences/by-nc/4.0/).

1. Introduction

Multiple criteria decision analysis (MCDA) denotes to analyzing decision making situations that encounter with multiple and conflicting criteria. It refers to analyzing tools and systematic approaches that empower the decision maker (DM). They help DM to represent his/her preferences and consider all subjective and objective conditions to assess the decision elements [1, 2]. MCDA methods provide a quantitative infrastructure to model the assessments of criteria and alternatives [3]. MCDA methods are classified into three major categories including multi-attribute value theory (MAVT), multi-objective mathematical programming, and outranking methods [4, 5, 6, 7, 8, 9].

(1) MAVT approach concerns hierarchical decision making problems with the overall goal on the top level and the criteria on the lowest level. MAVT involves in the weighting the criteria and assessment of alternatives with respect to criteria. Lastly, MAVT approach ranks the criteria through calculations of criteria weights and alternative assessments.

(2) Multi-objective mathematical programming methods focus on finding the aspiration points. These points can generate non-dominated solutions by scalarizing functions through reaching down from ideal solution [10, 11, 12].

(3) Outranking methods form two main steps: building the outranking relations, and exploiting the outranking relations [13]. At first, outranking methods systematically compare criteria. The binary comparisons construct the concordance and discordance sets which are linguistic comparisons. Then, these comparisons lead to numerical concordance and discordance indices.

There are different outranking methods including ELECTRE family, PROMETHEE family, QUALIFLES [14], ORESTE [15, 16], MELCHIOR [17], PRAGMA [18], MAPPACC [18], and TACTIC [19]. ELECTRE family includes ELECTRE I [20, 21], ELECTRE II [22, 23], ELECTRE III [24, 25], ELECTRE IV [26], ELECTRE IS [27], ELECTRE TRI [28], and ELECTREGKMS [29]. PROMETHEE family contains PROMETHEE I [30], and PROMETHEE II [30].

The PROMETHEE family of outranking methods includes the PROMETHEE I for partial ranking of the alternatives and the PROMETHEE II for complete ranking of the alternatives, the PROMETHEE III for ranking based on interval, the PROMETHEE IV for complete or partial ranking of the alternatives when the set of viable solutions is continuous, the PROMETHEE V for problems with segmentation constraints, the PROMETHEE VI for the human brain representation [31].

ELECTRE family methods are the most widely used outranking approaches [29, 32]. ELECTRE I outranking method is the first ELECTRE method which was introduced by Benayoun, Roy, and Sussman [20] and Roy [21]. Decision makers use ELECTRE I to construct a partial prioritization and choose a set of promising alternatives [4]. ELECTRE II is used for ranking the alternatives based on the determination of concordance and discordance matrices for each criterion and alternative pair. In ELECTRE III, an outranking degree is established, representing an outranking creditability between two alternatives which makes this method more sophisticated. ELECTRE III is based on the principle of fuzzy logic and uses the preference and indifference thresholds while determining the concordance and discordance indices [13].

Outranking methods are extended to deal with imprecision and uncertainty of decision making process. Roy [33] and Siskos, Lochard, and Lombard [34] developed fuzzy outranking method. This method uses fuzzy concordance and discordance relation. Sevkli [35] applied fuzzy ELECTRE method for supplier selection. Ertay and Kahraman [36] evaluated the design requirements by utilizing fuzzy outranking methods.

Afterwards, different extensions of fuzzy sets have been applied to outranking methods. Interval type-2 fuzzy sets have been employed in the ELECTRE method [37, 38]. Chen [39] implemented an ELECTRE based group decision making method using interval type-2 fuzzy sets. Devi and Yadav [40] developed an ELECTRE group decision making based on intuitionistic fuzzy sets for plant location selection. See Vahdani et al. [41], Li, Lin, and Chen [42], and Wue and Chen [43] developed ELECTRE method extension in an intuitionistic fuzzy environment. Hatami-Marbini and Tavana [4] applied fuzzy group decision making to extend ELECTRE I methodology. See Hatami-Marbini and Tavana [4] for more details.

Torra [44] proposed hesitant fuzzy sets (HFS) as a generalization of fuzzy sets. HFS are quite suitable in decision making when experts have to assess a set of alternatives. HFS were introduced to the literature for the common difficulty that often appears when the membership degree of an element must be established and the difficulty is not because of an error margin (as in intuitionistic fuzzy sets) or due to some possibility distribution (as in type-2 fuzzy sets), but rather because there are some possible values that cause hesitancy about which one would be the right one [45].

These new sets permit DM to represent his/her hesitation in decision making process. The mathematical background of hesitation is modeled by different membership functions. HFS establishes a relationship between the envelopes of fuzzy sets and generates new fuzzy sets with combined membership functions.

As decision making encounters linguistic information, we need to use computing with words (CWW) processes [46, 47, 48]. HFS provides an appropriate infrastructure to develop the concept of hesitant linguistic term set which is a useful CWW process. The experts involved in decision problems under uncertainty cannot easily provide a single term as an expression of his/her knowledge such as poor, good, very good, etc. because they might consider several terms at the same time or looking for a more complex linguistic term such as at least medium poor, at most very high, etc. Thus, Rodriguez, Martinez, and Herrara [49], and Liu and Rodriguez [50] proposed the hesitant fuzzy linguistic term set (HFLTS) and its application in decision making. This method could increase the ability and flexibility of linguistic elicitation of DM. HFLTS enables DM to represent the uncertainty in his/her assessments of actions, criteria, alternatives, etc. [49]. This would be implemented by context-free grammars based on the fuzzy linguistic tools to compare the expressions.

In this study, the ELECTRE I method is extended based on HFLTS. This new approach could be applied to decision making problems and deal with the uncertainty and imprecision of multiple criteria decision making. This method contains two phases: HFLTS phase and ELECTRE I phase. Firstly, HFLTS method considers DM’s hesitation between several values to evaluate the decision elements and manage the uncertainty of decision making. These hesitant linguistic terms are expressed by numerical fuzzy values. Next, in ELECTRE phase, we apply ELECTRE I method to analyze the decision making. This would be accomplished through building the outranking relations, and then exploiting the outranking relations. Finally, we can rank the alternatives and make an appropriate decision.

It is worthy to mention that ELECTRE method has an important pitfall [4]. DM must initiate this method by precise measurement of the performance ratings and the weights of criteria [5]. On the other hand, DM may prefer to express his/her judgment by linguistic expressions [51, 52, 53]. This approach takes this weakness into account and modifies it by utilizing linguistic expressions in decision making process and allows DM to evaluate the real world problem with imprecise and vague ratings.

This paper is organized as follows: Section 2 includes needed definitions and arithmetic operations about HFLTS and ELECTRE I method. In Section 3, the step-by-step methodology of hesitant fuzzy linguistic term set ELECTRE I (HFLTS-ELECTRE I) is proposed. Next, a real case study is added in Section 4. Lastly, the conclusions and future works are presented in Section 5.

2. Preliminaries

In this section, the basic definitions and arithmetic operations are given.

2.1. Comparative Linguistic Expressions

In many decision making situations, DM deals with uncertain and imprecise information to express his/her judgment. Zadeh [54] presented fuzzy sets and then developed the concept of linguistic variables [51]. Based on these concepts, Rodriguez, Martinez, and Herrera [49] developed the concepts of fuzzy hesitant linguistic terms set (HFLTS).

Definition 1. [51]

A linguistic variable has five characteristics and is defined as (H,T(H),U,G,M), where H is the name of variable, T(H) is the term set of H, i.e., the collection of its linguistic values, U is the universe of discourse, G is a syntactic rule which generates the terms in T(H), and M is a semantic rule which associates with each linguistic value X its meaning, M(X) denotes a fuzzy subset of U.

Most of decision making situations deal with comparative judgments. DM provides relative assessments by comparing the alternatives. Therefore, in order to apply the concept of HFLTS in decision making problems, we should elicit the comparative linguistic expressions [55].

Definition 2. [49, 50]

A hesitant fuzzy linguistic term set (HFLTS) is shown by HS, and is an ordered finite subset of the consecutive linguistic terms set S = {S0,S1,…,Sg}.

Definition 3.

The envelope of an HFLTS is a linguistic interval that its upper bound and lower bound determine the limits of the envelope as follows:

env(HS)=[HS,HS+],HSHS+HS+=max{Si}SiHSHS=min{Si}SiHS

As mentioned above, DM applies comparative linguistic expressions to assess the criteria and alternatives of a common decision making problem. It is appropriate to represent the comparative expressions by fuzzy membership functions [50]. Thus, trapezoidal membership function is used and ordered weighted averaging (OWA) operator is established [56] in order to compute the elements of fuzzy membership function. OWA weights would represent the decision maker’s hesitation in linguistic terms. There are different approaches to obtain the OWA weights, but in this study Filev and Yager’s [56] approach is applied.

Definition 4.

Let OWA operator maps from dimension n as follows:

OWA:RnROWA(a1,a2,,an)=j=1nwjbj
where a1, a2, …, an is the aggregated set of arguments and bj is the jth largest argument of the aggregated set, and the associated weighting vector W = (w1,w2,…,wn)T satisfies wi ∈ [0,1] ∀i = 1,2, …,n
i=1nwi=1

In order to calculate the parameters b and c, OWA operator is implemented as follows:

b=OWAWs(aMi,aMi+1,,aMj)
c=OWAWt(aMi,aMi+1,,aMj)
where Ws and Wt represent the form of OWA weighting vectors for computing b and c, respectively. s, t = 1,2 and st or s = t.

Definition 5.

Let α be the parameter in the unit interval (0, 1). This parameter is used to calculate the first type of weights of OWA operator as follows:

W1=(w11,w21,,wn1)Tw11=α,w21=α(1α),w31=α(1α)2,,wn11=α(1α)n2,wn1=(1α)n1

Also, the second type of weights of OWA operator is as follows:

W2=(w12,w22,,wn2)Tw12=α,w22=α(1α),w32=α(1α)2,,wn12=α(1α)n2,wn2=(1α)n1

2.2. Fuzzy arithmetic operations

Based on the previous definitions, proposed method employs trapezoidal fuzzy numbers and arithmetic operations of trapezoidal fuzzy numbers should be calculated. The required concepts are defined below.

Definition 6.

A fuzzy number à = (a,b,c,d) is a trapezoidal fuzzy number shown in Fig. 1, and its membership function is given by:

μ(x)={0,x<a,xaba,axb,1,b<x<c,dxdc,cxd,0,x>d,
Fig. 1.

Trapezoidal fuzzy number, Ã = (a,b,c,d)

Definition 7.

Let A and B be two positive trapezoidal fuzzy numbers, where à = (a,b,c,d) and B˜=(p,q,r,s) . The arithmetic operations between à and B˜ is as follows:

  • Addition operation:

    A˜B˜=(a+p,b+q,c+r,d+s)

  • Subtraction operation:

    A˜B˜=(as,br,cq,dp)

  • Multiplication operation:

    A˜B˜(ap,bq,cr,ds)

  • Division operation:

    A˜B˜(as,br,cq,dp)

  • Multiplication operation:

    u.A˜=A˜.u={(u·a,u·b,u·c,u·d)ifu0(u·d,u·c,u·b,u·a)ifu<0

  • Division operation:

    A˜u={(au,bu,cu,du)ifu0(du,cu,bu,au)ifu0

2.3 Basic definitions of ELECTRE I

As mentioned earlier, ELECTRE I is an outranking method that is characterized due to analyzing binary relations of the alternatives. Some definitions regarding outranking relations are provided as follows [4, 5, 57, 58, 59]:

Definition 8.

In ELECTRE I method, DM indicates the preferences by binary outranking relations. Let x and y be two alternatives. “x S y” means “x is at least as good as y”.

Definition 9.

The concept of outranking relation is based on two concepts including the concordance and the discordance. The statement of “x S y” provides insights into these concepts: The concordance concept: For an outranking “x S y” to be validated, a sufficient majority of the criteria should be in favor of this assertion. On the other hand, when the concordance condition holds, discordance concept implies that none of the criteria in the minority should oppose too strongly to the assertion “x S y”.

3. Methodology

As mentioned above, we propose a new fuzzy ELECTRE I method based on HFLTS. This method facilitates decision making process due to utilizing linguistic terms and context-free grammars for evaluating the importance weights of the criteria and performance ratings. Since DMs compare the criteria and alternatives to present relative assessments, comparative linguistic expressions are needed for the assessments of alternatives and criteria. Therefore, HFLTS method contributes to determine fuzzy numbers for hesitant judgments. Then, fuzzy ELECTRE I outranking method is established to complete the decision making process. Our methodology includes two main phases: HFLTS operations and the application of fuzzy ELECTRE I outranking method.

I. First phase: HFLTS operations

Step 1. Determine the semantics and syntax of linguistic terms set and context-free grammar GH.

In this phase, DM should determine the decision making conditions to apply CWW process. The first step is to determine the semantic and syntax of comparative linguistic terms set. Definition 1 also lets us to determine the context-free grammar GH. An example of context-free grammar is provided below [50]:

GH={VN,VT,I,P}
where VN shows non-terminal symbols and VT shows terminal symbols as follows:
VN={primary term,binary relation,unary term,conjunction}VT={at least,at most,between,and,S0,S1,,Sg}IVN

P indicates the production rules for the context-free grammar GH. Bracket signs shows optimal elements and the symbol “|” indicates alternative elements. In this example, the production rules are as follows:

P={primary term::=s0|s1||sg,binary relation::=between,unary relation::=at least|at most|lower than|greater than,conjunction::=and}

This set specifies the fuzzy partition of the HFLTS phase. This partitioning lets us to represent the fuzzy envelopes of the assessments. Three values are considered including left value (L), middle value (M), and right value (R) for the partitioning. As you see in Fig. 2, each triangular represents three values L, M, and R. In the proposed approach, assessments are represented by trapezoidal fuzzy numbers.

Fig. 2.

Fuzzy partitioning.

Step 2. Gather the assessments of performance and criteria.

In this step, DM evaluates the performance ratings with respect to criteria in square matrices. Also, the weights of criteria are calculated by DM’s assessments of criteria with respect to goal.

Step 3. Transform the assessments into HFLTS by transformation function EGH.

Let S = {s0,s1,…,sg} be the linguistic terms set. DM evaluates the performance and criteria by comparative linguistic terms based on the defined context-free grammar GH. The transformation of the evaluations into HFLTS will be carried out by transformation function EGH as follows:

EGH(si)={si|siS}={si},EGH(betweensiandsj)={sk|sisksjandskS}={si,si+1,,sj},EGH(at leastsi)={sj|sjsiandsjS}={si,si+1,,sg},EGH(atmostsi)={sj|sjsiandsjS}={s0,s1,,si}

Step 4 : Obtain an aggregate set for each assessment.

In this step, the fuzzy envelope is obtained through calculating the elements of a trapezoidal fuzzy number. We calculate these elements by aggregate set as follows: Let A be the aggregate set for EGH (between si and sj)

Abetween={aLi,aMi,aLi+1,aRi,aMi+1,aLi+2,aRi+1,,aLj,aRj1,aMj,aRj}
Aatleast={aLi,aMi,aLi+1,aRi,aMi+1,aLi+2,aRi+1,,aLg,aRg1,aMg,aRg}
Aat most={aL0,aM0,aL1,aR0,aM1,aL2,aR1,,aLi,aRi1,aMi,aRi}

As you see in Fig. 2, aLk+1=aMk=aRk1 for k = 1,2, …,g − 1. So that, the Eqs. (15), (16), and (17) are simplified as follows:

Abetween={aLi,aMi,aMi+1,,aMj,aRj}Aatleast={aLi,aMi,aMi+1,,aMg,aRg}Aat most={aL0,aM0,aM1,,aMi,aRi}

Step 5. Obtain fuzzy envelope for each assessment.

T˜between=(a,b,c,d),
a

= minAbetween=min{aLi,aMi,aMi+1,,aMj,aRj}=aLi ,

d

= maxAbetween=max{aLi,aMi,aMi+1,,aMj,aRj}=aRj ,

b

= OWAW2(aMi,aMi+1,,aMj) ,

c

= OWAW1(aMi,aMi+1,,aMj) .

As in Definition 5, we can compute the weights of OWA operator to calculate element b, and c for “between” comparative linguistic terms set. The first type of weights of OWA operator is used to calculate b and the second type of weights of OWA operator to calculate c as follows:

  1. I)

    If i+j is odd, then

    b=OWAW2(aMi,aMi+1,,aMi+j12),c=OWAW1(aMj,aMj1,,aMi+j+12),W2=(w12,w22,,wji+1/22)Tandα1=g(ji)g1
    where
    w12=α1ji12,w22=(1α1)α1ji32,,wji122=(1α1)α1,wji+122=(1α1)b=w12·aMi+j12,w22·aMi+j32,,wji122·aMi+1,wji+122·aMi
    and W1=(w11,w21,,wji+1/21)T and α2=1α1=(ji)1g1

    where

    w11=α2,w21=α2(1α2),,wji122=α2(1α2)ji32,wji+122=(1α2)ji12c=w11·aMi+j+12,w21·aMi+j12,,wji121·aMi+1,wji+121·aMi

  2. II)

    If i+j is even, then

    b=OWAW2(aMi,aMi+1,,aMi+j2),c=OWAW1(aMj,aMj1,,aMi+j2),W2=(w12,w22,,wji+2/22)T
    where
    w12=α1ji2,w22=(1α1)α1ji22,,wji22b=w12·aMi+j2,w22·aMi+j22,,wji22·aMi+1,wji+222·aMi
    and W1=(w11,w21,,wji+2/21)T where
    w11=α2,w21=α2(1α2),,wji21=α2(1α2)ji22,wji+221=(1α2)ji2c=w11·aMi+j2,w21·aMi+j22,,wji21·aMj1,wji+221·aMj

Finally, the HFLTS envelope will be formed as T˜between=(aLi,b,c,aRj) . See Fig. 3.

Fig. 3.

Fuzzy Envelope for comparative linguistic expression “between”

T˜at least=(a,b,c,d),
a

= minAatleast=min{aLi,aMi,aMi+1,,aMg,aRg}=aLi ,

d

= maxAat least=max{aLi,aMi,aMi+1,,aMg,aRg}=aRg

b

= OWAW2(aMi,aMi+1,,aMg)

c

= OWAW2(aMi,aMi+1,,aMg)

As can be seen in Definition 5, the weights of OWA operator to compute element b for “at least” comparative linguistic terms set are as follows:

n=gi+1,W2=(w12,w22,,wn2)Tandα=ig
where
w12=αgi,w22=(1α)αgi1,,wgi2=α(1α)wgi+12=1αb=w12·aMg,w22·aMg1,,wgi2·aMi+1,wgi+12·aMi,b=αgi·aMg,(1α)αgi1·aMg1,,α(1α)·aMi+1,(1α).aMic=aMg

Finally, the HFLTS envelope is formed as T˜at least=(aLi,b,aMg,aLg) . See Fig. 4.

Fig. 4.

Fuzzy Envelope for comparative linguistic expression “at least”.

T˜at most=(a,b,c,d),
a

= minAat most=min{aL0,aM0,aM1,,aMj,aRj}=aL0 ,

d

= maxAat most=max{aL0,aM0,aM1,,aMj,aRj}=aRj

b

= OWAW1(aM0,aM1,,aMj)

c

= OWAW1(aM0,aM1,,aMj)

As you see in Definition 5, the weights of OWA operator to compute element c for “at most” comparative linguistic terms set are as follows:

n=i+1,W1=(w12,w22,,wn2)Tandα=ig
where
w11=α,w21=α(1α),,wi1=α(1α)i1,wi+11=(1α)ic=w11·aMi,w21·aMi1,,wi1·aM1,wi+11·aM0,c=α·aMi,α(1α)·aMi1,,α(1α)i1·aM1,(1α)i·aM0b=aM0

Finally, the HFLTS envelope is formed as T˜at most=(aL0,aM0,c,aRi) . See Fig. 5.

Fig. 5.

Fuzzy Envelope for comparative linguistic expression “at most”

II. Second phase: ELECTRE I outranking method

Step 6. Defuzzify and rank the alternatives for each criterion.

In this step, the trapezoidal fuzzy values, which DM evaluated the alternatives in previous steps, will be defuzzified. The center of gravity defuzzification method is applied to rank alternatives with respect to each criterion.

Then, the utility matrix Ũ is built. By taking utility matrix into account, assessment matrices will be formed to compare alternatives with respect to each criterion. The trapezoidal fuzzy elements of final utility matrix form the performance ratings as in Eq. (18).

U˜=[t˜12t˜1nt˜21t˜2nt˜m1t˜m2]

The decision matrix Ĩ is also built to compare the importance of criteria. Next, the weight of each criterion could be obtained through computing the average of the elements in each row using summation and division rules of trapezoidal fuzzy numbers. See Definition 7 for more details. Finally, the weight matrix will be obtained as Eq. (20).

I˜=[F˜12F˜1nF˜21F˜2nF˜m1F˜m2]
W˜=[W˜1,W˜2,,W˜m]T
where W˜1=j=2nF˜1jn,W˜2=j=2nF˜2jn , and W˜m=j=2nF˜mjn

Since the values in the matrices are between 0 and 1, there is no need to obtain a normalized matrix. The weight matrix W˜ and the utility matrix Ũ are multiplied to construct weighted utility matrix as follows:

V˜=[v˜11v˜12v˜1nv˜21v˜22v˜2nv˜m1v˜m2v˜mn]
where
v˜ij

= (vijl,vijp,vijq,viju)=(wjltijl,wjptijp,wjqtijq,wjutiju) ,

i = 1,2, …, m, and j = 1,2, …, n.

Step 7. Calculate the fuzzy concordance indices and discordance indices.

The fuzzy concordance matrix is built as Eq. (21). The elements of this matrix are measured by pairwise comparison of alternatives. In this regard, two alternatives and are considered. Using af and ag are considered. Using Definition 8, the outranking relation agSagf will be defined. The arguments in favor of this statement “af outranks ag” is measured to obtain the fuzzy concordance indices as follows:

C˜=[c˜g1c˜1(m1)c˜1mc˜g1c˜gfc˜g(m1)c˜gmc˜m1c˜mfc˜m(m1)]c˜gf=(cgfl,cgfp,cgfq,cgfu)
where cgfl=k:xikxjkwklW,cgfp=k:xikxjkwkpW,cgfq=k:xikxjkwkqW , and cgfu=k:xikxjkwkuW and W=i=1mw˜i .

In addition, the discordance set is formed to construct the discordance matrix. The concept of discordance set is in contrast to the concept of concordance set. If DM has doubt upon the statement “af outranks ag” and DM is against the assertion “ag is at least as good as af”, we can represent the discordance matrix and the arguments as follows [60]:

D=[dg1d1(m1)d1mdg1dgfdg(m1)dgmdm1dmfdm(m1)]

where dgf=h°|v˜ghv˜fh|h|v˜ghv˜fh| , if v˜gh<v˜fh , g = 1,2, …, m, and f = 1,2, …, n, and |v˜ghv˜fh|=|(vghlvfhu,vghpvfhq,vghqvfhp,vghuvfhl)| .

The discordance indices are obtained between alternatives i and j with respect to the criterion h°. ih and jh are ranked by Yager’s centroid index [61]. Then, the defuzzification method, center of gravity will be employed to calculate the discordance index.

Since this method calculates the outranking relations without defining concordance and discordance thresholds, a combination method to unite fuzzy concordance and discordance matrices is proposed. A modified version of Aouam and Chang’s [62] formula is applied to calculate the elements of matrix D′. Hence, D′ is formed by subtracting each element of the discordance matrix from 1. Then, the fuzzy global matrix Z˜ will be calculated through peer to peer multiplication of the elements of the matrices C˜ and D′ (using Eq. (11)) as follows:

Z˜=C˜D
where
D=[1dg11d1(m1)1d1m1dg11dgf1dg(m1)1dgm1dm11dmf1dm(m1)]

Step 8. Compare the alternatives by fuzzy Kernel diagram.

In the final step, the outranking relations is exploited to compare the alternatives by the elements of matrix Z˜ . In this regard, graph G = (V,J) is defined, where V is the set of nodes and J is the set of arcs. One node is considered for each alternative and an arc between two nodes of alternatives af and ag. Here, we propose fuzzy Kernel diagram. This diagram could depict uncertainty in comparison of alternatives. We draw solid arc if ag absolutely outranks af, that means ag is absolutely preferred to af as shown in Fig. 6 a. We proposed to draw long dashed arc if ag strongly outranks af as shown in Fig. 6 b, and dashed arc if ag weakly outranks af as shown in Fig. 6 c. If ag and af are incomparable, we consider no arc between them as shown in Fig. 6 d. In addition, Fig. 6 e illustrates the indifference between ag and af.

Fig. 6.

Outranking relations represented by fuzzy Kernel diagram.

In order to determine the arcs, the elements of fuzzy global matrix Z˜ should be defuzzified. In this step, center of gravity defuzzification method is applied and zgf represents the defuzzified values. The intervals of the defuzzified values are defined for associated binary preferences are presented below:

{Incomparable0zgf0.25Weak Preference0.25<zgf0.5Strong Preference0.5<zgf0.75Absolute Preference0.75<zgf1

In summary, you can find the associated steps of our combined method as follows:

  1. I.

    First phase: HFLTS operations

    • Step 1.

      Determine the semantics and syntax of linguistic terms set and context-free grammar GH.

    • Step 2.

      Gather the assessments of performance and criteria.

    • Step 3.

      Transform the assessments into HFLTS by transformation function EGH.

    • Step 4 :

      Obtain aggregate set for each assessment.

    • Step 5.

      Obtain fuzzy envelope for each assessment.

  2. II.

    Second phase: ELECTRE I outranking method

    • Step 6.

      Rank alternatives for each criterion.

    • Step 7.

      Calculate fuzzy concordance and discordance indices.

    • Step 8.

      Compare alternatives by fuzzy Kernel diagram.

4. An Illustrative Example

In this part, a simple MCDA problem is added to provide a step by step solution in order to illustrate the proposed method. The problem is the selection of best supplier among three alternatives by considering three attributes including price, quality, and delivery. The best way of illustrating the relations between the criteria and alternatives is to use a hierarchy often applied in Analytic Hierarchy Process (AHP) [63]. The hierarchy of decision making is provided below:

Fig. 7.

The hierarchy of MCDA supplier selection problem.

The steps of the method are presented below:

Step 1. At first, the semantic and syntax of linguistic term set S and context-free grammar GH should be defined. We determine the linguistic term set S as:

S={no importance(ni),very low importance(vli),low importance(li),mediumimportance(mi),high importance(hi),very high importance(vhi),absolute importance(ai)}

Step 2. In this step, the importance of criteria is assessed with respect to the goal as Table 1. This will obtain the weights of criteria in next steps. Then, the performance of alternatives with respect to each criterion will be evaluated as you see in Tables 2, 3, and 4.

Price Quality Delivery
Price - At most li At most mi
Quality At least hi - At least hi
Delivery At least mi At most li -
Table 1

The linguistic assessment of criteria with respect to the goal.

Supplier 1 Supplier 2 Supplier 3
Supplier 1 - At most vli At most li
Supplier 2 At least vhi - At least hi
Supplier 3 At least hi At most li -
Table 2

The linguistic assessment of alternatives with respect to price.

Supplier 1 Supplier 2 Supplier 3
Supplier 1 - At most li At most li
Supplier 2 At least hi - Between li and mi
Supplier 3 At least hi Between mi and hi -
Table 3

The linguistic assessment of alternatives with respect to quality

Supplier 1 Supplier 2 Supplier 3
Supplier 1 - At most li At most mi
Supplier 2 At least hi - At least hi
Supplier 3 At least mi At most li -
Table 4

The linguistic assessment of alternatives with respect to delivery.

Steps 3, 4, and 5. Firstly, each of comparative linguistic expressions is transformed into corresponding HFLTS. Then, the aggregate set is formed for each assessment. This aggregate set lets us to obtain fuzzy envelopes. For more details, follow the provided example below:

EGH(at most li)={ni,vli,li}

By using Eq. (17), simplified aggregate set is obtained:

Aat most li={aLni,aMni,aMvli,aMli,aRli}

We form fuzzy envelope similar to the presented fuzzy partitioning in Fig. 5.

a=minAat most li=min{aLni,aMni,aMvli,aMli,aRli}=aLni=0,d=maxAat most li=max{aLni,aMni,aMvli,aMli,aRli}=aRli=0.5,b=OWAW1(aLni,aMni,aMvli,aMli,aRli),c=OWAW1(aLni,aMni,aMvli,aMli,aRli).n=i+1=2+1=3,α=ig=26=13W1=(w12,w22,,wn2)T
where
w11=1/3,w21=1/3(11/3),w31=(11/3)2,W1=(13,29,49)T,c=w11·aMli,w21·aMvli,wi+11·aMni,c=13·aMli+29·aMvli+49·aMni=13·0.33+29·0.17+49·0=0.15b=aMni=0T˜at most li=(0,0,0.15,0.5)

Eventually, a trapezoidal fuzzy number will be obtained for each of evaluations as shown in Tables 5, 6, 7, and 8.

Price Quality Delivery
Price - (0, 0, 0.15, 0.5) (0, 0, 0.35, 0.67)
Quality (0.5, 0.85, 1, 1) - (0.5, 0.85, 1,1)
Delivery (0.33, 0.65, 1, 1) (0, 0, 0.15, 0.5) -
Table 5

The numerical assessments of the criteria with respect to the goal.

Supplier 1 Supplier 2 Supplier 3
Supplier 1 - (0, 0, 0.028, 0.33) (0, 0, 0.15, 0.5)
Supplier 2 (0.67, 0.972, 1, 1) - (0.5, 0.85, 1, 1)
Supplier 3 (0.5, 0.85, 1, 1) (0, 0, 0.15, 0.5) -
Table 6

The numerical assessments of the alternatives with respect to price.

Supplier 1 Supplier 2 Supplier 3
Supplier 1 - (0, 0, 0.15, 0.5) (0, 0, 0.35, 0.67)
Supplier 2 (0.5, 0.85, 1, 1) - (0.17, 0.33, 0.5, 0.67)
Supplier 3 (0.5, 0.85, 1, 1) (0.33, 0.5, 0.67, 0.83) -
Table 7

The numerical assessments of the alternatives with respect to quality.

Supplier 1 Supplier 2 Supplier 3
Supplier 1 - (0, 0, 0.15, 0.5) (0, 0, 0.35, 0.67)
Supplier 2 (0.5, 0.85, 1, 1) - (0.5, 0.85, 1, 1)
Supplier 3 (0.33, 0.65, 1, 1) (0, 0, 0.15, 0.5) -
Table 8

The numerical assessments of the alternatives with respect to delivery.

As mentioned before, the weights of criteria should be calculated using assessment of criteria with respect to goal as shown in Table 5. Similarly, the average performance of three suppliers is obtained by calculating the average amount of each row in Tables 6, 7, and 8. Thus, the utility matrix, which represents the average performance ratings, and the weights of criteria are shown in Table 9. The weighted utility matrix is given in Table10.

Price Quality Delivery
Weights (0, 0, 0.25, 0.585) (0.5, 0.85, 1, 1) (0.165, 0.325, 0.575, 0.75)

Supplier 1 (0, 0, 0.89, 0.415) (0, 0, 0.15, 0.5) (0, 0, 0.25, 0.585)
Supplier 2 (0.585, 0.911, 1, 1) (0.335, 0.59, 0.75, 0.835) (0.5, 0.85, 1, 1)
Supplier 3 (0.25, 0.425, 0.575, 0.75) (0.415,0.675,0.835,0.915) (0.165,0.325,0.575,0.75)
Table 9

The utility matrix (average performance ratings) and criteria weights.

Price Quality Delivery
Supplier 1 (0, 0, 0.222, 0.242) (0, 0, 0.15, 0.5) (0, 0, 0.143, 0.438)
Supplier 2 (0, 0, 0.25, 0.585) (0.167,0.501,0.75,0.835) (0.082,0.276,0.575,0.75)
Supplier 3 (0, 0, 0.143, 0.438) (0.207, 0.573, 0.835,0.915) (0.027,0.105,0.330,0.562)
Table 10

The weighted utility matrix.

Step 6. In this step, the alternatives are ranked with respect to each criterion as presented in Tables 11, 12, and 13. We apply center of gravity defuzzification method to rank the fuzzy numbers. Considered fuzzy numbers are average amount of row elements of performance ratings shown in Tables. 6, 7, and 8. This ranking will help us to obtain fuzzy concordance and discordance indices in next step. For instance, center of gravity defuzzifies (0, 0, 0.089, 0.415) as follows:

00.089xdx+0.0890.415(3.07x2+1.27x)dx00.089dx+0.0890.415(3.07x+1.27)dx=0.143

Average Defuzzified Value
Supplier 1 (0, 0, 0.089, 0.415) 0.143
Supplier 2 (0.585, 0.911, 1, 1) 0.855
Supplier 3 (0.25, 0.425, 0.575, 0.75) 0.488

Ranking Supplier 2 ≻ Supplier 3 ≻ Supplier 1
Table 11

Ranking of alternatives with respect to price.

Average Defuzzified Value
Supplier 1 (0, 0, 0.15, 0.5) 0.178
Supplier 2 (0.335, 0.59, 0.75, 0.835) 0.597
Supplier 3 (0.415, 0.675, 0.835, 0.915) 0.700

Ranking Supplier 3 ≻ Supplier 2 ≻ Supplier 1
Table 12

Ranking of alternatives with respect to quality.

Average Defuzzified Value
Supplier 1 (0, 0, 0.25, 0.585) 0.220
Supplier 2 (0.5, 0.85, 1, 1) 0.820
Supplier 3 (0.165, 0.325, 0.575, 0.75) 0.480

Ranking Supplier 2 ≻ Supplier 3 ≻ Supplier 1
Table 13

Ranking of alternatives with respect to delivery.

Step 7. In this step, fuzzy concordance and discordance indices are computed to construct fuzzy concordance and discordance matrices. By taking into account the ranking of alternatives, the fuzzy concordance indices are obtained through adding criteria weights each of which the alternative ranking is higher. Then, it will be divided by the total amount of criteria weights. Fuzzy concordance table is presented in Table 14:

Supplier 1 Supplier 2 Supplier 3
Supplier 1 - (0,0,0,0)w˜com.+w˜accl.+w˜des. (0,0,0,0)w˜com.+w˜accl.+w˜des.
Supplier 2 w˜com.+w˜accl.+w˜des.w˜com.+w˜accl.+w˜des. - w˜com.+w˜des.w˜com.+w˜accl.+w˜des.
Supplier 3 w˜com.+w˜accl.+w˜des.w˜com.+w˜accl.+w˜des. w˜accl.w˜com.+w˜accl.+w˜des. -
Table 14

Fuzzy concordance table.

From Table 14, we have

C˜=(0,0,0,0)(0,0,0,0)(1,1,1,1)(0.25,0.27,0.45,0.57)(1,1,1,1)(0.43,0.55,0.73,0.75)

Based on the Eq. (22), the discordance matrix is obtained using Table 10 as follows:

D=[1100.23100.768]

Based on Eq. (24), the supplement of discordance matrix will be as follows:

D=[0010.76910.232]

Finally, the fuzzy global matrix Z˜ is built through multiplication of matrices C˜ and D′ as follows:

Z˜=[(0,0,0,0)(0,0,0,0)(1,1,1,1)(0.192,0.207,0.346,0.438)(1,1,1,1)(0.099,0.127,0.169,0.174)]

After defuzzification, the preferences of the alternatives will be obtained by using Eq. (25). Therefore, supplier 2 is absolutely preferred to supplier 1, supplier 3 is absolutely preferred to supplier 1, and supplier 2 is weakly preferred to supplier 3. Supplier 3 is also incomparable to supplier 2.

Step 8. Consequently, the fuzzy Kernel graph is drawn based on the Z˜ matrix as Fig. 8. According to Fig. 6, solid arc is drawn from suppliers 2 and 3 to supplier 1 to represent absolutely preferred nodes, and dashed line from supplier 2 to supplier 3. Eventually, the fuzzy Kernel graph in Fig. 9 depicts that supplier 2 outranks other two suppliers using HFLTS-ELECTRE I method.

Fig. 8.

Tat most li

Fig. 9.

Fuzzy Kernel graph

5. Conclusion

In this study, ELECTRE I outranking method using hesitant fuzzy linguistic term set is proposed. Since decision making involves uncertainty and imprecision, DMs prefer to render linguistic judgments about decision alternatives and criteria. Our proposed method facilitates decision making process due to using comparative linguistic expressions to evaluate the alternatives and criteria and it could be a DM-friendly MCDA method. Fuzzy sets theory and subsequent developments empowers us to numerically interpret the uncertainty of MCDA problems. Ultimately, ELECTRE I method is established to investigate binary relations of alternatives and criteria in order to compare them by comparative linguistic terms.

In this proposed HFLTS-ELECTRE I method, the combination matrix of concordance and discordance remains fuzzy. Our method offers an evolution in outranking relations through elimination of concordance and discordance thresholds form ELECTRE I method. Hence, we propose fuzzy Kernel diagram to represent the preferences of actions. This new diagram uses dashed lines to depict decision maker’s uncertainty toward actions. Since fuzzy Kernel diagram is implemented, defuzzification of the fuzzy concordance matrix is not needed and we continue the calculations by fuzzy values until the final step. Thus, the output clearly illustrates decision maker’s assessments and his/her preferences. Future works could focus on extensions of HFLTS for other members of ELECTRE family or outranking methods. In addition, extension of fuzzy Kernel diagram could represent decision maker’s preferences in detail.

References

[3]RT Clemen, Making hard decisions: an introduction to decision analysis, 2nd ed, Duxbury Press at Wadsworth Publishing Company, Belmont, 1996.
[5]J Figueira, S Greco, and M Ehrgott, Multiple criteria decision analysis: state of the art surveys, Springer, New York, 2005.
[17]JP Leclercq, Propositions d’extension de la notion de dominance en presence de relations d’ordre sur les pseudo-criteres: MELCHIOR, Revue Belge de Recherche Operationelle de Statistique et d’hformatique, Vol. 24, No. 1, 1984, pp. 32-46.
[20]R Benayoun, B Roy, and B Sussman, ELECTRE: Une méthode pour guider le choix en présence de points de vue multiples, in Note de travail 49, SEMA-METRA international, direction scientifique (Paris, 1966).
[22]B Roy and P Bertier, La Methode ELECTRE II: Use Methode de Classement en Presence de Criteres Multiples., Note de Travail No. 142, Direction Scientifique, Group Metra, Paris, 1971.
[24]B Roy, Electre III Un algorithme de classements fonde sur une representation floue en presence de criteres multiples, Cahairs du CERO, Vol. 20, No. 1, 1978, pp. 3-24.
[27]B Roy and JC Hugonnard, Classement des prolongements de lignes de metro en banlieue parisienne (presentation d’une methode multicritere originale), Cahiers du CERO, Vol. 24, No. 2,3,4, 1982, pp. 153-171.
[28]W Yu, ELECTRE TRI: Aspects methodologiques et manuel d’utilisation, Document du LAMSADE 74, Universite Paris-Dauphine, 1992.
[32]B Roy and D Bouyssou, Aide Multicritere a la Decision: Methodes et Cas, Economica, Paris, 1993.
[34]JL Siskos, J Lochard, and J Lombardo, A multicriteria decision-making methodology under fuzziness: application to the evaluation of radiological protection in nuclear power plants, TIMS Studies in the Management Sciences, Vol. 20, 1984, pp. 261-283.
[58]DJ Dubois, Fuzzy sets and systems: theory and applications, Academic Press, New York, 1980.
[60]KP Yoon and CL Hwang, Multiple Attribute Decision Making, An Introduction, Sage University Paper series on Quantitative Applications in the Social Sciences, Sage, Thousand Oaks, 1995.
[63]TL Saaty, The Analytic Hierarchy Process, McGraw Hill, New York, 1980.
Journal
International Journal of Computational Intelligence Systems
Volume-Issue
9 - 1
Pages
153 - 167
Publication Date
2016/01/01
ISSN (Online)
1875-6883
ISSN (Print)
1875-6891
DOI
10.1080/18756891.2016.1146532How to use a DOI?
Copyright
© 2016. the authors. Co-published by Atlantis Press and Taylor & Francis
Open Access
This is an open access article under the CC BY-NC licence (http://creativecommons.org/licences/by-nc/4.0/).

Cite this article

TY  - JOUR
AU  - Ali Fahmi
AU  - Cengiz Kahraman
AU  - Ümran Bilen
PY  - 2016
DA  - 2016/01/01
TI  - ELECTRE I Method Using Hesitant Linguistic Term Sets: An Application to Supplier Selection
JO  - International Journal of Computational Intelligence Systems
SP  - 153
EP  - 167
VL  - 9
IS  - 1
SN  - 1875-6883
UR  - https://doi.org/10.1080/18756891.2016.1146532
DO  - 10.1080/18756891.2016.1146532
ID  - Fahmi2016
ER  -