You are viewing a javascript disabled version of the site. Please enable Javascript for this site to function properly.
Go to headerGo to navigationGo to searchGo to contentsGo to footer
In content section. Select this link to jump to navigation

A 2-order additive fuzzy measure identification method based on intuitionistic fuzzy sets and its application in credit evaluation

Abstract

In view of the problem that it is difficult to quantitatively assess the interactivity between attributes in the identification process of 2-order additive fuzzy measure, this work uses the intuitionistic fuzzy sets (IFSs) to describe and deal with the interactivity between attributes. Firstly, the interactivity between attributes is defined by the supermodular game theory. On this basis, the experts employ the intuitionistic fuzzy number (IFN) to assess the interactivity between attributes, Secondly, the opinions of all experts are aggregated by using the intuitionistic fuzzy weighted average operator (IFWA). Finally, based on the aggregated results, the intuitionistic fuzzy interaction degree between attributes is defined and calculated by the score function of IFN. Thus, a 2-order additive fuzzy measure identification method based on IFSs is further proposed. Based on the proposed method, using the Choquet fuzzy integral as nonlinear integration operator, a multi-attribute decision making (MADM) process is presented. Taking the credit evaluation of the big data listed companies in China as an application example, the feasibility and effectiveness of the proposed method is verified by the analysis results of application example.

1Introduction

Due to the interaction among attributes, such as complementarity and repeatability, the classical weighted arithmetic mean method is often invalid in the process of multi-attribute decision making (MADM) [1]. In order to overcome this limitation, and effectively describe and deal with any interaction among attributes, Sugeno [2] proposed the concept of fuzzy measure and fuzzy integral. Since then, after nearly 20 years of development, the theory of fuzzy measures had been formed [3, 4]. However, in the practical application process, when there are n attributes, the general fuzzy measure generally needs to determine 2n - 2 parameters [4, 5]. This complexity greatly limits its practical application ability.

In order to reduce the complexity, Sugeno [2] proposed the λ fuzzy measure, that only needs to determine n parameters, however, the λ fuzzy measure can only represent a certain kind of interaction among attributes, which weakens its performance ability. Thereafter, in order to face with the complexity of discrete fuzzy measures, Grabisch [6] proposed the k-order additive fuzzy measures, including usual additive measures and fuzzy measures. Every discrete fuzzy measure is a k-order additive fuzzy measure for a unique k. The k-order additive fuzzy measures cover all fuzzy measures with any complexity from classical additive measure (k = 1) to general fuzzy measure (k = n). Among them, the 2-order additive fuzzy measure not only merely needs to determine n (n + 1)/2 parameters, but also merely involves the relative importance of attributes and the interaction between attributes, which well solves the contradiction between complexity and performance ability, so it has been widely used [7].

However, because the interaction between attributes is difficult to explain and understand [8], combined with the uncertainty of decision-makers’ cognition, the decision-makers often cannot give accurate quantitative assessment on the interactivity between attributes, but with a certain fuzziness and hesitation, which brings difficulties to decision-making. The existing 2-order additive fuzzy measure identification methods use the subjective methods [9–15] and the objective methods [7, 16, 17] to describe and deal with the interactivity between attributes. The subjective methods are more explanatory than the objective methods, so they have been applied more widely. However, the above methods cannot reflect the fuzziness and hesitation of decision-makers in the assessment of interactivity between attributes.

In [18], Atanassov proposed the concept of intuitionistic fuzzy set (IFS). IFS is an extension of traditional fuzzy set. It considers the information of membership degree, non-membership degree and hesitation degree, so it is more flexible and practical than traditional fuzzy set in dealing with fuzziness and uncertainty [19]. IFS can better and more completely reflect the fuzziness and hesitation of the decision-maker’s judgment of the objective things, so it has been successfully applied in the field of uncertain MADM [20–22].

Based on this observation, in view of the problem that it is difficult to quantitatively assess the interactivity between attributes, the present work uses the IFSs to describe and deal with the interactivity between attributes. As a result, a 2-order additive fuzzy measure identification method based on IFSs is proposed. Based on the proposed method, using the Choquet fuzzy integral as nonlinear integration operator [23], a MADM process is further presented.

This paper is structured as follows: Section 2 introduces the preparatory knowledge employed in this work. Section 3 proposes the 2-order additive fuzzy measure identification method based on IFSs. Section 4 presents the MADM process based on the proposed method. Section 5 describes the analysis results of application example. Section 6 discusses the results obtained and Section 7 concludes this paper.

2Preparatory knowledge

This section introduces the related definitions of IFSs, 2-order additive fuzzy measure and Choquet fuzzy integral reported in the literatures. This is the basis of Section 3 and Section 4.

Definition 1 [18]. Let X be a non-empty set, then A ={ 〈 x, μA (x) , νA (x) 〉 |x ∈ X  } is called the IFS, where μA (x) and νA (x) are the membership degree and non-membership degree of element x in X belonging to A respectively, then μA : X → [0, 1], νA : X → [0, 1] holds, and satisfies 0 ⩽ μA (x) + νA (x) ⩽1, ∀x ∈ X. In addition, πA (x) =1 - μA (x) - νA (x) is called the hesitation degree of element x in X belonging to A, and satisfies 0 ⩽ πA (x) ⩽1, ∀x ∈ X. For convenience, let α = (μα, να) be an IFN, where μα ∈ [0, 1], να ∈ [0, 1], 0 ⩽ μα + να ⩽ 1, and let Θ be the set of all IFNs.

Definition 2 [18]. Let a1 = (μ1, ν1) and a2 = (μ2, ν2) be two IFNs on a given domain; let λ be a real number, and λ > 0, then the operational rules of IFNs are as follows:

1) a1 ⊕ a2 = (μ1 + μ2 - μ1μ2, ν1ν2); 2) a1 ⊗ a2 = (μ1μ2, ν1 + ν2 - ν1ν2); 3) λa1=(1-(1-μ1)λ,ν1λ) ; 4) a1λ=(μ1λ,1-(1-ν1)λ) .

Definition 3 [24]. Let αj = (μαj, ναj) (j = 1, 2,  ⋯ , n) be a set of IFNs, wj is the weight of αj, where wj ∈ [0, 1], and j=1nwj=1 . Let IFWA : Θn → Θ, if

(1)
IFWAw(α1,α2,,αn)=(1-j=1n(1-μαj)wj,j=1nναjwj)

then IFWA is called the intuitionistic fuzzy weighted average operator.

Definition 4 [25]. Let U ={ x1, x2, ⋯ , xn } be a given domain, A=i=1nμA(xi),νA(xi)/xi is an IFS, then the intuitionistic fuzzy entropy of A is

(2)
E(A)=1ni=1nmin(μA(xi),νA(xi))+πA(xi)max(μA(xi),νA(xi))+πA(xi)

Definition 5 [26, 27]. Let α = (μα, να) be an IFN, where μα ∈ [0, 1], να ∈ [0, 1], 0 ⩽ μα + να ⩽ 1, when the hesitation degree is not considered, the score function of IFN α can be expressed as [26]

(3)
S(α)=μα-να,S(α)[-1,1]

When considering the hesitation degree, the score function of IFN α can be expressed as [27]

(4)
SL(α)=μα+μα(1-μα-να),SL(α)[0,1]

Definition 6 [2, 6]. Let X ={ x1, x2, ⋯ , xn } be a set of attributes, let X* ={ 1, 2, ⋯ , n } be a set of subscripts of attributes. P (X) is the power set of X. If the set function g : P (X) → [0, 1] satisfies the following two conditions:

1) g (∅) =0, g (X) =1; 2) If K ∈ P (X), T ∈ P (X), K ⊆ T, then g (K) ⩽ g (T).

then g is called a fuzzy measure on P (X) [2].

Grabisch [6] proposed the k-order additive fuzzy measures based on pseudo-Boolean function and Möbius transformation. On this basis, the 2-order additive fuzzy measure is then defined as

(5)
g(K)=iK*mi+{i,j}K*mij,KX
where mi is the Möbius transformation coefficient of xi (i = 1, 2, ⋯ , n), which is an overall importance; mij is the Möbius transformation coefficient of {xi, xj} (i, j = 1, 2, ⋯ , n; i ≠ j), which represents the extent of interaction between xi and xj.

Definition 7 [13]. Let X ={ x1, x2, ⋯ , xn } be a set of attributes, W ={ w1, w2, ⋯ , wn } is the weight set of X, the Möbius transformation coefficients of xi and {xi, xj} are respectively

(6)
{mi=wi/Pmij=ξijwiwj/P,i,j=1,2,,n;ij
where P = ∑iX*wi + ∑{i,j}⊂X*ξijwiwj is the sum of the importance of all xi and {xi, xj}, ξij is the interaction degree between xi and xj, ξij ∈ [-1, 1].

Definition 8 [23]. Let f be a nonnegative function defined on X, F is a σ- algebra composed of subsets of X (when X is finite, F is the power set P (X) of X), g is a fuzzy measure defined on F, then the Choquet fuzzy integral of function f on set X for fuzzy measure g is defined as (c)fdg=0g(Fα)dα , where Fα ={ x|f (x) ⩾ α, x ∈ X  }, α ∈ [0, ∞]; 0g(Fα)dα is the Riemann integral.

When X is a finite set, the elements in X are rearranged as {x(1), x(2), ⋯ , x(n) }, which makes f (x(1)) ⩽ f (x(2)) ⩽ ⋯ ⩽ f (x(n)). Let H = (c) ∫fdg, then the Choquet fuzzy integral has the following simplified formula:

(7)
H=(c)fdg=i=1n[f(x(i))-f(x(i-1))]g(X(i))
where X(i) ={ x(i), x(i+1), ⋯ , x(n) }, (i) = (1) , (2) ,  ⋯ , (n); f (x(0)) =0.

3A 2-order additive fuzzy measure identification method based on IFSs

Based on IFSs, this section proposes a 2-order additive fuzzy measure identification method. In addition, the correctness of the proposed method is proved theoretically.

Let A ={ A1, A2, ⋯ , Am } be a finite set of alternatives, and C ={ C1, C2, ⋯ , Cn } be a set of attributes to compare the alternatives, let C* ={ 1, 2, ⋯ , n } be a set of subscripts of attributes. The weight vector of attributes is WC = (w1, w2, ⋯ , wn), where wi ∈ [0, 1], and i=1nwi=1 . Let D ={ D1, D2, ⋯ , Dp } be a set of experts, the weight vector of experts is WD = (w1, w2, ⋯ , wp), where wl ∈ [0, 1], and l=1pwl=1 . Based on IFSs, the identification process of 2-order additive fuzzy measure is as follows:

Step 1: Defines the interactivity between attributes.

According to the supermodular game theory, the interactivity between attributes is defined as follows:

Definition 9. Let any two attributes Ci and Cj (ij) in attribute set C have partial order relation, the supremum Ci ∨ Cj and the infimum Ci ∧ Cj are in C, then C is called a sub-lattice [28]. Let f be a real-valued function defined on the sub-lattice C, C ⊆ Rn. For ∀Ci, Cj ∈ C, when f (Ci ∨ Cj) + f (Ci ∧ Cj) > f (Ci) + f (Cj), f is a supermodular function [28], then it is said that there is complementarity between Ci and Cj (ij); when f (Ci ∨ Cj) + f (Ci ∧ Cj) < f (Ci) + f (Cj), f is a submodular function [28], then it is said that there is repeatability between Ci and Cj (ij); particularly, when f (Ci ∨ Cj) + f (Ci ∧ Cj) = f (Ci) + f (Cj), then it is said that there is independence between Ci and Cj (ij).

Step 2: Calculate the individual assessment result of each expert on the interactivity between Ci and Cj (ij).

According to Definition 9, the expert Dl (l = 1, 2, ⋯ , p) employs the IFN α = (μα, να) to assess the interactivity between Ci and Cj (ij) ( Cn2 pairs in total). When the expert thinks that there is complementarity between Ci and Cj (ij), μα > 0.5, and the larger μα is, the stronger complementarity is. When the expert thinks that there is repeatability between Ci and Cj (ij), να > 0.5, and the larger να is, the stronger repeatability is. In addition, when the expert thinks that Ci and Cj (ij) are independent of each other, μα = να = 0.5. Thus, the individual assessment result ξ˜ijl=(μξ˜ijl,νξ˜ijl) (l = 1, 2, ⋯ , p) of the interactivity between Ci and Cj (ij) is then given.

Step 3: Calculate the group assessment result of p experts on the interactivity between Ci and Cj (ij).

Given the weight vector WD of experts, using Equation (1) to aggregate the opinions of p experts, the group assessment result ξ˜ij of the interactivity between Ci and Cj (ij) is further calculated as

(8)
ξ˜ij=(μξ˜ij,νξ˜ij)=(1-l=1p(1-μξ˜ijl)wl,l=1pνξ˜ijlwl)

Step 4: Determine the intuitionistic fuzzy interaction degree S(ξ˜ij) between Ci and Cj (ij).

With Equation (3), the intuitionistic fuzzy interaction degree S(ξ˜ij) between Ci and Cj (ij) is defined as follows:

Definition 10. Let ξ˜ij=(μξ˜ij,νξ˜ij) be the group assessment result of the interactivity between Ci and Cj (ij), according to Equation (3), if S(ξ˜ij) satisfies

(9)
S(ξ˜ij)=μξ˜ij-νξ˜ij,S(ξ˜ij)[-1,1]
then S(ξ˜ij) is called the intuitionistic fuzzy interaction degree between Ci and Cj (ij). If there is complementarity between Ci and Cj (ij), then S(ξ˜ij)>0 , and the larger S(ξ˜ij) is, the stronger complementarity is. If there is repeatability between Ci and Cj (ij), then S(ξ˜ij)<0 , and the smaller S(ξ˜ij) is, the stronger repeatability is. If Ci and Cj (ij) are independent of each other, then S(ξ˜ij)=0 .

Thus, according to Definition 10, the intuitionistic fuzzy interaction degree S(ξ˜ij) between Ci and Cj (ij) can be determined.

Step 5: Calculate the Möbius transformation coefficients mi and mij of attributes.

According to the weight wi of the i-th attribute and the intuitionistic fuzzy interaction degree S(ξ˜ij) between Ci and Cj (ij), using Equation (6), the Möbius transformation coefficients mi and mij of attributes can be calculated as

(10)
{mi=wi/Pmij=S(ξ˜ij)wiwj/P,i,j=1,2,,n;ij
where P=iC*wi+{i,j}C*S(ξ˜ij)wiwj is the sum of the importance of all Ci and {Ci, Cj} (ij).

Step 6: Identify the 2-order additive fuzzy measure gK.

According to the Möbius transformation coefficients mi and mij of attributes, using Equation (5), the 2-order additive fuzzy measure gK can be identified as

(11)
gK=g(K)=iK*mi+{i,j}K*mij,KC

Theorem 1. The fuzzy measure identified by steps 1 to 6 is a 2-order additive fuzzy measure.

To prove that the fuzzy measure identified by steps 1 to 6 is a 2-order additive fuzzy measure, it is only necessary to prove that the determined Möbius transformation coefficients satisfy the following constrained conditions [6]:

1) m (∅) =0; 2) mi ⩾ 0, ∀i ∈ C*; 3) ∑iC*mi + ∑{i,j}⊂C*mij = 1; 4) mi + ∑jK*∖{i}mij ⩾ 0, ∀K ⊂ C.

Proof.

  • 1) m (∅) =0, obviously holds.

  • 2) Because S(ξ˜ij)=S(ξ˜ji) , and i=1nwi=1 , P can be further written as

P=1+12i=1nj=1,ijnS(ξ˜ij)wiwj=1+12i=1n[wij=1,ijnS(ξ˜ij)wj]

Since S(ξ˜ij)[-1,1] , and wj ∈ [0, 1], we have

-wjS(ξ˜ij)wjwj,ij

Sum the two sides of the above inequality to j, we obtain

-(1-wi)j=1,ijnS(ξ˜ij)wj1-wi

Multiply both sides of the above inequality by wi, we can get

-(wi-wi2)wij=1,ijnS(ξ˜ij)wjwi-wi2

Sum the two sides of the above inequality to i, the following inequality can be given

-(1-i=1nwi2)i=1n[wij=1,ijnS(ξ˜ij)wj]1-i=1nwi2

Therefore, we have

12+12i=1nwi21+12i=1n[wij=1,ijnS(ξ˜ij)wj]32-12i=1nwi2

That is to say

0<12+12i=1nwi2P

Because P > 0, and wi ⩾ 0, we can get mi=wiP0 .

  • 3) iC*mi+{i,j}C*mij=iC*wiP+{i,j}C*S(ξ˜ij)wiwjP=1P[iC*wi+{i,j}C*S(ξ˜ij)wiwj]=1 , obviously holds.

  • 4) mi+jK*{i}mij=wiP+jK*{i}S(ξ˜ij)wiwjP=wiP[1+jK*{i}S(ξ˜ij)wj]

Since S(ξ˜ij)[-1,1] , and wj ∈ [0, 1], we have -(1-wi)jK*{i}S(ξ˜ij)wj1-wi . Thus, the following inequality can be given

wiP[1+jK*{i}S(ξ˜ij)wj]wiP[1-(1-wi)]wi2P0

Hence, we get mi + ∑jK*∖{i}mij ⩾ 0, ∀K ⊂ C. Q.E.D.

4A MADM process based on the proposed method

Using the Choquet fuzzy integral as nonlinear integration operator, this section presents a MADM process based on the proposed method.

Step 1: Calculate the group intuitionistic fuzzy assessment value x˜jik .

Let Ci ={ Ci1, Ci2, ⋯ , Cik, ⋯ , Cini } (k = 1, 2, ⋯ , ni) be a set of sub-attributes under the i-th attribute. According to the type of sub-attribute (such as positive, negative and neutral), the expert Dl (l = 1, 2, ⋯ , p) employs the IFN to assess the performance of the alternative Aj (j = 1, 2, ⋯ , m) under the sub-attribute Cik (i = 1, 2, ⋯ , n), the individual intuitionistic fuzzy assessment value x˜jikl=(μx˜jikl,νx˜jikl) of the k-th sub-attribute of the alternative Aj under the i-th attribute is then given.

Given the weight vector WD of experts, using Equation (1) to aggregate the opinions of p experts, the group intuitionistic fuzzy assessment value x˜jik of the k-th sub-attribute of the alternative Aj under the i-th attribute is further calculated as

(12)
x˜jik=(μx˜jik,νx˜jik)=(1-l=1p(1-μx˜jikl)wl,l=1pνx˜jiklwl)

Step 2: Calculate the overall weight wi of the i-th attribute and the relative weight wik of the k-th sub-attribute under the i-th attribute.

By using Equation (2), the intuitionistic fuzzy entropy Eik of the k-th sub-attribute under the i-th attribute is calculated as

(13)
Eik=1mj=1mmin(μx˜jik,νx˜jik)+πx˜jikmax(μx˜jik,νx˜jik)+πx˜jik

According to the definition of information entropy, if the information entropy of an attribute is smaller, the information quantity provided by the attribute is larger, and the weight of the attribute is also greater [29]. Therefore, the entropy weight wik of the k-th sub-attribute under the i-th attribute can be expressed as [30]

(14)
wik=1-Eiki=1nni-i=1nk=1niEik

That is to say, the overall weight vector of sub-attributes is Wi = (wi1, wi2, ⋯ , wik, ⋯ , wini), i = 1, 2, ⋯ , n. The overall weight wi of the i-th attribute is then obtained as

(15)
wi=k=1niwik

It can be proved that i=1nwi=1 holds. Therefore, the overall weight vector of attributes is WC = (w1, w2, ⋯ , wi, ⋯ , wn). The relative weight wik of the k-th sub-attribute under the i-th attribute is then obtained as

(16)
wik=wikwi

Step 3: Calculate the intuitionistic fuzzy assessment score SL(x˜ji) of the alternative Aj under the i-th attribute.

According to x˜jik=(μx˜jik,νx˜jik) and wik , with Equation (1), the intuitionistic fuzzy assessment value x˜ji of the alternative Aj under the i-th attribute is calculated as

(17)
x˜ji=(μx˜ji,νx˜ji)=(1-k=1ni(1-μx˜jik)wik,k=1niνx˜jikwik)

Hence, with Equation (4), the intuitionistic fuzzy assessment score SL(x˜ji) of the alternative Aj under the i-th attribute is further calculated as

(18)
SL(x˜ji)=μx˜ji+μx˜ji(1-μx˜ji-νx˜ji),SL(x˜ji)[0,1]

Step 4: Identify the 2-order additive fuzzy measure gK.

Given the weight vector WD of experts, the intuitionistic fuzzy interaction degree S(ξ˜ij) between Ci and Cj (ij) can be determined by steps 1 to 4 in Section 3. Thus, according to wi and S(ξ˜ij) , the 2-order additive fuzzy measure gK can be further identified by step 5 and step 6 in Section 3.

Step 5: Calculate the Choquet fuzzy integral values and the ranking of alternatives.

Reordering the intuitionistic fuzzy assessment score SL(x˜ji) (i = 1, 2, ⋯ , n) of the alternative Aj (j = 1, 2, ⋯ , m) under the i-th attribute from small to large, the SL(x˜j(i)) can be then obtained. Substituting the SL(x˜j(i)) and the 2-order additive fuzzy measure gK into Equation (7), the Choquet fuzzy integral value Hj of the alternative Aj can be calculated. Simultaneously, the ranking of alternatives can be given, where the larger Hj is, the better the alternative Aj is.

5Application example

This section takes the credit evaluation of the big data listed companies in China as an application example to illustrate the feasibility and effectiveness of the proposed method (see Section 3).

5.1Credit evaluation Index system and sample data

Considering the characteristics of big data enterprises [31], and following the principles of selecting indicators, such as scientificalness, objectivity, systematization, functionality, dynamics, relative independence, feasibility (or operability), comparability and so on, a credit evaluation index system for big data enterprises was constructed, as shown in Table 1.

Table 1

Credit evaluation index system for big data enterprises

Primary IndicatorsSecondary Indicators
Debt Paying Ability C1Current Ratio C11; Quick Ratio C12; Asset-liability Ratio C13
Operational Capability C2Turnover Rate of Accounts Receivable C21; Turnover Rate of Total Assets C22; Inventory Turnover C23
Profitability C3Profit Margin of Main Business C31; Return on Equity C32; Return on Total Assets C33
Growth Capability C4Net Profit Growth Rate C41; Growth Rate of Main Business C42
Technological Innovation Capability C5Development Expenditure C51; Growth Rate of Intangible Assets C52; Number of Invention Patent Applications Announced C53
Industry Growth C6Network Attention of Industry C61; Industry Average Net Profit Growth Rate C62

We selected the big data listed companies in the Growth Enterprise Market (GEM) in China – Wangsu Science & Technology Co., Ltd. (300017), Beijing Lanxum Technology Co., Ltd. (300010) and Wuhan Tianyu Information Industry Co., Ltd. (300205) to form a set of alternatives, denoted by A = A1, A2, A3. Where the alternative A1 and A2 belong to the software service industry, and the alternative A3 belongs to the electronic components industry. The sample data were the section data of 2016, and the original data were shown in Table 2. Where the original data of the Number of Invention Patent Applications Announced were from Tian Yan Cha website, the original data of the Network Attention of Industry were from Baidu Index website, and the rest of the original data were from East Money website.

Table 2

Original data

Secondary Indicators (Unit, Types)A1A2A3
C11 (—, *)6.731.672.48
C12 (—, *)6.491.381.76
C13 (%, −)15.0828.5233.54
C21 (time, +)5.482.813.40
C22 (time, +)0.730.370.97
C23 (time, +)17.322.782.88
C31 (%, +)28.0715.940.28
C32 (%, +)17.085.211.79
C33 (%, +)20.525.830.27
C41 (%, +)50.52130.2512.38
C42 (%, +)51.6784.109.71
C51 (10 thousand yuan, +)560598680
C52 (%, +)60.5751.24-6.73
C53 (piece, +)83024
C61 (time, +)159159457
C62 (%, +)30.7730.7733.42

Note: ding172“+” indicates positive index, “−” indicates negative index, “*” indicates neutral index; ding173The Network Attention of Industry is the overall daily average of Baidu Search Index with the industry name as the key word.

5.2Process and results of credit evaluation

Step 1: Calculate the group intuitionistic fuzzy assessment value x˜jik of the k-th sub-attribute of the alternative Aj under the i-th attribute.

According to the type of sub-attribute (such as positive, negative and neutral), three experts employed the IFN to assess the performance of the alternative Aj under the sub-attribute Cik respectively, the individual intuitionistic fuzzy assessment value x˜jikl of the k-th sub-attribute of the alternative Aj under the i-th attribute was then given.

Adopting the cycle mutual evaluation method [32], the weight vector of experts was calculated as WD = (0.3976, 0.3012, 0.3012) (see Appendix A for full calculation principle and process).

Using Equation (12) to aggregate the opinions of three experts, the group intuitionistic fuzzy assessment value x˜jik of the k-th sub-attribute of the alternative Aj under the i-th attribute was further obtained, as shown in Table 3.

Table 3

Group intuitionistic fuzzy assessment values of sub-attributes

Sub-attributesA1A2A3
C11(0.399, 0.429)(0.868, 0.108)(0.800, 0.150)
C12(0.248, 0.664)(0.868, 0.108)(0.704, 0.183)
C13(0.551, 0.266)(0.720, 0.178)(0.763, 0.164)
C21(0.868, 0.108)(0.610, 0.226)(0.704, 0.183)
C22(0.704, 0.183)(0.308, 0.565)(0.763, 0.164)
C23(0.950, 0.050)(0.566, 0.255)(0.566, 0.255)
C31(0.704, 0.183)(0.500, 0.300)(0.050, 0.950)
C32(0.924, 0.070)(0.459, 0.350)(0.248, 0.664)
C33(0.950, 0.050)(0.566, 0.255)(0.050, 0.950)
C41(0.763, 0.164)(0.950, 0.050)(0.551, 0.266)
C42(0.763, 0.164)(0.950, 0.050)(0.459, 0.350)
C51(0.596, 0.235)(0.868, 0.108)(0.050, 0.950)
C52(0.924, 0.070)(0.868, 0.108)(0.050, 0.950)
C53(0.924, 0.070)(0.050, 0.950)(0.566, 0.255)
C61(0.610, 0.226)(0.610, 0.226)(0.924, 0.070)
C62(0.763, 0.164)(0.763, 0.164)(0.885, 0.097)

Step 2: Determine the weights of attributes.

Based on Table 3, with Equations (13) and (14), the intuitionistic fuzzy entropy Eik and entropy weight wik of the k-th sub-attribute under the i-th attribute were calculated. In addition, with Equation (15), the overall weight vector WC of attributes was also obtained. With Equation (16), the relative weight wik of the k-th sub-attribute under the i-th attribute was then obtained. The above calculation results were shown in Table 4.

Table 4

Weights of attributes

AttributesOverall WeightsSub-attributesIntuitionistic Fuzzy EntropyEntropy WeightsRelative Weights
C10.170C110.4450.0520.304
C120.3190.0630.373
C130.4120.0550.323
C20.171C210.3380.0620.361
C220.4250.0540.314
C230.4060.0550.325
C30.181C310.3760.0580.321
C320.4530.0510.282
C330.2290.0720.397
C40.121C410.3160.0640.528
C420.3890.0570.472
C50.226C510.2430.0710.312
C520.0940.0840.374
C530.2390.0710.314
C60.131C610.3630.0590.453
C620.2310.0720.547

Step 3: Calculate the intuitionistic fuzzy assessment value x˜ji of the alternative Aj under the i-th attribute.

Based on Tables 3 and 4, using Equation (17), the intuitionistic fuzzy assessment value x˜ji of the alternative Aj under the i-th attribute was calculated as x˜11=(0.405,0.433) , x˜12=(0.876,0.099) , x˜13=(0.900,0.083) , x˜14=(0.763,0.164) , x˜15=(0.872,0.102) , x˜16=(0.703,0.189) ; x˜21=(0.832,0.127) , x˜22=(0.517,0.313) , x˜23=(0.517,0.294) , x˜24=(0.950,0.050) , x˜25=(0.755,0.214) , x˜26=(0.703,0.189) ; x˜31=(0.756,0.166) , x˜32=(0.688,0.197) , x˜33=(0.111,0.859) , x˜34=(0.510,0.302) , x˜35=(0.257,0.629) , x˜36=(0.905,0.083) .

Step 4: Calculate the intuitionistic fuzzy assessment score SL(x˜ji) of the alternative Aj under the i-th attribute.

According to the intuitionistic fuzzy assessment value x˜ji of the alternative Aj under the i-th attribute, with Equation (18), the intuitionistic fuzzy assessment score SL(x˜ji) of the alternative Aj under the i-th attribute was calculated as SL(x˜11)=0.4709 , SL(x˜12)=0.8976 , SL(x˜13)=0.9151 , SL(x˜14)=0.8192 , SL(x˜15)=0.8949 , SL(x˜16)=0.7789 ; SL(x˜21)=0.8662 , SL(x˜22)=0.6048 , SL(x˜23)=0.6148 , SL(x˜24)=0.9500 , SL(x˜25)=0.7786 , SL(x˜26)=0.7789 ; SL(x˜31)=0.8147 , SL(x˜32)=0.7671 , SL(x˜33)=0.1140 , SL(x˜34)=0.6054 , SL(x˜35)=0.2868 , SL(x˜36)=0.9156 .

Step 5: Determine the intuitionistic fuzzy interaction degree S(ξ˜ij) between Ci and Cj (ij).

According to step 2 in Section 3, three experts employed the IFN α = (μα, να) to assess the interactivity between Ci and Cj (ij) ( C62 pairs in total) respectively. Thus, the individual assessment result ξ˜ijl of the interactivity between Ci and Cj (ij) was then given, as shown in Table 5.

Table 5

Intuitionistic fuzzy interaction degrees between attributes

Attributes Ci and Cj (ij)Individual Assessment Results of InteractivityGroup Assessment Results of InteractivityIntuitionistic Fuzzy Interaction Degrees
Expert 1Expert 2Expert 3
{C1, C2}(0.50, 0.30)(0.65, 0.20)(0.50, 0.30)(0.551, 0.266)0.285
{C1, C3}(0.65, 0.20)(0.80, 0.15)(0.80, 0.15)(0.750, 0.168)0.582
{C1, C4}(0.50, 0.30)(0.65, 0.20)(0.50, 0.30)(0.551, 0.266)0.285
{C1, C5}(0.50, 0.30)(0.65, 0.20)(0.50, 0.30)(0.551, 0.266)0.285
{C1, C6}(0.50, 0.30)(0.50, 0.50)(0.50, 0.30)(0.500, 0.350)0.150
{C2, C3}(0.65, 0.20)(0.80, 0.15)(0.80, 0.15)(0.750, 0.168)0.582
{C2, C4}(0.50, 0.30)(0.50, 0.30)(0.65, 0.20)(0.551, 0.266)0.285
{C2, C5}(0.50, 0.30)(0.50, 0.30)(0.65, 0.20)(0.551, 0.266)0.285
{C2, C6}(0.50, 0.30)(0.50, 0.50)(0.50, 0.30)(0.500, 0.350)0.150
{C3, C4}(0.65, 0.20)(0.80, 0.15)(0.65, 0.20)(0.704, 0.183)0.521
{C3, C5}(0.50, 0.30)(0.65, 0.20)(0.65, 0.20)(0.596, 0.235)0.361
{C3, C6}(0.50, 0.30)(0.65, 0.20)(0.65, 0.20)(0.596, 0.235)0.361
{C4, C5}(0.50, 0.30)(0.65, 0.20)(0.65, 0.20)(0.596, 0.235)0.361
{C4, C6}(0.20, 0.65)(0.15, 0.80)(0.15, 0.80)(0.170, 0.736)–0.566
{C5, C6}(0.50, 0.30)(0.50, 0.50)(0.50, 0.30)(0.500, 0.350)0.150

According to step 3 in Section 3, given the weight vector WD = (0.3976, 0.3012, 0.3012) of the experts, using Equation (8) to aggregate the opinions of three experts, the group assessment result ξ˜ij of the interactivity between Ci and Cj (ij) was obtained (see Table 5). Taking ξ˜12 as an example, we had

ξ˜12=(1-(1-0.50)0.3976×(1-0.65)0.3012×(1-0.50)0.3012,0.300.3976×0.200.3012×0.300.3012)=(0.551,0.266).

With Equation (9), the intuitionistic fuzzy interaction degree S(ξ˜ij) between Ci and Cj (ij) was calculated, as shown in Table 5. Taking S(ξ˜12) as an example, we had S(ξ˜12)=0.551-0.266=0.285 .

Step 6: Calculate the Möbius transformation coefficients mi and mij of attributes.

Based on Tables 4 and 5, using Equation (10), the Möbius transformation coefficients mi and mij of attributes were calculated, as shown in Table 6, where P = 1.125. Taking m12 as an example, we had

Table 6

Calculation results of Möbius transformation coefficients

M öbius Transformation CoefficientsCoefficient ValuesMöbius Transformation CoefficientsCoefficient ValuesMöbius Transformation CoefficientsCoefficient Values
m10.151m130.016m260.003
m20.152m140.005m340.010
m30.161m150.010m350.013
m40.108m160.003m360.008
m50.201m230.016m450.009
m60.116m240.005m46–0.008
m120.007m250.010m560.004

m12=(0.285×0.170×0.171)/1.125=0.007.

Step 7: Identify the 2-order additive fuzzy measure gK.

Based on Table 6, using Equation (11), the 2-order additive fuzzy measure gK was calculated, as shown in Table 7. Taking g{1,2} as an example, we had

Table 7

Calculation results of 2-order additive fuzzy measure

KgKKgKKgKKgK
{ø}0.000{3, 4}0.279{2, 3, 4}0.452{1, 3, 4, 5}0.683
{1}0.151{3, 5}0.375{2, 3, 5}0.553{1, 3, 4, 6}0.570
{2}0.152{3, 6}0.285{2, 3, 6}0.456{1, 3, 5, 6}0.683
{3}0.161{4, 5}0.317{2, 4, 5}0.484{1, 4, 5, 6}0.599
{4}0.108{4, 6}0.216{2, 4, 6}0.376{2, 3, 4, 5}0.684
{5}0.201{5, 6}0.321{2, 5, 6}0.486{2, 3, 4, 6}0.571
{6}0.116{1, 2, 3}0.503{3, 4, 5}0.501{2, 3, 5, 6}0.684
{1, 2}0.310{1, 2, 4}0.428{3, 4, 6}0.395{2, 4, 5, 6}0.600
{1, 3}0.328{1, 2, 5}0.531{3, 5, 6}0.503{3, 4, 5, 6}0.621
{1, 4}0.264{1, 2, 6}0.433{4, 5, 6}0.430{1, 2, 3, 4, 5}0.874
{1, 5}0.362{1, 3, 4}0.451{1, 2, 3, 4}0.631{1, 2, 3, 4, 6}0.753
{1, 6}0.271{1, 3, 5}0.552{1, 2, 3, 5}0.737{1, 2, 3, 5, 6}0.871
{2, 3}0.329{1, 3, 6}0.455{1, 2, 3, 6}0.633{1, 2, 4, 5, 6}0.776
{2, 4}0.265{1, 4, 5}0.483{1, 2, 4, 5}0.658{1, 3, 4, 5, 6}0.806
{2, 5}0.363{1, 4, 6}0.375{1, 2, 4, 6}0.543{2, 3, 4, 5, 6}0.807
{2, 6}0.271{1, 5, 6}0.485{1, 2, 5, 6}0.657{1, 2, 3, 4, 5, 6}1.000

g{1,2}=0.151+0.152+0.007=0.310.

Step 8: Calculate the Choquet fuzzy integral values and the ranking of alternatives.

Take the alternative A1 as an example. Reordering the intuitionistic fuzzy assessment score SL(x˜1i) (i = 1, 2, ⋯ , 6) of alternative A1 under the i-th attribute from small to large, we can get

SL(x˜11)<SL(x˜16)<SL(x˜14)<SL(x˜15)<SL(x˜12)<SL(x˜13).

It can be denoted by

SL(x˜1(1))<SL(x˜1(2))<SL(x˜1(3))<SL(x˜1(4))<SL(x˜1(5))<SL(x˜1(6)).

Substituting the SL(x˜1(i)) and the 2-order additive fuzzy measure gK into Equation (7), the Choquet fuzzy integral value of the alternative A1 was calculated as

H1=(0.4709-0.0000)×1.000+(0.7789-0.4709)×0.807+(0.8192-0.7789)×0.684+(0.8949-0.8192)×0.553+(0.8976-0.8949)×0.329+(0.9151-0.8976)×0.161=0.7927.

Similarly, we can also obtain H2 = 0.7652 and H3 = 0.5157. Since H1 > H2 > H3, then the ranking of alternatives was A1 ≻ A2 ≻ A3.

That is to say, the credit status of the alternative A1 was relatively good, and the credit status of the alternative A3 was relatively poor.

5.3Comparative analysis

For comparison, we also replaced the intuitionistic fuzzy interaction degrees between attributes (see Table 5) in step 5 of Section 5.2 with the scoring interaction degrees between attributes (see Table 8) determined by the interaction degrees scoring method [13]. Ceteris paribus, the Choquet fuzzy integral values of alternatives were calculated as H1=0.7946 , H2=0.7686 , H3=0.5244 . Since H1>H2>H3 , then the ranking of alternatives was A1 ≻ A2 ≻ A3.

Table 8

Scoring interaction degrees between attributes

Attributes Ci and Cj (ij)Individual Scoring Interaction DegreesGroup Scoring Interaction Degrees
Expert 1Expert 2Expert 3
{C1, C2}0.200.300.200.23
{C1, C3}0.300.500.500.42
{C1, C4}0.200.300.200.23
{C1, C5}0.200.300.200.23
{C1, C6}0.200.000.200.14
{C2, C3}0.300.500.500.42
{C2, C4}0.200.200.300.23
{C2, C5}0.200.200.300.23
{C2, C6}0.200.000.200.14
{C3, C4}0.300.500.300.36
{C3, C5}0.200.300.300.26
{C3, C6}0.200.300.300.26
{C4, C5}0.200.300.300.26
{C4, C6}–0.30–0.50–0.50–0.42
{C5, C6}0.200.000.200.14

Note: The weighted arithmetic mean method was employed to aggregate the opinions of three experts.

We further investigated the discrimination of the MADM method based on the intuitionistic fuzzy interaction degrees (hereinafter referred to as Method 1) for alternatives, as well as the MADM method based on the scoring interaction degrees (hereinafter referred to as Method 2). Adopting the algorithm of discrimination given in [33], the discrimination of Method 1 for alternatives was calculated as

ρ=0.7927-0.76520.7927+0.7927-0.51570.7927+0.7652-0.51570.7652=0.7102.

Simultaneously, the discrimination of Method 2 for alternatives was calculated as

ρ=0.7946-0.76860.7946+0.7946-0.52440.7946+0.7686-0.52440.7686=0.6903.

Although the ranking result of Method 2 is consistent with that of Method 1, the discrimination of Method 1 for alternatives is higher than that of Method 2. Which means that the decision-making result of Method 1 is better than that of Method 2.

6Discussion

From the results and analysis of the previous section, we observed that Method 1 was able to obtain the higher discrimination value than Method 2 (Method 1 was 0.7102, and Method 2 was 0.6903), and the slightly lower Choquet fuzzy integral mean value than Method 2 (Method 1 was 0.6912, and Method 2 was 0.6959).

Compared with Method 2, Method 1 can obtain the higher discrimination value. Since the variance of interaction degrees determined by Method 1 (its value was equal to 0.0684) was higher than that of Method 2 (its value was equal to 0.0355). According to Equation (6), the variances of mi and mij increase with the increase of the variance of interaction degrees. From Equations (5) and (7), we can further see that the variance of g (K) and the variance of H also increase correspondingly. Thus, according to the algorithm of discrimination [33], the discrimination value becomes larger.

Compared with Method 2, Method 1 can obtain the slightly lower Choquet fuzzy integral mean value. Since the average value of interaction degrees determined by Method 1 (its value was equal to 0.2718) was higher than that of Method 2 (its value was equal to 0.2087). According to Equation (6), when the average value of interaction degrees increases, the mean value of P increases, meanwhile, the mean value of mi decreases and the mean value of mij increases. However, because |mi| ⪢ |mij|, from Equation (5), we can further see that the mean value of g (K) also decreases correspondingly. Thus, according to Equation (7), the mean value of H becomes smaller.

The variance or average value of interaction degrees determined by Method 1 was higher than that of Method 2, which was closely related to that the experts used the IFN to assess the interactivity between attributes, so it can fully express their opinions.

7Conclusion

In this paper, the proposed method defines the interactivity between attributes by using the supermodular game theory, so that the interaction between attributes is easier to explain and understand, which lays a solid foundation for experts to qualitatively assess the interactivity between attributes. The proposed method allows the experts to assess the interactivity between attributes by using the IFN, which better preserves the assessment information of experts and embodies the fuzziness and hesitation of interactivity assessment. Furthermore, the proposed method uses the IFWA to aggregate the opinions of all experts, which considers the weights of experts, thereby ensuring the rationality of decision-making. In addition, the proposed method uses the score function of IFN to define and calculate the intuitionistic fuzzy interaction degree between attributes, so the transformation from qualitative description to quantitative characterization is finally realized. Therefore, using the IFSs, this work successfully solves the problem that it is difficult to quantitatively assess the interactivity between attributes in the identification process of 2-order additive fuzzy measure.

This work proposed a 2-order additive fuzzy measure identification method based on IFSs. Obviously, on the one hand, compared with the objective method describing and dealing with the interactivity between attributes (as in [7, 16] and [17]), the proposed method has subjectivity. On the other hand, compared with the subjective method describing and dealing with the interactivity between attributes (as in [9–14] and [15]), the proposed method has fuzziness and hesitation.

Future application example analysis will consider increasing the number of samples and experts to improve the persuasiveness of application example analysis results and the stability of interactivity assessment. Furthermore, because the sample time may affect the analysis results of application example, it is necessary to update the sample data from 2016 to 2020 in the future. In addition, Method 1 should be compared with other methods except Method 2, such as the diamond pairwise comparisons method [9, 10, 14], proportional scaling method [11], multicriteria correlation preference information method [12], etc.

According to the supermodular game theory, Definition 9 gives the definition of interactivity between attributes. However, this definition is still relatively general. The detailed theoretical analysis of the connotation of interactivity between attributes needs to be completed in the future. The existing 2-order additive fuzzy measure identification methods use the subjective methods and the objective methods to describe and deal with the interactivity between attributes. Thus, combining the results of subjective methods and objective methods is one of the main research directions in the future. In addition, constructing a Group Decision Support System will reduce the complexity of Method 1 in practical application.

Acknowledgments

To the Regional Project of National Natural Science Foundation of China (71861003) and the Innovative Exploration and New Academic Seedlings Project of Guizhou University of Finance and Economics (Guizhou-Science Cooperation Platform Talents [2018] 5774-016) for their support.

References

[1] 

Marichal J.L. , An axiomatic approach of the discrete Choquet integral as a tool to aggregate interacting criteria, IEEE Transactions on Fuzzy Systems 8: (6) ((2000) ), 800–807.

[2] 

Sugeno M. , Theory of Fuzzy Integrals and Its Application, Ph.D. Dissertation, Tokyo Institute of Technology ((1974) ).

[3] 

Murofushi T. and Sugeno M. , A theory of fuzzy measures: Representations, the Choquet integral, and null sets, Journal of Mathematical Analysis and Applications 159: (2) ((1991) ), 532–549.

[4] 

Wang Z.Y. and Klir G.J. , Fuzzy Measure Theory, Plenum Publishing Corporation, New York, ((1992) ).

[5] 

Ishii K. and Sugeno M. , A Model of Human Evaluation Process Using Fuzzy Measure, International Journal of Man-Machine Studies 22: (1) ((1985) ), 19–38.

[6] 

Grabisch M. , K-order Additive Discrete Fuzzy Measures and Their Representation, Fuzzy Sets and Systems 92: (2) ((1997) ), 167–189.

[7] 

Wu J.Z. and Zhang Q. , 2-order additive fuzzy measures identification method based on maximum entropy principle, Systems Engineering and Electronics 32: (11) ((2010) ), 2346–2351.

[8] 

Mayag B. and Bouyssou D. , Necessary and possible interaction between criteria in a 2-additive Choquet integral model, European Journal of Operational Research 283: (1) ((2020) ), 308–320.

[9] 

Takahagi E. , A fuzzy Measure Identification Method by Diamond Pairwise Comparisons and ø(s) Transformation, Fuzzy Optimization and Decision Making 7: (3) ((2008) ), 219–232.

[10] 

Wu J.Z. and Zhang Q. , Multicriteria decision making method based on 2-order additive fuzzy measures, Systems Engineering-Theory & Practice 30: (7) ((2010) ), 1229–1237.

[11] 

Zhang Z. , Fang K. and Yang M. , Method for simulation credibility evaluation based on 2-additive fuzzy measures, Control and Decision 28: (1) ((2013) ), 147–151.

[12] 

Wu J.Z. , Yang S.L. , Zhang Q. and Ding S. , 2-additive Capacity Identification Methods from Multicriteria Correlation Preference Information, IEEE Transactions on Fuzzy Systems 23: (6) ((2015) ), 2094–2106.

[13] 

Chang Z.P. and Cheng L.S. , Grey Fuzzy Integral Correlation Degree Decision Model, Chinese Journal of Management Science 23: (11) ((2015) ), 105–111.

[14] 

Lin S.W. and Jerusalem M.A. , Hybrid method for assessing marketable fashion designs: diamond pairwise comparison combined with the Choquet integral, Textile Research Journal 88: (13) ((2018) ), 1531–1542.

[15] 

Kadaifci C. , Asan U. and Bozdag E. , A new 2-additive Choquet integral based approach to qualitative cross-impact analysis considering interaction effects, Technological Forecasting and Social Change 158: ((2020) ), 120131.

[16] 

Chang Z.P. , Cheng L.S. and Liu J.S. , Multi-Attribute Decision Making Method Based on Mahalanobis-Taguchi System and 2-additive Choquet integral, Journal of Industrial Engineering and Engineering Management 30: (1) ((2016) ), 133–139.

[17] 

Chang Z.P. , Cheng L.S. , Lu S.Z. and Wang X.Z. , Exit decision model of public rental housing based on 2-order additive fuzzy measures, Journal of Nanjing University of Science and Technology (Nature Science) 41: (1) ((2017) ), 132–138.

[18] 

Atanassov K.T. , Intuitionistic Fuzzy Sets, Fuzzy Sets and Systems 20: (1) ((1986) ), 87–96.

[19] 

Xu Z.S. , Intuitionistic fuzzy information aggregation theory and its application, Science Press, Beijing ((2008) ).

[20] 

Tan C.Q. , TOPSIS Multiple Attribute Decision Making Based on Interval-valued Intuitionistic Fuzzy Sets, Fuzzy Systems and Mathematics 24: (1) ((2010) ), 92–97.

[21] 

Wang K.R. and Yang H. , Multiple Attribute Decision Making Method Based on Intuitionistic Fuzzy Sets, Fuzzy Systems and Mathematics 24: (3) ((2010) ), 114–118.

[22] 

Verma R. and Merigó J.M. , Multiple attribute group decision making based on 2-dimension linguistic intuitionistic fuzzy aggregation operators, Soft Computing 24: ((2020) ), 17377–17400.

[23] 

Murofushi T. and Sugeno M. , An Interpretation of Fuzzy Measures and the Choquet Integral as an Integral with Respect to a Fuzzy Measure, Fuzzy Sets and Systems 29: (2) ((1989) ), 201–227.

[24] 

Xu Z.S. , Intuitionistic Fuzzy Aggregation Operators, IEEE Transactions on Fuzzy Systems 15: (6) ((2007) ), 1179–1187.

[25] 

Wang Y. and Lei Y.J. , A technique for constructing intuitionistic fuzzy entropy, Control and Decision 22: (12) ((2007) ), 1390–1394.

[26] 

Chen S.M. and Tan J.M. , Handling Multi-criteria Fuzzy Decision-making Problems Based on Vague Set Theory, Fuzzy Sets and Systems 67: (2) ((1994) ), 163–172.

[27] 

Liu H.W. , Vague Set Methods of Multicriteria Fuzzy Decision Making, Systems Engineering-Theory & Practice 24: (5) ((2004) ), 103–109.

[28] 

Topkis D.M. , Minimizing a submodular function on a lattice, Operations Research 26: (2) ((1978) ), 305–321.

[29] 

Qiu W.H. , Manage decision and entropy, Mechanical Industry Press, Beijing, ((2002) ).

[30] 

Meng J.N. , Fang N. , Liu B.S. and Xue B. , Sustainability Evaluation Model for Urban Infrastructure Projects Based on Intuitionistic Fuzzy Sets, Fuzzy Systems and Mathematics 29: (1) ((2015) ), 167–174.

[31] 

Zhang M. and Wang X.P. , Evaluation on Risk of Intellectual Property Pledge Financing for Big Data Enterprises, Science and Technology Management Research 39: (21) ((2019) ), 61–66.

[32] 

Li Y.H. , Evaluation and selection of strategic emerging industries based on Fuzzy AHP considering the weight of experts: An empirical analysis of Tangshan, China Collective Economy (6) ((2017) ), 60–62.

[33] 

Ke H.F. , Chen Y.G. and Xia B. , An Algorithm of Multiple Criteria Decision-Making Based on Similarity to Ideal Grey Relational Projection, Acta Electronica Sinica 35: (9) ((2007) ), 1757–1761.

Appendices

Appendix A

In step 1 of Section 5.2, the weight vector of experts was calculated by using the cycle mutual evaluation method [32]. Its calculation principle is as follows:

Suppose there are p experts, bij is the mutual evaluation weight given by expert i to expert j, then the mutual evaluation weight matrix is B = (bijp×p, where 0 ⩽ bij ⩽ 1, and j=1pbij=1 . Then the weight of expert j given by the group is qj=p-1i=1pbij . Therefore, the final weights of experts are determined through the cycle evaluation. Let t be the number of mutual evaluations, when t = 1, the weight vector of experts is q1 = (q1, q2, ⋯ , qp); when t > 1, the weight vector of experts is qt = qt-1B, until qt converges to the stable value.

In this paper, the mutual evaluation weight matrix was given by three experts, it was

B=[0.500.250.250.380.310.310.280.360.36]

When t = 1, the weight vector of experts was q1 = (0.3866, 0.3067, 0.3067); when t = 5, the weight vector of experts converged to the stable value q5 = (0.3976, 0.3012, 0.3012). Therefore, the final weight vector of experts was WD = (0.3976, 0.3012, 0.3012).