Next Article in Journal
A Top-Down Interactive Visual Analysis Approach for Physical Simulation Ensembles at Different Aggregation Levels
Next Article in Special Issue
Dombi Aggregation Operators of Linguistic Cubic Variables for Multiple Attribute Decision Making
Previous Article in Journal
AAC Double Compression Audio Detection Algorithm Based on the Difference of Scale Factor
Previous Article in Special Issue
Pythagorean Fuzzy Interaction Muirhead Means with Their Application to Multi-Attribute Group Decision-Making
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Green Supplier Assessment Method for Manufacturing Enterprises Based on Rough ANP and Evidence Theory

1
College of Mechatronic Engineering, North Minzu University, Yinchuan 750021, China
2
The 713th Research Institute of China Shipbuilding Industry Corporation, Zhengzhou 450052, China
*
Author to whom correspondence should be addressed.
Information 2018, 9(7), 162; https://doi.org/10.3390/info9070162
Submission received: 10 May 2018 / Revised: 25 June 2018 / Accepted: 28 June 2018 / Published: 2 July 2018

Abstract

:
Within the context of increasingly serious global environmental problems, green supplier assessment has become one of the key links in modern green supply chain management. In the actual work of green supplier assessment, the information of potential suppliers is often ambiguous or even absent, and there are interrelationships and feedback-like effects among assessment indexes. Additionally, the thinking of experts in index importance judgment is always ambiguous and subjective. To handle the uncertainty and incompleteness in green supplier assessment, we propose a green supplier assessment method based on rough ANP and evidence theory. The uncertain index value is processed by membership degree. Trapezoidal fuzzy number is adopted to express experts’ judgment on the relative importance of the indexes, and rough boundary interval is used to integrate the judgment opinions of multiple experts. The ANP structure is built to deal with the interrelationship and feedback-like effects among indexes. Then, the index weight is calculated by ANP method. Finally, the green suppliers are assessed by a trust interval, based on evidence theory. The feasibility and effectiveness of the proposed method is verified by an application of a bearing cage supplier assessment.

1. Introduction

With the increasing awareness of global environmental protection and the increasing number of related environmental regulations, manufacturing enterprises are facing more stringent environmental requirements. Nowadays, green supply chains have become inevitable choices for manufacturing enterprises who wish to deal with environmental problems. Green supply chain management includes many links, such as green supplier assessment, green product design, green production, green marketing and waste recovery [1,2,3]. Green supplier assessment is in the upstream of the whole supply chain, and its effect on environmental protection and cost saving can be transmitted to every part of the downstream through the supply chain.
In the process of green supply chain management, various factors make the relationship between suppliers and manufacturing enterprises complicated and vague. However, competitors constantly adjust their strategies, and the supply chain must constantly improve to adapt to the complex environment which changes rapidly. In this context, green supplier assessment plays a very important role in reducing costs, and improving product quality and market competitiveness. Through effective assessment and supervisions of suppliers, problems can be found and solved in time, and the green and healthy development of the entire supply chain can be promoted [4,5,6,7,8,9]. It can be seen that green supplier assessment plays a decisive role in green supply chain management, which directly determines the competitiveness of the entire supply chain. The rise of the internet has provided convenience for manufacturing enterprises to assess green suppliers, but enterprises cannot quickly choose suppliers which meet their needs in the face of so many uneven suppliers. Considering finances and the effective utilization of resources, how to assess green supplier quickly and effectively becomes the key problem in modern green supply chain.
As seen in the present literature, there are many significant works on green supplier assessment. On the whole, the existing research mainly includes the following aspects.
(1)
Supplier assessment models based on Analytic Hierarchy Process (AHP) [10] or Analytic Network Process (ANP) [11]. Noci [12] used AHP to evaluate supplier’s environmental efficiency. Lee et al. [13] used the Delphi technique to distinguish the evaluation criterion difference between the traditional supplier and the green supplier, and then used Fuzzy Analytic Hierarchy Process (FAHP) to solve the green supplier selection process. Hsu and Hu [14] contained the interdependence between components of decision structure and used ANP for green supplier selection which reflected a more realistic result.
(2)
Supplier assessment model based on mathematical Programming. Yeh and Chuang [15] put forward a mathematical programming model of green partner selection, which includes four goals: cost, time, product quality, and green score. They adopted two multi-objective genetic algorithms (MOGA) to find a set of Pareto optimal solutions, and used weighted summation to generate more solutions. Yousefi et al. [16] used Dynamic Data Envelopment Analysis (D-DEA) and scenario-based robust model for supplier selection. In this supplier selection model, the shortcomings of the DEA model (the benchmarks were determined based on previous performance) were overcome, and the disadvantages of the D-DEA model (the decision unit couldn’t get a unified efficiency score) were avoided.
(3)
Supplier assessment model based on Technique for Order Preference by Similarity to an Ideal Solution (TOPSIS). Awasthi et al. [17] proposed a three-step method for green supplier evaluation: identification standard, expert score, evaluate expert score by fuzzy TOPSIS. The fuzzy TOPSIS method integrated profit and cost standard, and this method is suitable for the situation of lack of partial quantitative information. Kannan and Jabbour [18] used fuzzy TOPSIS to solve the problem of green supplier selection, and applies three types of fuzzy TOPSIS method to sort green suppliers.
(4)
Other hybrid models for supplier assessment. Gandhi et al. [5] proposed a combined approach using AHP and Decision-Making and Trial Evaluation Laboratory (DEMATEL) for evaluating success factors in implementation of green supply chain management and gave a case study in Indian manufacturing industries. Chatterjee et al. [19] combined DEMATEL and ANP in a rough context, and then proposed a rough DEMATEL-ANP (R’AMATEL) method to evaluate the performance of suppliers for green supply chain implementation in electronics industry. Wu et al. [20] used the Continuous Ordered Weighted Averaging (COWA) operator to transform the trapezoid fuzzy number into the exact real number to select the green supplier, and make a sensitivity analysis according to the degree of risk of decision-maker to rank the suppliers. Luo and Peng [21] proposed a multi-level supplier evaluation and selection model. In this model, AHP is used to determine the weights, and then TOPSIS is used for supplier evaluation. Kuo and Lin [22] integrated ANP and DEA and proposed a green supplier evaluation method. The interdependence between standards were considered by ANP, which allowed users to choose their own weight preferences to limit weights, expanded DEA method, and allowed more flexible number of decision units. Shi et al. [23] used the improved attribute reduction algorithm based on rough set to reduce the index of the green supplier evaluation index system, and then evaluated the data by RBF neural network training. Akman [24] identified the suppliers that should be included in the green supplier development plan through the C mean clustering algorithm and the VIKOR method. This method can be used to solve the problem of supplier classification and evaluation.
The study of green supplier assessment, which is of great theoretical and practical significance, has been a hot topic all along. However, the existing research has obvious shortcomings and the research gaps are mainly as follows:
(1)
The information of the suppliers to be assessed is not clear enough in the actual work of green supplier assessment. There is often no information sharing between manufacturing enterprise and suppliers, and the information between them is often ambiguous or even absent. The deterministic assessing model can no longer meet the needs of the increasingly complex decision-making environment.
(2)
Green supplier assessment is a complex decision problem and its indexes are interrelated. When calculating the index weight, the core idea of traditional AHP is to divide the index system into isolated and hierarchical levels. Only the upper level elements’ dominating effect on the lower level elements is considered and the elements in the same level are deemed to be independent of each other. However, the relationships among the indexes are often interdependent and sometimes provide feedback-like effects in green supplier assessment. Therefore, traditional AHP cannot solve the complex relationships among indexes to obtain the weight in green supplier assessment.
(3)
The accurate number is used to describe the relative importance of the indexes in the expression of experts’ judgment in most of the existing research, which cannot reflect the ambiguities and subjectivity of the actual thinking. It is more reasonable to use fuzzy numbers [25] to express the experts’ judgment. After introducing fuzzy numbers to express experts’ judgment, analyzing and processing the imprecise and inconsistent information becomes a difficult problem.
To fill the research gaps in green supplier assessment, a green supplier assessment method for manufacturing enterprises based on rough ANP and evidence theory is proposed. We process the uncertain index value by membership degree, adopt trapezoidal fuzzy number to express experts’ judgment on the relative importance of the indexes, use rough boundary interval to integrate the judgment opinions of multiple experts, set up the ANP structure to deal with the interrelationship and the feedback-like effects among indexes and then calculate the index weight by ANP method, and finally solve the incomplete information problem by evidence theory and assess the green suppliers by trust interval.
The rest of this paper is organized as follows. Section 2 establishes the index system of green supplier assessment; Section 3 uses membership degree method to process the uncertain index value of suppliers to be assessed; Section 4 adopts rough ANP to calculate the index weight; Section 5 gives the green supplier assessment procedure based on evidence theory; Section 6 provides an application case of bearing cage supplier assessment and discusses the feasibility and effectiveness of the proposed method for green supplier assessment. We conclude this paper in Section 7.

2. Index System

The first and very important segment of green supplier assessment is the establishment of a complete and overall index system. The attribute of the supplier’s product is the main representative of its ability, and the comprehensive ability of supplier can provide a strong support to its product. The comprehensive ability of supplier mainly includes internal competitiveness, external competitiveness and cooperation ability. Internal competitiveness of a supplier can be subdivided into its innovation capacity, manufacturing capacity and agility capacity. Furthermore, a supplier is not isolated and is inevitably restricted by its external competitiveness. External competitiveness of a supplier mainly includes its economic environment, geographical environment, social environment and legal environment. Additionally, cooperation ability of a supplier is affected by its technical compatibility degree, cultural compatibility degree, information platform compatibility degree, and reputation.
As shown in Figure 1, we establish the index system of green supplier assessment.
The green supplier assessment objective (AO) includes four first-level indexes: product attribute (C1), internal competitiveness (C2), external competitiveness (C3) and cooperation ability (C4).
  • C1 is decomposed into four second-level indexes: cost (C1,1), quality (C1,2), service (C1,3) and flexibility (C1,4). C1,1 and C1,2 belong to quantitative type, and C1,3 and C1,4 belong to qualitative type.
  • C2 is decomposed into three second-level indexes: innovation capacity (C2,1), manufacturing capacity (C2,2) and agility capacity (C2,3). C2,1, C2,2 and C2,3 all belong to qualitative type.
  • C3 is decomposed into four second-level indexes: economic environment (C3,1), geographical environment (C3,2), social environment (C3,3) and legal environment (C3,4). C3,1, C3,2, C3,3 and C3,4 all belong to qualitative type.
  • C4 is decomposed into four second-level indexes: technical compatibility degree (C4,1), cultural compatibility degree (C4,2), information platform compatibility degree (C4,3) and reputation (C4,4). C4,1, C4,2, C4,3 and C4,4 all belong to qualitative type.

3. Index Value Processing

For different types of indexes, different methods are used to get their values. To quantitative type index (i.e., C1,1 and C1,2), its value is obtained directly. To qualitative type index (i.e., C1,3, C1,4, C2,1, C2,2, C2,3 C3,1, C3,2, C3,3, C3,4, C4,1, C4,2, C4,3 and C4,4), its value, which is a score, is given by manager. If an index value can be accurately determined, it is a point value. If an index value is relatively fuzzy, it is an interval value. If an index value is completely unknown, it is a null value.
The suppliers to be assessed are x 1 , x 2 , , x M . For the supplier x r ( r = 1 ,   2 ,   ,   M ), the value on the index C j , l ( j = 1 ,   2 ,   ,   N and l = 1 ,   2 ,   ,   n j ) is represented as v r , ( j , l ) . Then, the normalized index value v r , ( j , l ) is calculated as follows. If the index belongs to benefit-type, v r , ( j , l ) = v r , ( j , l ) / max { v 1 , ( j , l ) , v 2 , ( j , l ) , , v M , ( j , l ) } . If the index belongs to cost-type, v r , ( j , l ) = min { v 1 , ( j , l ) , v 2 , ( j , l ) , , v M , ( j , l ) } / v r , ( j , l ) . Here, the interval index value is replaced with its left and right ends.
We set five comment levels which are very bad (G1), bad (G2), middle (G3), good (G4), very good (G5). Furtherly, G1 and G5 are the comment level corresponding to the lowest normalized index value v ( j , l ) ( G 1 ) = min { v 1 , ( j , l ) , v 2 , ( j , l ) , , v M , ( j , l ) } and the highest normalized index value v ( j , l ) ( G 5 ) = max { v 1 , ( j , l ) , v 2 , ( j , l ) , , v M , ( j , l ) } , respectively. Similarly, the interval index value is replaced with its left and right ends. Then, a number sequence v ( j , l ) ( G 1 ) , v ( j , l ) ( G 2 ) , v ( j , l ) ( G 3 ) , v ( j , l ) ( G 4 ) , v ( j , l ) ( G 5 ) is obtained.
It is assumed that the corresponding numbers of the five comment levels are π 1 = 0.1 ,   π 2 = 0.3 ,   π 3 = 0.5 ,   π 4 = 0.7 and π 5 = 0.9 . β u represents the membership degree of the index value to the comment level G u . On the index C j , l , the utility value of the supplier x r is represented as τ r , ( j , l ) . The normalized index value of the supplier x r is a point-value χ 1 , an interval value [ χ 1 , χ 2 ] or a null value. When v p χ 1 v p + 1 or v p χ 1 χ 2 v p + 1 ( p = 1 ,   2 ,   3 ,   4 ), τ r , ( j , l ) = β p π p + β p + 1 π p + 1 . When v p χ 1 v p + 1 and v p + 1 χ 2 v p + 2 ( p = 1 , 2 , 3 ), τ r , ( j , l ) = β p π p + β p + 1 π p + 1 + β p + 2 π p + 2 . When v p χ 1 v p + 1 and v o χ 2 v o + 1 ( p = 1 ,   2 ,   3 ,   4 , o = 1 ,   2 ,   3 ,   4 and o > p + 1 ), τ r , ( j , l ) = β p π p + β p + 1 π p + 1 + + β o π o + β o + 1 π o + 1 .

4. Index Weight Calculating

In ANP [11], the system elements are divided into two parts: (1) The first is called the control layer, including the problem objective and decision criteria. All decision criteria are considered independent of each other and are governed only by the problem objective. There can be no decision criteria in the control layer, but at least one objective; (2) The second part is the network layer, which is composed of all the elements that are controlled by the control layer, and its internal network structure is interacted.
Therefore, we set up the ANP structure of green supplier assessment as shown in Figure 2. The control layer only has one element: the green supplier assessment objective (AO), and the network layer has four element groups: product attribute (C1), internal competitiveness (C2), external competitiveness (C3) and cooperation ability (C4). Each element group affects each other and contains different elements. The elements in the same element group also affect each other. For example, internal competitiveness (C2) is affected by product attribute (C1), external competitiveness (C3) and cooperation ability (C4), and innovation capacity (C2,1), manufacturing capacity (C2,2) and agility capacity (C2,3) also affect each other.
According to the ANP structure of green supplier assessment shown in Figure 2, the control layer has the element AO and the network layer has the element groups C1, C2, …, CN (here, N = 4). The element group Ci (i = 1, 2, …, N) contains the elements C i , 1 , C i , 2 , , C i , n i . The control layer element AO is taken as the criterion and the element C j , l (l = 1, 2, …, n j ) in C j (j = 1, 2, …, N) is taken as the sub-criterion. Based on the influence of the elements in Ci on C j , l , the indirect dominance comparison of the elements in Ci are conducted.
Here, the influence of the elements in Ci on C j , l are assessed according to the personal experience and subjective judgment of experts, so using exact numbers to describe the influence of the elements in Ci on C j , l is unreasonable. In contrast, the fuzzy number can reflect the inherent uncertainty of the expert’s preference. At the same time, when integrating the opinions of multiple experts, the assessment of the influence of the elements in Ci on C j , l by experts is obviously with indiscernibility. The rough boundary interval in rough sets theory [26,27,28,29] can describe the indiscernibility as a set boundary area instead of a membership function, which can better integrate the assessment of multiple experts.
Based the index weight obtaining method in Reference [30], we design a rough ANP method to determine the index weight in green supplier assessment. The specific process of the designed rough ANP is as follows.
Step 1: Under the control layer element AO, we conduct the indirect dominance comparison of the elements in Ci according to their influence on C j , l .
There are q experts participating in the indirect dominance comparison of the elements in Ci. The fuzzy reciprocal judgement matrix E k , i , ( j , l ) = ( e g , h k , i , ( j , l ) ) n i × n i given by the expert k (k = 1, 2, …, q) is as follows:
E k , i , ( j , l ) = C j , l C i , 1 C i , 2 C i , n i C i , 1 1 e 1 , 2 k , i , ( j , l ) e 1 , n i k , i , ( j , l ) C i , 2 e 2 , 1 k , i , ( j , l ) 1 e 2 , n i k , i , ( j , l ) C i , n i e n i , 1 k , i , ( j , l ) e n i , 2 k , i , ( j , l ) 1
where e g , h k , i , ( j , l ) represents the indirect dominance score of the element C i , h compared to the element C i , g giver by the expert k, here g,h = 1, 2, …, ni and gh. e g , h k , i , ( j , l ) is a trapezoidal fuzzy number and e g , h k , i , ( j , l ) = ( a g , h k , i , ( j , l ) , b g , h k , i , ( j , l ) , c g , h k , i , ( j , l ) , d g , h k , i , ( j , l ) ) . a g , h k , i , ( j , l ) ,   b g , h k , i , ( j , l ) ,   c g , h k , i , ( j , l ) and d g , h k , i , ( j , l ) are all positive real numbers and a g , h k , i , ( j , l ) b g , h k , i , ( j , l ) c g , h k , i , ( j , l ) d g , h k , i , ( j , l ) . Then the consistency of the matrix given by each expert is verified. If it is qualified, the next step will be carried out; otherwise, this step will be returned.
Then, E k , i , ( j , l ) is split into A k , i , ( j , l ) = ( a g , h k , i , ( j , l ) ) n i × n i , B k , i , ( j , l ) = ( b g , h k , i , ( j , l ) ) n i × n i , C k , i , ( j , l ) = ( c g , h k , i , ( j , l ) ) n i × n i and D k , i , ( j , l ) = ( d g , h k , i , ( j , l ) ) n i × n i . Based on A 1 , i , ( j , l ) , A 2 , i , ( j , l ) , , A q , i , ( j , l ) , the rough group decision matrix A i , ( j , l ) = ( a g , h i , ( j , l ) ) n i × n i is constructed where a g , h i , ( j , l ) = { a g , h 1 , i , ( j , l ) , a g , h 2 , i , ( j , l ) , , a g , h q , i , ( j , l ) } , g,h = 1, 2, …, ni and gh.
The rough boundary interval of a g , h k , i , ( j , l ) a g , h i , ( j , l ) (k = 1, 2, …, q) is R N ( a g , h k , i , ( j , l ) ) = [ a g , h k , i , ( j , l ) , , a g , h k , i , ( j , l ) , + ] where a g , h k , i , ( j , l ) , , a g , h k , i , ( j , l ) , + are the rough lower limit and rough upper limit of a g , h k , i , ( j , l ) in the set a g , h i , ( j , l ) .
Thus, the rough boundary interval of a g , h i , ( j , l ) can be expressed as R N ( a g , h i , ( j , l ) ) = { [ a g , h 1 , i , ( j , l ) , , a g , h 1 , i , ( j , l ) , + ] , [ a g , h 2 , i , ( j , l ) , , a g , h 2 , i , ( j , l ) , + ] , , [ a g , h q , i , ( j , l ) , , a g , h q , i , ( j , l ) , + ] } . According to the calculation rule of rough boundary interval, the mean form of R N ( a g , h i , ( j , l ) ) can be obtained as A v g _ R N ( a g , h i , ( j , l ) ) = [ a g , h i , ( j , l ) , , a g , h i , ( j , l ) , + ] = [ k = 1 q a g , h k , i , ( j , l ) , / q , k = 1 q a g , h k , i , ( j , l ) , + / q ] where a g , h i , ( j , l ) , , a g , h i , ( j , l ) , + are the rough lower limit and rough upper limit of the set a g , h i , ( j , l ) .
The rough judgement matrix is constructed as E A i , ( j , l ) = ( A v g _ R N ( a g , h i , ( j , l ) ) ) n i × n i . Then, E A i , ( j , l ) is split into the rough lower limit matrix E A i , ( j , l ) , = ( a g , h i , ( j , l ) , ) n i × n i and the rough upper limit matrix E A i , ( j , l ) , + = ( a g , h i , ( j , l ) , + ) n i × n i .
The eigenvectors corresponding to the maximum eigenvalues of E A i , ( j , l ) , and E A i , ( j , l ) , + are V A i , ( j , l ) , = [ v a 1 i , ( j , l ) , , v a 2 i , ( j , l ) , , , v a n i i , ( j , l ) , ] T and V A i , ( j , l ) , + = [ v a 1 i , ( j , l ) , + , v a 2 i , ( j , l ) , + , , v a n i i , ( j , l ) , + ] T respectively, where v a h i , ( j , l ) , , v a h i , ( j , l ) , + are the value of V A i , ( j , l ) , and V A i , ( j , l ) , + on the h (h = 1, 2, …, ni) dimension. Then, a set G A i , ( j , l ) = { g a 1 i , ( j , l ) , g a 2 i , ( j , l ) , , g a n i i , ( j , l ) } can be obtained where g a h i , ( j , l ) = ( | v a h i , ( j , l ) , | + | v a h i , ( j , l ) , + | ) / 2 . Similarly, we also get the other sets: G B i , ( j , l ) = { g b 1 i , ( j , l ) , g b 2 i , ( j , l ) , , g b n i i , ( j , l ) } , G C i , ( j , l ) = { g c 1 i , ( j , l ) , g c 2 i , ( j , l ) , , g c n i i , ( j , l ) } and G D i , ( j , l ) = { g d 1 i , ( j , l ) , g d 2 i , ( j , l ) , , g d n i i , ( j , l ) } .
Thereupon, the eigenvector w i , ( j , l ) = [ w 1 i , ( j , l ) , w 2 i , ( j , l ) , , w n i i , ( j , l ) ] T is obtained, where w h i , ( j , l ) = [ ( g d h i , ( j , l ) ) 2 + g d h i , ( j , l ) g c h i , ( j , l ) + ( g c h i , ( j , l ) ) 2 ] [ ( g a h i , ( j , l ) ) 2 + g a h i , ( j , l ) g b h i , ( j , l ) + ( g b h i , ( j , l ) ) 2 ] 3 ( g d h i , ( j , l ) + g c h i , ( j , l ) g a h i , ( j , l ) g b h i , ( j , l ) ) , h = 1, 2, …, ni. Then, the eigenvector w i , ( j , l ) = [ w 1 i , ( j , l ) , w 2 i , ( j , l ) , , w n i i , ( j , l ) ] T is normalized as ω i , ( j , l ) = [ ω 1 i , ( j , l ) , ω 2 i , ( j , l ) , , ω n i i , ( j , l ) ] T where ω h i , ( j , l ) = w h i , ( j , l ) / h = 1 n i w h i , ( j , l ) .
Step 2: We represent Ω i , j = ( ω h i , ( j , l ) ) n i × n j as follows:
Ω i , j = [ ω 1 i , ( j , 1 ) ω 1 i , ( j , 2 ) ω 1 i , ( j , n j ) ω 2 i , ( j , 1 ) ω 2 i , ( j , 2 ) ω 2 i , ( j , n j ) ω n i i , ( j , 1 ) ω n i i , ( j , 2 ) ω n i i , ( j , n j ) ]
where the column vector ω i , ( j , l ) is the normalized influence degree sorting vector of the elements C i , 1 , C i , 2 , , C i , n i in C i on the element C j , l in C j . If the elements in C j is not affected by the elements in C i , Ω i , j = 0 .
So we get the hyper-matrix Ω under the control layer element AO as follows:
C 1 C 2 C N C 1 , 1 C 1 , n 1 C 2 , 1 C 2 , n 2 C N , 1 C N , n N Ω = C 1 , 1 C 1 C 1 , n 1 C 2 , 1 C 2 C 2 , n 2 C N , 1 C N C N , n N [      Ω 1 , 1          Ω 1 , 2               Ω 1 , N          Ω 2 , 1          Ω 2 , 2               Ω 2 , N                              Ω N , 1          Ω N , 2               Ω N , N    ]
Step 3: The sub-block Ω i , j of Ω is column normalized, but Ω isn’t column normalized. To solve this problem, we conduct the indirect dominance comparison of the element groups C 1 , C 2 , , C N according to their influence on C j ( j = 1 ,   2 ,   ,   N ) under the control layer element AO. Here, we adopt a similar approach to Step 1 and get the relative importance matrix Ψ = ( ψ i , j ) N × N of element groups as follows:
Ψ = [ ψ 1 , 1 ψ 1 , 2 ψ 1 , N ψ 2 , 1 ψ 2 , 2 ψ 2 , N ψ N , 1 ψ N , 2 ψ N , N ]
where the column vector ψ , j = [ ψ 1 , j , ψ 2 , j , , ψ N , j ] T is the normalized influence degree sorting vector of the element groups C 1 , C 2 , , C N on C j .
The weighted form of the hyper-matrix Ω is Ω ¯ , where Ω ¯ i , j = ψ i , j Ω i , j ( i = 1 , 2 , , N , j = 1 , 2 , , N ) .
Step 4: We do the square operation of the weighted hyper-matrix Ω ¯ until the result converges to a stable limit hyper-matrix Ω ¯ as follows:
Ω ¯ = lim t Ω ¯ t
Any column of Ω ¯ is the limit relative ranking vector of all elements in the network layer. So we can get the weight vector of the index system shown in Figure 1 as follows:
θ = [ θ 1 , 1 , θ 1 , 2 , , θ 1 , n 1 , θ 2 , 1 , θ 2 , 2 , , θ 2 , n 2 , θ 3 , 1 , θ 3 , 2 , , θ 3 , n 3 , θ 4 , 1 , θ 4 , 2 , , θ 4 , n 4 ] T
where θ j , 1 , θ j , 2 , , θ j , n j are the weights of C j , 1 , C j , 2 , , C j , n j , j = 1 ,   2 ,   ,   N .
Furtherly, the weight of C j ( j = 1 ,   2 ,   ,   N ) is obtained as follows:
θ j = l = 1 n j θ j , l

5. Green Supplier Assessment

5.1. Related Concepts of Evidence Theory

We assume that there are M suppliers to be assessed. Based on evidence theory [31,32], the set of the suppliers to be assessed is defined as the identification framework Θ = { x 1 , x 2 , , x M } . All possible sets in Θ are represented by the power set 2 Θ . If each element in Θ is incompatible with each other, the number of elements in 2 Θ is 2 M . Then, a set function m a s s : 2 Θ [ 0 , 1 ] , which satisfies m a s s ( ϕ ) = 0 and φ Θ m a s s ( φ ) = 1 , is defined. The set function m a s s is known as the basic probability distribution function on Θ . Here, φ represents a supplier to be assessed. m a s s ( φ ) is the basic probability distribution value of φ and represents the trust degree for Φ , Any supplier to be assessed satisfying the condition “ m a s s ( φ ) > 0 ” is called a focal element.
For φ , φ 1 , φ 2 , , φ n Θ , the fusion rule of the basic probability distribution functions m a s s 1 , m a s s 2 , , m a s s n on Θ is as follows:
m a s s ( φ ) = 1 K φ 1 φ 2 φ n = m a s s 1 ( φ 1 ) m a s s 2 ( φ 2 ) m a s s n ( φ n )
where m a s s = m a s s 1 m a s s 2 m a s s n .
The normalization constant K is defined as follows:
K = φ 1 φ 2 φ n m a s s 1 ( φ 1 ) m a s s 2 ( φ 2 ) m a s s n ( φ n )
The total trust degree of φ on Θ can be expressed as the belief function B e l ( φ ) = ϕ Θ , ϕ φ m a s s ( ϕ ) , and the uncertainty degree of φ on Θ can be expressed as the plausible function P l ( φ ) = ϕ Θ , ϕ φ m a s s ( ϕ ) . For a supplier φ on Θ , B e l ( φ ) represents the sum of the possibility measurement for all subsets of φ , and P l ( φ ) represents the sum of the uncertainty measurement for all subsets of φ . The confirmation degree of φ is represented by the trust interval [ B e l ( φ ) , P l ( φ ) ] .
Therefore, B e l ( φ ) reflects the sum of exact reliability which the evidences support φ , and P l ( φ ) reflects the sum of reliability which the evidences do not negate φ . As a result, B e l ( φ ) and P l ( φ ) can be considered as the minimum and maximum probability bounds respectively, so [ B e l ( φ ) , P l ( φ ) ] can form the trust interval.
Based on the above analysis, assessing the suppliers by trust interval is more reliable than by maximum belief function or maximum plausible function [30,31]. We assume that the supplier φ s is better than supplier φ y with a degree of ε s , y ( 0 ε s , y 1 ). If the trust intervals of φ s and φ y are [ B e l ( φ s ) , P l ( φ s ) ] and [ B e l ( φ y ) , P l ( φ y ) ] respectively, ε s , y is defined as follows:
ε s , y = max { 0 , P l ( φ s ) B e l ( φ y ) } max { 0 , B e l ( φ s ) P l ( φ y ) } ( P l ( φ s ) B e l ( φ s ) ) + ( P l ( φ y ) B e l ( φ y ) )
The decision rules based on trust interval are as follows:
  • If ε s , y > 0.5 , φ s is better than φ y (recorded as φ s φ y );
  • If ε s , y < 0.5 , φ s is worse than φ y (recorded as φ s φ y );
  • If ε s , y = 0.5 , φ s and φ y are equal (recorded as φ s = φ y );
  • For three suppliers φ s , φ y and φ z , if ε s , y > 0.5 and ε y , z > 0.5 , φ s is better than φ y and φ y is better than φ z , so φ s φ y φ z .

5.2. Green Supplier Assessing Procedure

On the index C j , l , the weighted basic probability distribution value of the focal element φ s (s < 2N) is represented as m a s s C j , l ( φ s ) . In this paper, we use m a s s C j , l ( φ s ) as the evidence input of green supplier assessment.
The utility value of each focal element except Θ can be calculated through index value processing. The special focal element Θ can indicate the uncertainty of the expert on an index. If we don’t consider the influence of Θ , the green supplier assessment problem will become a simple probability distribution problem, and the advantages of evidence theory will not be applied. Meanwhile, the trust degree of expert in each index is different and the uncertainty of index is reflected by the probability distribution of Θ . Thus the probability distribution value of Θ on different indexes should also be treated differently.
In the green supplier assessment problem, the index weight is obviously not fixed in the case of different requirements. When costs need to be reduced, C11 is more important than other indexes and its weight must be higher than other indexes, and the basic probability distribution value of Θ on C11 should be smaller than other indexes. Therefore, the index weight calculated by rough ANP is introduced to adjust the preference of experts and solve the probability distribution problem of Θ on different indexes. Then, the weighted basic probability distribution value of each focal element on each index is obtained as m a s s C j , l ( φ s ) .
We take a weighted normalization treatment for the basic probability distribution values of all focal elements and calculate m a s s C j , l ( φ s ) , as follows:
{ m a s s C j , l ( φ s ) = θ j , l τ s , ( j , l ) s = 1 l 1 τ s , ( j , l ) , φ s Θ m a s s C j , l ( φ s ) = 1 θ j , l , φ s = Θ
Based on the two-level index system shown in Figure 1, we establish a two-order fusion evidence theory model for green supplier assessment. The procedure is shown in Figure 3.

6. Case Study

6.1. An Application Case of Bearing Cage Supplier Assessment

For a bearing manufacturing enterprise, there are three bearing cage suppliers to be assessed. The set of suppliers to be assessed is { φ 1 , φ 2 , φ 3 } . The best bearing cage supplier need to be selected after green supplier assessment. The index value of quantitative type (C11 and C12) is obtained directly from the enterprise resources planning system (ERP) of the bearing manufacturing enterprise, and the index value of qualitative type (other indexes) is obtained by the method of expert’s scoring (i.e., 0.1, 0.3, 0.5, 0.7 and 0.9). The index values of three bearing cage suppliers are shown in Table 1. The units of the index value on C 1 , 1 and C 1 , 2 are RMB and mm (error value) respectively.
Then, the normalized index values of three bearing cage suppliers are shown in Table 2.
Five comment levels: very bad (G1), bad (G2), middle (G3), good (G4) and very good (G5) are set. Taking the normalized index values v 1 , ( 1 , 1 ) = 0.0023 , v 2 , ( 1 , 1 ) = 0.0669 and v 3 , ( 1 , 1 ) = 1.0000 for an example, we can get v ( 1 , 1 ) ( G 1 ) = 0.0023 and v ( 1 , 1 ) ( G 5 ) = 1.0000 , so v ( 1 , 1 ) ( G 2 ) = 0.2517 , v ( 1 , 1 ) ( G 3 ) = 0.5011 and v ( 1 , 1 ) ( G 4 ) = 0.7506 . The corresponding numbers of the five comment levels are π 1 = 0.1 ,   π 2 = 0.3 ,   π 3 = 0.5 ,   π 4 = 0.7 and π 5 = 0.9 . To the normalized index value v 2 , ( 1 , 1 ) = 0.0669 , the membership degrees are β 1 = 0.7410 and β 2 = 0.2590 , so τ 2 , ( 1 , 1 ) = β 1 π 1 + β 2 π 2 = 0.0648 . The utility values of three bearing cage suppliers are shown in Table 3.
Next, the index weight is calculated by rough ANP. There are four experts: expert 1, expert 2, expert 3 and expert 4.
Taking the indirect dominance comparison of the elements in C 1 according to their influence on C 2 , 1 under the control layer element AO for an example, the fuzzy reciprocal judgement matrices given by the four experts are shown as follows:
E 1 , 1 , ( 2 , 1 ) = C 2 , 1 C 1 , 1 C 1 , 2 C 1 , 3 C 1 , 4 C 1 , 1 ( 1 , 1 , 1 , 1 ) ( 1 , 11 / 9 , 13 / 7 , 7 / 3 ) ( 1 , 11 / 9 , 13 / 7 , 7 / 3 ) ( 3 / 2 , 13 / 7 , 3 , 4 ) C 1 , 2 ( 3 / 7 , 7 / 13 , 9 / 11 , 1 ) ( 1 , 1 , 1 , 1 ) ( 1 , 11 / 9 , 13 / 7 , 7 / 3 ) ( 1 , 11 / 9 , 13 / 7 , 7 / 3 ) C 1 , 3 ( 3 / 7 , 7 / 13 , 9 / 11 , 1 ) ( 3 / 7 , 7 / 13 , 9 / 11 , 1 ) ( 1 , 1 , 1 , 1 ) ( 1 , 11 / 9 , 13 / 7 , 7 / 3 ) C 1 , 4 ( 1 / 4 , 1 / 3 , 7 / 13 , 2 / 3 ) ( 3 / 7 , 7 / 13 , 9 / 11 , 1 ) ( 3 / 7 , 7 / 13 , 9 / 11 , 1 ) ( 1 , 1 , 1 , 1 ) ,
E 2 , 1 , ( 2 , 1 ) = C 2 , 1 C 1 , 1 C 1 , 2 C 1 , 3 C 1 , 4 C 1 , 1 ( 1 , 1 , 1 , 1 ) ( 1 , 1 , 1 , 1 ) ( 1 , 11 / 9 , 13 / 7 , 7 / 3 ) ( 1 , 11 / 9 , 13 / 7 , 7 / 3 ) C 1 , 2 ( 1 , 1 , 1 , 1 ) ( 1 , 1 , 1 , 1 ) ( 1 , 11 / 9 , 13 / 7 , 7 / 3 ) ( 1 , 11 / 9 , 13 / 7 , 7 / 3 ) C 1 , 3 ( 3 / 7 , 7 / 13 , 9 / 11 , 1 ) ( 3 / 7 , 7 / 13 , 9 / 11 , 1 ) ( 1 , 1 , 1 , 1 ) ( 1 , 1 , 1 , 1 ) C 1 , 4 ( 3 / 7 , 7 / 13 , 9 / 11 , 1 ) ( 3 / 7 , 7 / 13 , 9 / 11 , 1 ) ( 1 , 1 , 1 , 1 ) ( 1 , 1 , 1 , 1 ) ,
E 3 , 1 , ( 2 , 1 ) = C 2 , 1 C 1 , 1 C 1 , 2 C 1 , 3 C 1 , 4 C 1 , 1 ( 1 , 1 , 1 , 1 ) ( 3 / 2 , 13 / 7 , 3 , 4 ) ( 1 , 1 , 1 , 1 ) ( 1 , 11 / 9 , 13 / 7 , 7 / 3 ) C 1 , 2 ( 1 / 4 , 1 / 3 , 7 / 13 , 2 / 3 ) ( 1 , 1 , 1 , 1 ) ( 1 / 4 , 1 / 3 , 7 / 13 , 2 / 3 ) ( 3 / 7 , 7 / 13 , 9 / 11 , 1 ) C 1 , 3 ( 1 , 1 , 1 , 1 ) ( 3 / 2 , 13 / 7 , 3 , 4 ) ( 1 , 1 , 1 , 1 ) ( 1 , 11 / 9 , 13 / 7 , 7 / 3 ) C 1 , 4 ( 3 / 7 , 7 / 13 , 9 / 11 , 1 ) ( 1 , 11 / 9 , 13 / 7 , 7 / 3 ) ( 3 / 7 , 7 / 13 , 9 / 11 , 1 ) ( 1 , 1 , 1 , 1 ) ,
E 4 , 1 , ( 2 , 1 ) = C 2 , 1 C 1 , 1 C 1 , 2 C 1 , 3 C 1 , 4 C 1 , 1 ( 1 , 1 , 1 , 1 ) ( 1 , 11 / 9 , 13 / 7 , 7 / 3 ) ( 3 / 2 , 13 / 7 , 3 , 4 ) ( 7 / 3 , 3 , 17 / 3 , 9 ) C 1 , 2 ( 3 / 7 , 7 / 13 , 9 / 11 , 1 ) ( 1 , 1 , 1 , 1 ) ( 1 , 11 / 9 , 13 / 7 , 7 / 3 ) ( 3 / 2 , 13 / 7 , 3 , 4 ) C 1 , 3 ( 1 / 4 , 1 / 3 , 7 / 13 , 2 / 3 ) ( 3 / 7 , 7 / 13 , 9 / 11 , 1 ) ( 1 , 1 , 1 , 1 ) ( 1 , 11 / 9 , 13 / 7 , 7 / 3 ) C 1 , 4 ( 1 / 9 , 3 / 17 , 1 / 3 , 3 / 7 ) ( 1 / 4 , 1 / 3 , 7 / 13 , 2 / 3 ) ( 3 / 7 , 7 / 13 , 9 / 11 , 1 ) ( 1 , 1 , 1 , 1 ) .
We check the consistency of E 1 , 1 , ( 2 , 1 ) ,   E 2 , 1 , ( 2 , 1 ) ,   E 3 , 1 , ( 2 , 1 ) and E 4 , 1 , ( 2 , 1 ) and all of them are qualified.
Then, E k , 1 , ( 2 , 1 ) ( k = 1 ,   2 ,   3 ,   4 ) is split into A k , 1 , ( 2 , 1 ) = ( a g , h k , 1 , ( 2 , 1 ) ) 4 × 4 , B k , 1 , ( 2 , 1 ) = ( b g , h k , 1 , ( 2 , 1 ) ) 4 × 4 , C k , 1 , ( 2 , 1 ) = ( c g , h k , 1 , ( 2 , 1 ) ) 4 × 4 and D k , 1 , ( 2 , 1 ) = ( d g , h k , 1 , ( 2 , 1 ) ) 4 × 4 . For example, A 1 , 1 , ( 2 , 1 ) = ( a g , h 1 , 1 , ( 2 , 1 ) ) 4 × 4 is as follows:
A 1 , 1 , ( 2 , 1 ) = [ 1 1 1 3 / 2 3 / 7 1 1 1 3 / 7 3 / 7 1 1 1 / 4 3 / 7 3 / 7 1 ] .
Based on A 1 , 1 , ( 2 , 1 ) ,   A 2 , 1 , ( 2 , 1 ) ,   A 3 , 1 , ( 2 , 1 ) ,   A 4 , 1 , ( 2 , 1 ) , the rough group decision matrix A 1 , ( 2 , 1 ) = ( a g , h 1 , ( 2 , 1 ) ) 4 × 4 is constructed as follows:
A 1 , ( 2 , 1 ) = [ { 1 , 1 , 1 , 1 } { 1 , 1 , 3 / 2 , 1 } { 1 , 1 , 1 , 3 / 2 } { 3 / 2 , 1 , 1 , 7 / 3 } { 3 / 7 , 1 , 1 / 4 , 3 / 7 } { 1 , 1 , 1 , 1 } { 1 , 1 , 1 / 4 , 1 } { 1 , 1 , 3 / 7 , 3 / 2 } { 3 / 7 , 3 / 7 , 1 , 1 / 4 } { 3 / 7 , 3 / 7 , 3 / 2 , 3 / 7 } { 1 , 1 , 1 , 1 } { 1 , 1 , 1 , 1 } { 1 / 4 , 3 / 7 , 3 / 7 , 1 / 9 } { 3 / 7 , 3 / 7 , 1 , 1 / 4 } { 3 / 7 , 1 , 3 / 7 , 3 / 7 } { 1 , 1 , 1 , 1 } ] .
For the partition a 1 , 4 1 , 1 , ( 2 , 1 ) = 3/2 in the element a 1 , 4 1 , ( 2 , 1 ) = { 3 / 2 , 1 , 1 , 7 / 3 } , its upper approximation set is {3/2, 7/3} and lower approximation set is {3/2, 1, 1}, then L _ ( a 1 , 4 1 , 1 , ( 2 , 1 ) ) = (3/2 + 1 + 1)/3 = 1.17, L ¯ ( a 1 , 4 1 , 1 , ( 2 , 1 ) ) = (3/2 + 7/3)/2 = 1.92 and R N ( a 1 , 4 1 , 1 , ( 2 , 1 ) ) = [1.17, 1.92]. Similarly, R N ( a 1 , 4 2 , 1 , ( 2 , 1 ) ) = R N ( a 1 , 4 3 , 1 , ( 2 , 1 ) ) = [1, 1.46], R N ( a 1 , 4 4 , 1 , ( 2 , 1 ) ) = [1.46, 2.33]. So R N ( a 1 , 4 1 , ( 2 , 1 ) ) = {[1.17, 1.92], [1, 1.46], [1, 1.46], [1.46, 2.33]}.
According to the calculation rule of rough boundary interval, the mean form of R N ( a 1 , 4 1 , ( 2 , 1 ) ) is obtained as A v g _ R N ( a 1 , 4 1 , ( 2 , 1 ) ) = [1.16, 1.79]. Then, the rough judgement matrix E A 1 , ( 2 , 1 ) = ( A v g _ R N ( a g , h 1 , ( 2 , 1 ) ) ) 4 × 4 is constructed as follows:
E A 1 , ( 2 , 1 ) = [ [ 1 , 1 ] [ 1.03 , 1.22 ] [ 1.03 , 1.22 ] [ 1.16 , 1.79 ] [ 0.38 , 0.74 ] [ 1 , 1 ] [ 0.67 , 0.95 ] [ 0.72 , 1.17 ] [ 0.38 , 0.74 ] [ 0.50 , 0.90 ] [ 1 , 1 ] [ 1 , 1 ] [ 0.23 , 0.38 ] [ 0.38 , 0.74 ] [ 0.46 , 0.68 ] [ 1 , 1 ] ] .
E A 1 , ( 2 , 1 ) is split into the rough lower limit matrix E A 1 , ( 2 , 1 ) , = ( a g , h 1 , ( 2 , 1 ) , ) 4 × 4 and the rough upper limit matrix E A 1 , ( 2 , 1 ) , + = ( a g , h 1 , ( 2 , 1 ) , + ) 4 × 4 . The eigenvectors corresponding to the maximum eigenvalues of E A 1 , ( 2 , 1 ) , and E A 1 , ( 2 , 1 ) , + are V A 1 , ( 2 , 1 ) , = [0.71, 0.44, 0.45, 0.30]T and V A 1 , ( 2 , 1 ) , + = [0.65, 0.49, 0.47, 0.34]T respectively. Then, we get a set G A 1 , ( 2 , 1 ) = {0.68, 0.47, 0.46, 0.32}. Similarly, we also get the other sets: G B 1 , ( 2 , 1 ) = {0.73, 0.51, 0.66, 0.58}, G C 1 , ( 2 , 1 ) = {0.82, 0.67, 0.73, 0.69} and G D 1 , ( 2 , 1 ) = {0.95, 0.77, 0.83, 0.75}.
Thereupon, we obtain the normalized eigenvector ω 1 , ( 2 , 1 ) = [0.30, 0.23, 0.25, 0.22]T.
Similarly, we obtain the normalized eigenvector ω 1 , ( 2 , 2 ) = [0.28, 0.41, 0.17, 0.14]T and ω 1 , ( 2 , 3 ) = [0.33, 0.34, 0.13, 0.20]T. So we obtain Ω 1 , 2 = ( ω h 1 , ( 2 , l ) ) 4 × 4 as follows:
Ω 1 , 2 = [ 0.30 0.28 0.33 0.23 0.41 0.34 0.25 0.17 0.13 0.22 0.14 0.20 ] .
After the similar calculating, we get the hyper-matrix Ω under the control layer element G as follows:
Ω = [ 0.31 0.19 0.22 0.33 0.30 0.28 0.33 0.18 0.30 0.41 0.38 0.22 0.19 0.33 0.40 0.24 0.25 0.33 0.31 0.23 0.41 0.34 0.32 0.33 0.22 0.19 0.31 0.17 0.18 0.10 0.29 0.25 0.24 0.18 0.25 0.17 0.13 0.23 0.10 0.20 0.22 0.40 0.33 0.27 0.29 0.16 0.31 0.21 0.18 0.22 0.14 0.20 0.27 0.27 0.17 0.21 0.07 0.31 0.22 0.21 0.24 0.35 0.48 0.20 0.13 0.11 0.51 0.43 0.29 0.35 0.36 0.28 0.51 0.40 0.35 0.25 0.30 0.31 0.41 0.29 0.40 0.24 0.16 0.13 0.45 0.40 0.33 0.30 0.49 0.35 0.51 0.35 0.21 0.39 0.58 0.49 0.25 0.41 0.58 0.20 0.24 0.39 0.19 0.11 0.30 0.33 0.28 0.31 0.25 0.17 0.29 0.18 0.10 0.33 0.30 0.15 0.44 0.32 0.30 0.33 0.32 0.20 0.25 0.31 0.44 0.20 0.25 0.30 0.18 0.19 0.25 0.18 0.28 0.39 0.33 0.25 0.28 0.09 0.25 0.16 0.20 0.32 0.27 0.18 0.31 0.24 0.26 0.30 0.11 0.20 0.10 0.24 0.35 0.19 0.23 0.31 0.25 0.33 0.31 0.20 0.36 0.12 0.10 0.20 0.14 0.22 0.30 0.18 0.41 0.15 0.29 0.19 0.33 0.23 0.25 0.28 0.38 0.41 0.25 0.33 0.28 0.20 0.41 0.13 0.30 0.17 0.22 0.36 0.18 0.19 0.30 0.21 0.12 0.28 0.16 0.40 0.19 0.19 0.17 0.30 0.21 0.08 0.17 0.22 0.25 0.10 0.19 0.17 0.19 0.31 0.10 0.31 0.22 0.29 0.25 0.33 0.51 0.14 0.37 0.31 0.32 0.22 0.30 0.28 0.20 ] .
Then we conduct the indirect dominance comparison of the element groups C 1 , C 2 , C 3 , C 4 according to their influence on C j ( j = 1 ,   2 ,   3 ,   4 ) under the control layer element AO. As a result, we get the relative importance matrix Ψ = ( ψ i , j ) 4 × 4 of element groups as follows:
Ψ = [ 0.51 0.38 0.42 0.29 0.12 0.28 0.20 0.21 0.23 0.19 0.15 0.10 0.14 0.15 0.23 0.40 ] .
The weighted form of the hyper-matrix Ω is Ω ¯ , where Ω ¯ i , j = ψ i , j Ω i , j ( i = 1 , 2 , 3 , 4 and j = 1 , 2 , 3 , 4 ). The square operation of the weighted hyper-matrix Ω ¯ is done continuously and Ω ¯ 4 converges to a stable limit hyper-matrix. So we get the weight vector of the index system shown in Figure 1 as follows: θ = [0.12, 0.11, 0.10, 0.09, 0.06, 0.06, 0.07, 0.05, 0.04, 0.04, 0.04, 0.06, 0.05, 0.05, 0.06]T. In the case of ensuring the weight proportional relationship among C 1 , 1 , C 1 , 2 , C 1 , 3 , C 1 , 4 , we expand the weight vector of C 1 , 1 , C 1 , 2 , C 1 , 3 , C 1 , 4 to [0.95, 0.87, 0.79, 0.71]T. Similarly, the weight vector of C 2 , 1 , C 2 , 2 , C 2 , 3 is expanded to [0.81, 0.81, 0.95]T, the weight vector of C 3 , 1 , C 3 , 2 , C 3 , 3 , C 3 , 4 is expanded to [0.95, 0.76, 0.76, 0.76]T, and the weight vector of C 4 , 1 , C 4 , 2 , C 4 , 3 , C 4 , 4 is expanded to [0.95, 0.79, 0.79, 0.95]T. Furtherly, the weight vector of C 1 , C 2 , C 3 , C 4 is obtained as [0.42, 0.19, 0.17, 0.22]T. The relative weight vector of the index C 1 , C 2 , C 3 , C 4 is [0.95, 0.43, 0.38, 0.50]T.
Based on evidence theory, the set of the three suppliers to be assessed is defined as the identification framework: Θ = { φ 1 ,   φ 2 ,   φ 3 } . According to Equation (11), the weighted basic probability distribution values of all focal elements are calculated based on the utility values of three bearing cage suppliers shown in Table 3 and the above relative weight vector as follows:
(1)
m a s s C 1 , 1 ( φ 1 ) = 0.0825, m a s s C 1 , 1 ( φ 2 ) = 0.1252, m a s s C 1 , 1 ( φ 3 ) = 0.7423, m a s s C 1 , 1 ( Θ ) = 0.0500;
(2)
m a s s C 1 , 2 ( φ 1 ) = 0.4121, m a s s C 1 , 2 ( φ 2 ) = 0.4121, m a s s C 1 , 2 ( φ 3 ) = 0.0458, m a s s C 1 , 2 ( Θ ) = 0.1300;
(3)
m a s s C 1 , 3 ( φ 1 ) = 0.3456, m a s s C 1 , 3 ( φ 2 ) = 0.3999, m a s s C 1 , 3 ( φ 3 ) = 0.0444, m a s s C 1 , 3 ( Θ ) = 0.2100;
(4)
m a s s C 1 , 4 ( φ 1 ) = 0.0374, m a s s C 1 , 4 ( φ 2 ) = 0.3363, m a s s C 1 , 4 ( φ 3 ) = 0.3363, m a s s C 1 , 4 ( Θ ) = 0.2900;
(5)
m a s s C 2 , 1 ( φ 1 ) = 0.1272, m a s s C 2 , 1 ( φ 2 ) = 0.2820, m a s s C 2 , 1 ( φ 3 ) = 0.4008, m a s s C 2 , 1 ( Θ ) = 0.1900;
(6)
m a s s C 2 , 2 ( φ 1 ) = 0.3141, m a s s C 2 , 2 ( φ 2 ) = 0.4463, m a s s C 2 , 2 ( φ 3 ) = 0.0496, m a s s C 2 , 2 ( Θ ) = 0.1900;
(7)
m a s s C 2 , 3 ( φ 1 ) = 0.3684, m a s s C 2 , 3 ( φ 2 ) = 0.5234, m a s s C 2 , 3 ( φ 3 ) = 0.0582, m a s s C 2 , 3 ( Θ ) = 0.0500;
(8)
m a s s C 3 , 1 ( φ 1 ) = 0.8550, m a s s C 3 , 1 ( φ 3 ) = 0.0950, m a s s C 3 , 1 ( Θ ) = 0.0500;
(9)
m a s s C 3 , 2 ( φ 1 ) = 0.0691, m a s s C 3 , 2 ( φ 2 ) = 0.6218, m a s s C 3 , 2 ( φ 3 ) = 0.0691, m a s s C 3 , 2 ( Θ ) = 0.2400;
(10)
m a s s C 3 , 3 ( φ 1 ) = 0.3683, m a s s C 3 , 3 ( φ 2 ) = 0.3508, m a s s C 3 , 3 ( φ 3 ) = 0.0409, m a s s C 3 , 3 ( Θ ) = 0.2400;
(11)
m a s s C 3 , 4 ( φ 1 ) = 0.2947, m a s s C 3 , 4 ( φ 2 ) = 0.4188, m a s s C 3 , 4 ( φ 3 ) = 0.0465, m a s s C 3 , 4 ( Θ ) = 0.2400;
(12)
m a s s C 4 , 1 ( φ 2 ) = 0.8550, m a s s C 4 , 1 ( φ 3 ) = 0.0950, m a s s C 4 , 1 ( Θ ) = 0.0500;
(13)
m a s s C 4 , 2 ( φ 1 ) = 0.6449, m a s s C 4 , 2 ( φ 2 ) = 0.0726, m a s s C 4 , 2 ( φ 3 ) = 0.0726, m a s s C 4 , 2 ( Θ ) = 0.2100;
(14)
m a s s C 4 , 3 ( φ 1 ) = 0.0790, m a s s C 4 , 3 ( φ 2 ) = 0.7110, m a s s C 4 , 3 ( Θ ) = 0.2100;
(15)
m a s s C 4 , 4 ( φ 1 ) = 0.3912, m a s s C 4 , 4 ( φ 2 ) = 0.5029, m a s s C 4 , 4 ( φ 3 ) = 0.0559, m a s s C 4 , 4 ( Θ ) = 0.0500.
Then we make { m a s s C 1 , 1 ( φ s ) ,   m a s s C 1 , 2 ( φ s ) ,   m a s s C 1 , 3 ( φ s ) ,   m a s s C 1 , 4 ( φ s ) } , { m a s s C 2 , 1 ( φ s ) ,   m a s s C 2 , 2 ( φ s ) ,   m a s s C 2 , 3 ( φ s ) } , { m a s s C 3 , 1 ( φ s ) ,   m a s s C 3 , 2 ( φ s ) ,   m a s s C 3 , 3 ( φ s ) ,   m a s s C 3 , 4 ( φ s ) } and { m a s s C 4 , 1 ( φ s ) ,   m a s s C 4 , 2 ( φ s ) ,   m a s s C 4 , 3 ( φ s ) ,   m a s s C 4 , 4 ( φ s ) } as the evidence input and execute the first evidence fusion respectively. Here, m a s s C 1 = m a s s C 1 , 1 m a s s C 1 , 2 m a s s C 1 , 3 m a s s C 1 , 4 , m a s s C 2 = m a s s C 2 , 1 m a s s C 2 , 2 m a s s C 2 , 3 , m a s s C 3 = m a s s C 3 , 1 m a s s C 3 , 2 m a s s C 3 , 3 m a s s C 3 , 4 , m a s s C 4 = m a s s C 4 , 1 m a s s C 4 , 2 m a s s C 4 , 3 m a s s C 4 , 4 . The basic probability distribution values m a s s C 1 ( φ s ) ,   m a s s C 2 ( φ s ) ,   m a s s C 3 ( φ s ) ,   m a s s C 4 ( φ s ) of all focal elements are calculated as follows:
(1)
m a s s C 1 ( φ 1 ) = 0.0491, m a s s C 1 ( φ 2 ) = 0.6871, m a s s C 1 ( φ 3 ) = 0.1938, m a s s C 1 ( Θ ) = 0.0700;
(2)
m a s s C 2 ( φ 1 ) = 0.2305, m a s s C 2 ( φ 2 ) = 0.7000, m a s s C 2 ( φ 3 ) = 0.0345, m a s s C 2 ( Θ ) = 0.0350;
(3)
m a s s C 3 ( φ 1 ) = 0.8165, m a s s C 3 ( φ 2 ) = 0.1122, m a s s C 3 ( φ 3 ) = 0.0033, m a s s C 3 ( Θ ) = 0.0679;
(4)
m a s s C 4 ( φ 1 ) = 0.0104, m a s s C 4 ( φ 2 ) = 0.9841, m a s s C 4 ( φ 3 ) = 0.0008, m a s s C 4 ( Θ ) = 0.0046.
With the consideration of the relative weight vector of the index C 1 ,   C 2 ,   C 3 ,   C 4 , we normalize the basic probability distribution values { m a s s C 1 ( φ 1 ) ,   m a s s C 1 ( φ 2 ) ,   m a s s C 1 ( φ 3 ) ,   m a s s C 1 ( Θ ) } , { m a s s C 2 ( φ 1 ) ,   m a s s C 2 ( φ 2 ) ,   m a s s C 2 ( φ 3 ) ,   m a s s C 2 ( Θ ) } , { m a s s C 3 ( φ 1 ) ,   m a s s C 3 ( φ 2 ) ,   m a s s C 3 ( φ 3 ) ,   m a s s C 3 ( Θ ) } and { m a s s C 3 ( φ 1 ) ,   m a s s C 3 ( φ 2 ) ,   m a s s C 3 ( φ 3 ) ,   m a s s C 3 ( Θ ) } . The weighted basic probability distribution values m a s s C 1 ( φ s ) ,   m a s s C 2 ( φ s ) ,   m a s s C 3 ( φ s ) ,   m a s s C 4 ( φ s ) of all focal elements are calculated as follows:
(1)
m a s s C 1 ( φ 1 ) = 0.0466, m a s s C 1 ( φ 2 ) = 0.6527, m a s s C 1 ( φ 3 ) = 0.1841, m a s s C 1 ( Θ ) = 0.1165;
(2)
m a s s C 2 ( φ 1 ) = 0.0991, m a s s C 2 ( φ 2 ) = 0.3010, m a s s C 2 ( φ 3 ) = 0.0148, m a s s C 2 ( Θ ) = 0.5851;
(3)
m a s s C 3 ( φ 1 ) = 0.3103, m a s s C 3 ( φ 2 ) = 0.0426, m a s s C 3 ( φ 3 ) = 0.0013, m a s s C 3 ( Θ ) = 0.6458;
(4)
m a s s C 4 ( φ 1 ) = 0.0052, m a s s C 4 ( φ 2 ) = 0.4921, m a s s C 4 ( φ 3 ) = 0.0004, m a s s C 4 ( Θ ) = 0.5023.
Then, we make m a s s C 1 ( φ s ) ,   m a s s C 2 ( φ s ) ,   m a s s C 3 ( φ s ) ,   m a s s C 4 ( φ s ) m ˜ 4 ( A i ) as the evidence input and execute the second evidence fusion. The basic probability distribution values m a s s ( φ s ) are calculated as follows:
(1)
m a s s ( φ 1 ) = 0.0042;
(2)
m a s s ( φ 2 ) = 0.4160;
(3)
m a s s ( φ 3 ) = 0.0006;
(4)
m a s s ( Θ ) = 0.5792;
Therefore, the belief function B e l ( φ s ) and plausible function P l ( φ s ) of the three suppliers are calculated, and then the trust interval [ B e l ( φ s ) ,   P l ( φ s ) ] are obtained as follows:
(1)
[ B e l ( φ 1 ) ,   P l ( φ 1 ) ] = [0.0042, 0.5834];
(2)
[ B e l ( φ 2 ) ,   P l ( φ 2 ) ] = [0.4160, 0.9952];
(3)
[ B e l ( φ 3 ) ,   P l ( φ 3 ) ] = [0.0006, 0.5798].
According to the decision rules based on trust interval in Section 5.1, we obtain the results as follows:
(1)
P(x1 > x2) = 0, so φ 1 φ 2 ;
(2)
P(x1 > x3) = 1, so φ 3 φ 1 .
Finally, we get the green supplier assessing results that are φ 3 φ 1 φ 2 and the best bearing cage supplier is φ 2 .

6.2. Discussion

In this paper, the index system, which contains four first-level indexes and fifteen second-level indexes, is established. The indexes in the index system are interrelated and sometimes provide feedback-like effects. Since the suppliers to be assessed are independent, we calculate the index weight by rough ANP. Then we process the uncertain index value by membership degree and get the utility value of a supplier on each index. At last, we solve the information incomplete problem in green supplier assessing by evidence theory.
From the case study in Section 6, the green supplier assessment result is φ 3 φ 1 φ 2 . Based on the Overall view of Table 1, it is also known that φ 2 is the best, φ 3 is the worst and φ 1 is middle. The details are as follows:
  • To φ 2 , the performance is the best on ten of the fifteen indexes (i.e., quality (C1,2), service (C1,3), flexibility (C1,4), manufacturing capacity (C2,2), agility capacity (C2,3), geographical environment (C3,2), legal environment (C3,4), technical compatibility degree (C4,1), information platform compatibility degree (C4,3) and reputation (C4,4)).
  • To φ 3 , the performance is the worst on twelve of the fifteen indexes (i.e., cost (C1,1), quality (C1,2), service (C1,3), manufacturing capacity (C2,2), agility capacity (C2,3), economic environment (C3,1), geographical environment (C3,2), social environment (C3,3), legal environment (C3,4), technical compatibility degree (C4,1), cultural compatibility degree (C4,2) and reputation (C4,4)).
  • To φ 1 , the performance is the worst on four indexes (i.e., flexibility (C1,4), innovation capacity (C2,1), geographical environment (C3,2) and information platform compatibility degree (C4,3)), and the performance is the best on five indexes (i.e., cost (C1,1), quality (C1,2), economic environment (C3,1), social environment (C3,3), cultural compatibility degree (C4,2) and reputation (C4,4)).
Thus, it can be seen that the assessing results are in accordance with the actual situation.
According to the utility value in Table 3, we compare the results of the proposed method with those based on Fuzzy Synthetic Evaluation (FSE) [33,34,35] and Fuzzy Analytic Hierarchy Process (FAHP) [36,37,38] as shown in Table 4.
As shown in Table 4, the results of the three methods, which all shows φ 2 is the best supplier, are generally consistent. However, there are differences in the ranking between φ 1 and φ 3 . By FSE method and FAHP method, the complex and interrelated index system is simplified and the index interrelationship information is partially lost. By the proposed method, the weights of the indexes are processed by rough ANP and the index interrelationship information is successfully translated into the hyper-matrix through comparison among indexes.
In addition, green supplier assessment based on rough ANP and evidence theory can provide group decision-making information for enterprise managers, and the index weight can accurately reflect which index has the greatest impact on the supplier selection, thus providing the decision-making basis for the enterprise to reduce the cost and improve the competitiveness. The above analysis and the comparison in Table 4 verify the feasibility and effectiveness of the proposed method for green supplier assessing.

7. Conclusions

In the context of increasingly serious global environmental problems, an ideal manufacturer requires efficient, green suppliers. To handle the uncertainty and incompleteness in green supplier assessment, we propose a green supplier assessment method based on rough ANP and evidence theory. To the best of our knowledge, this is the first attempt to deal with green supplier assessment in a manufacturing enterprise using the hybrid method of rough ANP and evidence theory. The most prominent advantage of the proposed method is that it overcomes the shortcomings of traditional AHP by considering the dependencies and uncertainty across the indexes and processing experts’ judgment on the relative importance of the indexes by fuzzy number and rough boundary interval. It can provide a simple and effective way for weight calculating. By comparing our method with FSE and FAHP approach, we have shown that the proposed method provides a systematic and optimal tool of decision-making for green supplier assessment.
The proposed method provides a way for simplified modeling of complex Multi-criteria Decision Making (MCDM) problems. The decision-making systems are becoming more and more complex nowadays, filled with imprecise and vague information. Evidence theory is adept in capturing such kind of uncertain information, and it provides us with a flexible and effective tool to deal with the green supplier selection problem under uncertain environment. Although the model has been verified on a small case including three potential suppliers and fifteen indexes, it is capable for solving more similar complex problems.
The method proposed in this paper could help us to reduce the risks of making poor investment decisions when dealing with complex networks of green suppliers. In future studies, we will demonstrate on the application of large-scale data sets and the consideration of experts’ reliability.

Author Contributions

Conceptualization, L.L. and H.W.; Methodology, L.L.; Validation, L.L.; Formal Analysis, L.L.; Investigation, L.L.; Writing-Original Draft Preparation, L.L. and H.W.

Funding

This research was funded by Ningxia Natural Science Fund, Grant No. NZ17113 and Ningxia first-class discipline and scientific research projects (electronic science and technology), Grant No. NXYLXK2017A07.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jiang, W.; Huang, C. A Multi-criteria Decision-making Model for Evaluating Suppliers in Green SCM. Int. J. Comput. Commun. Control 2018, 13, 337–352. [Google Scholar] [CrossRef]
  2. Mangla, S.K.; Kumar, P.; Barua, M.K. Flexible Decision Modeling for Evaluating the Risks in Green Supply Chain Using Fuzzy AHP and IRP Methodologies. Glob. J. Flex. Syst. Manag. 2015, 16, 19–35. [Google Scholar] [CrossRef]
  3. Banaeian, N.; Mobli, H.; Fahimnia, B.; Nielsen, I.E.; Omid, M. Green Supplier Selection Using Fuzzy Group Decision Making Methods: A Case Study from the Agri-Food Industry. Comput. Oper. Res. 2016, 89, 337–347. [Google Scholar] [CrossRef]
  4. Yu, F.; Yang, Y.; Chang, D. Carbon footprint based green supplier selection under dynamic environment. J. Clean. Prod. 2018, 170, 880–889. [Google Scholar] [CrossRef]
  5. Gandhi, S.; Mangla, S.K.; Kumar, P.; Kumar, D. A combined approach using AHP and DEMATEL for evaluating success factors in implementation of green supply chain management in Indian manufacturing industries. Int. J. Logist. Res. Appl. 2016, 19, 537–561. [Google Scholar] [CrossRef]
  6. Mangla, S.K.; Kumar, P.; Barua, M.K. Risk analysis in green supply chain using fuzzy AHP approach: A case study. Resour. Conserv. Recycl. 2015, 104, 375–390. [Google Scholar] [CrossRef]
  7. Mangla, S.K.; Kumar, P.; Barua, M.K. Prioritizing the responses to manage risks in green supply chain: An Indian plastic manufacturer perspective. Sustain. Prod. Consum. 2015, 1, 67–86. [Google Scholar] [CrossRef]
  8. Tang, X.; Wei, G. Models for Green Supplier Selection in Green Supply Chain Management with Pythagorean 2-Tuple Linguistic Information. IEEE Access 2018, 6, 18042–18060. [Google Scholar] [CrossRef]
  9. Mangla, S.K.; Govindan, K.; Luthra, S. Critical success factors for reverse logistics in Indian industries: A structural model. J. Clean. Prod. 2016, 129, 608–621. [Google Scholar] [CrossRef]
  10. Saaty, T.L. The Analytic Hierarchy Process; McGraw-Hill: New York, NY, USA, 1980. [Google Scholar]
  11. Saaty, T.L. Fundamentals of the analytic network process—Dependence and feedback in decision-making with a single network. J. Syst. Sci. Syst. Eng. 2004, 13, 129–157. [Google Scholar] [CrossRef]
  12. Noci, G. Designing ‘green’ vendor rating systems for the assessment of a supplier’s environmental performance. Eur. J. Purch. Supply Manag. 1997, 3, 103–114. [Google Scholar] [CrossRef]
  13. Lee, A.H.I.; Kang, H.Y.; Hsu, C.F.; Hung, H.C. A green supplier selection model for high-tech industry. Expert Syst. Appl. 2009, 36, 7917–7927. [Google Scholar] [CrossRef]
  14. Hsu, C.W.; Hu, A.H. Applying hazardous substance management to supplier selection using analytic network process. J. Clean. Prod. 2009, 17, 255–264. [Google Scholar] [CrossRef]
  15. Yeh, W.C.; Chuang, M.C. Using multi-objective genetic algorithm for partner selection in green supply chain problems. Expert Syst. Appl. 2011, 38, 4244–4253. [Google Scholar] [CrossRef]
  16. Yousefi, S.; Shabanpour, H.; Fisher, R.; Saen, R.F. Evaluating and ranking sustainable suppliers by robust dynamic data envelopment analysis. Measurement 2016, 83, 72–85. [Google Scholar] [CrossRef]
  17. Awasthi, A.; Chauhan, S.S.; Goyal, S.K. A fuzzy multicriteria approach for evaluating environmental performance of suppliers. Int. J. Prod. Econ. 2010, 126, 370–378. [Google Scholar] [CrossRef]
  18. Kannan, D.; Jabbour, C.J.C. Selecting green suppliers based on GSCM practices: Using fuzzy TOPSIS applied to a Brazilian electronics company. Eur. J. Oper. Res. 2014, 233, 432–447. [Google Scholar] [CrossRef]
  19. Chatterjee, K.; Pamucar, D.; Zavadskas, E.K. Evaluating the performance of suppliers based on using the R’AMATEL-MAIRCA method for green supply chain implementation in electronics industry. J. Clean. Prod. 2018, 184, 101–129. [Google Scholar] [CrossRef]
  20. Wu, J.; Cao, Q.W.; Li, H. A Method for Choosing Green Supplier Based on COWA Operator under Fuzzy Linguistic Decision-Making. J. Ind. Eng. Eng. Manag. 2010, 24, 61–65. [Google Scholar]
  21. Luo, X.X.; Peng, S.H. Research on the Vendor Evaluation and Selection Based on AHP and TOPSIS in Green Supply Chain. Soft Sci. 2011, 25, 53–56. [Google Scholar]
  22. Kuo, R.J.; Lin, Y.J. Supplier selection using analytic network process and data envelopment analysis. Int. J. Prod. Res. 2012, 50, 2852–2863. [Google Scholar] [CrossRef]
  23. Shi, L. Green Supplier Evaluation of RS-RBF Neural Network Model. Sci. Technol. Manag. Res. 2012, 32, 198–201. [Google Scholar]
  24. Akman, G. Evaluating suppliers to include green supplier development programs via fuzzy c-means and VIKOR methods. Comput. Ind. Eng. 2015, 86, 69–82. [Google Scholar] [CrossRef]
  25. Pamučar, D.; Petrović, I.; Ćirović, G. Modification of the Best-Worst and MABAC methods: A novel approach based on interval-valued fuzzy-rough numbers. Expert Syst. Appl. 2018, 91, 89–106. [Google Scholar] [CrossRef]
  26. Li, L.; Mo, R.; Chang, Z.; Zhang, H. Priority evaluation method for aero-engine assembly task based on balanced weight and improved TOPSIS. Comput. Integr. Manuf. Syst. 2015, 21, 1193–1201. (In Chinese) [Google Scholar]
  27. Wang, X.; Xiong, W. Rough AHP approach for determining the importance ratings of customer requirements in QFD. Comput. Integr. Manuf. Syst. 2010, 16, 763–771. (In Chinese) [Google Scholar]
  28. Vasiljević, M.; Fazlollahtabar, H.; Stević, Z.; Vesković, S. A rough multicriteria approach for evaluation of the supplier criteria in automotive industry. Decis. Mak. Appl. Manag. Eng. 2018, 1, 82–96. [Google Scholar] [CrossRef]
  29. Pamučar, D.; Stević, Ž.; Zavadskas, E.K. Integration of interval rough AHP and interval rough MABAC methods for evaluating university web pages. Appl. Soft Comput. 2018, 67, 141–163. [Google Scholar] [CrossRef]
  30. Zhu, G.N.; Hu, J.; Qi, J.; Gu, C.C.; Peng, Y.H. An integrated AHP and VIKOR for design concept evaluation based on rough number. Adv. Eng. Inform. 2015, 29, 408–418. [Google Scholar] [CrossRef]
  31. Shafer, G. A Mathematical Theory of Evidence; Princeton University Press: Princeton, NJ, USA, 1976. [Google Scholar]
  32. Sarabi-Jamab, A.; Araabi, B.N. How to decide when the sources of evidence are unreliable: A multi-criteria discounting approach in the Dempster–Shafer theory. Inf. Sci. 2018, 448, 233–248. [Google Scholar] [CrossRef]
  33. Li, R.; Jin, Y. The early-warning system based on hybrid optimization algorithm and fuzzy synthetic evaluation model. Inf. Sci. 2017, 435, 296–319. [Google Scholar] [CrossRef]
  34. Haider, H.; Hewage, K.; Umer, A.; Ruparathna, R.; Chhipi-Shrestha, G.; Culver, K.; Holland, M.; Kay, J.; Sadiq, R. Sustainability Assessment Framework for Small-sized Urban Neighbourhoods: An Application of Fuzzy Synthetic Evaluation. Sustain. Cities Soc. 2017, 36, 21–32. [Google Scholar] [CrossRef]
  35. Zhu, J. Evaluation of supplier strength based on fuzzy synthetic assessment method. In Proceedings of the International Conference on Test and Measurement, Hong Kong, China, 5–6 December 2009; pp. 247–250. [Google Scholar]
  36. Deepika, M.; Kannan, A.S.K. Global supplier selection using intuitionistic fuzzy Analytic Hierarchy Process. In Proceedings of the International Conference on Electrical, Electronics, and Optimization Techniques, Chennai, India, 3–5 March 2016; pp. 2390–2395. [Google Scholar]
  37. Torng, C.; Tseng, K.W. Using Fuzzy Analytic Hierarchy Process to Construct Green Suppliers Assessment Criteria and Inspection Exemption Guidelines. In Proceedings of the International Conference on Human-Computer Interaction; Springer: Berlin/Heidelberg, Germany, 2013; pp. 729–732. [Google Scholar]
  38. Labib, A.W. A supplier selection model: A comparison of fuzzy logic and the analytic hierarchy process. Int. J. Prod. Res. 2011, 49, 6287–6299. [Google Scholar] [CrossRef]
Figure 1. The index system.
Figure 1. The index system.
Information 09 00162 g001
Figure 2. The ANP structure of green supplier assessment.
Figure 2. The ANP structure of green supplier assessment.
Information 09 00162 g002
Figure 3. The green supplier assessment procedure.
Figure 3. The green supplier assessment procedure.
Information 09 00162 g003
Table 1. The index values of three bearing cage suppliers.
Table 1. The index values of three bearing cage suppliers.
SupplierC1,1C1,2C1,3C1,4C2,1C2,2C2,3C3,1C3,2C3,3C3,4C4,1C4,2C4,3C4,4
φ 1 6.4 × 1010.010.70.1[0.1, 0.3]0.50.70.50.50.70.7/[0.7, 0.9]0.50.7
φ 2 1.9 × 1030.010.90.70.50.70.9/0.7[0.5, 0.7]0.90.70.50.70.9
φ 3 2.8 × 1040.030.10.70.70.10.30.10.50.10.30.10.5/0.1
Table 2. The normalized index values of three bearing cage suppliers.
Table 2. The normalized index values of three bearing cage suppliers.
SupplierC1,1C1,2C1,3C1,4C2,1C2,2C2,3C3,1C3,2C3,3C3,4C4,1C4,2C4,3C4,4
φ 1 0.00231.00000.77780.1429[0.1429, 0.4286]0.71430.77781.00000.71431.00000.7778/[0.7778, 1.0000]0.71430.7778
φ 2 0.06691.00001.00001.00000.71431.00001.0000/1.0000[0.7143, 1.0000]1.00001.00000.55561.00001.0000
φ 3 1.00000.33330.11111.00001.00000.14290.33330.20000.71430.14290.33330.14290.5556/0.1111
Table 3. The utility values of three bearing cage suppliers.
Table 3. The utility values of three bearing cage suppliers.
SupplierC1,1C1,2C1,3C1,4C2,1C2,2C2,3C3,1C3,2C3,3C3,4C4,1C4,2C4,3C4,4
φ 1 0.10000.90000.77780.10000.28570.63330.63340.90000.10000.90000.6334/0.88890.10000.7000
φ 2 0.15180.90000.90000.90000.63330.90000.9000/0.90000.85720.90000.90000.10000.90000.9000
φ 3 0.90000.10000.10000.90000.90000.10000.10000.10000.10000.10000.10000.10000.1000/0.1000
Table 4. The comparison of assessing results of the proposed method, FSE and FAHP.
Table 4. The comparison of assessing results of the proposed method, FSE and FAHP.
SupplierRanking of Three Suppliers
Proposed MethodFSEFAHP
φ 1 233
φ 2 111
φ 3 322

Share and Cite

MDPI and ACS Style

Li, L.; Wang, H. A Green Supplier Assessment Method for Manufacturing Enterprises Based on Rough ANP and Evidence Theory. Information 2018, 9, 162. https://doi.org/10.3390/info9070162

AMA Style

Li L, Wang H. A Green Supplier Assessment Method for Manufacturing Enterprises Based on Rough ANP and Evidence Theory. Information. 2018; 9(7):162. https://doi.org/10.3390/info9070162

Chicago/Turabian Style

Li, Lianhui, and Hongguang Wang. 2018. "A Green Supplier Assessment Method for Manufacturing Enterprises Based on Rough ANP and Evidence Theory" Information 9, no. 7: 162. https://doi.org/10.3390/info9070162

APA Style

Li, L., & Wang, H. (2018). A Green Supplier Assessment Method for Manufacturing Enterprises Based on Rough ANP and Evidence Theory. Information, 9(7), 162. https://doi.org/10.3390/info9070162

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop