1 Introduction
Nowadays, organizations need to make decisions for different matters. Employing a suitable approach to make a correct decision is an ongoing concern of organizations. There are many types of methods to solve multiple criteria decision-making (MCDM) problems, which can be categorized into pairwise comparison-based methods (Doumpos and Zopounidis, 2004), distance-based methods, outranking methods (Liao and Wu, 2020) and utility function-based methods. Analytic hierarchy process (AHP) is the most well-known pairwise comparison-based method that has been utilized in a number of researches in recent decades (Saaty, 1979). In 2015, the best-worst method (BWM) was proposed by Rezaei (2015), which needs less comparisons on criteria and/or alternatives compared with the AHP (Rezaei, 2016). This method has received the attention of scholars in recent years because of its reasonable performance (Mi et al., 2019).
In real-world situations, the input data of decision-making problems include uncertainty and/or incompleteness (Zavadskas et al., 2010). In this regard, the methods which account for the lack of certainty in decision making have gained ever-increasing importance (Chen, 2018; Lin and Chen, 2019). Grey system theory has some features such as reality and vague aspects and can be used for small samples. This theory does not require distribution and membership function. An interval grey number is a number that belongs to the interval $[a,b]$ with $a<b$, which is not similar to the interval number within $[a,b]$. A grey number may appear in three cases of continuous, discrete, or random grey number, depending on the nature of its uncertainty (Liu et al., 2012). To compare grey theory with statistics, probability, interval and fuzzy theory, we can specify the type of their uncertainty (Ng and Deng, 2005).
The main objective of the current study is to present the grey BWM (GBWM). Specifically, the contributions of this study can be summarized as follows:
-
• We introduce the GBWM and the result calculated by this method is more reliable than the fuzzy BWM because the GBWM has a smaller inconsistency ratio compared with the fuzzy BWM.
-
• The GBWM uses a linear model that can present the global-optimum weights for MCDM problems, while the existing fuzzy BWM employs a non-linear model with local-optimum weights.
The current study is organized as follows: in Section 2, we provide the detailed knowledge regarding the BWM. In Section 3, the grey system theory and its components are examined, such as the grey linear programming, the basic concepts of grey numbers, and then we propose the GBWM with the linear-type model. In Section 4, a practical example with sensitivity analysis is conducted and the GBWM is compared with the fuzzy BWM in terms of performance. The paper ends with some conclusions and future directions.
2 Best Worst Method
AHP is one of the most well-known methods for MCDM. However, the BWM deduces more consistent weights based on the less comparisons compared with the AHP (Rezaei, 2015). In this method, the necessary criteria for decision-making are determined firstly. Then, the best and the worst criteria are specified. The next step is to compare the other criteria against the best and the worst criteria. A min–max mathematical model is formulated and solved. The ultimate output of this model would determine the weight of each criterion. Subsequently, criteria can be ranked according to their weights. The steps of the original BWM can be described briefly as follows Rezaei (2015):
Step 1: Determine a set of criteria for decision-making (this is done by the decision-maker): $\{{c_{1}},{c_{2}},\dots ,{c_{n}}\}$.
Step 2: Determine the best $(B)$ and the worst $(W)$ criteria based on their importance degrees.
Step 3: Determine the degree of preference of the best criterion to the other criteria, which is done using a number from 1 to 9. Then, the best-to-others (BO) vector is displayed as ${A_{B}}=({a_{B1}},{a_{B2}},{a_{B3}},\dots ,{a_{Bn}})$, where ${a_{Bn}}$ demonstrates the preference of criterion ${c_{B}}$ to criterion ${c_{n}}$.
Step 4: Determine the level of preference of each criterion to the worst criterion, which is done using a number from 1 to 9. Then, the others-to-worst (OW) vector is displayed as ${A_{W}}={({a_{1W}},{a_{2W}},{a_{3W}},\dots ,{a_{nW}})^{T}}$, where ${a_{jW}}$ indicates the preference of criterion ${c_{j}}$ to criterion ${c_{W}}$.
Step 5: The objective of this step is to determine the weights of criteria by minimizing the maximum absolute differences $\big|\frac{{W_{B}}}{{W_{j}}}-{a_{Bj}}\big|$ and $\big|\frac{{W_{j}}}{{W_{w}}}-{a_{jw}}\big|$ for all j. To solve this problem, a min–max Model (1) is formed:
(1)
\[\begin{aligned}{}& \operatorname{Min}\underset{j}{\max }\bigg\{\bigg|\frac{{W_{B}}}{{W_{j}}}-{a_{Bj}}\bigg|,\bigg|\frac{{W_{j}}}{{W_{w}}}-{a_{jw}}\bigg|\bigg\}\\ {} & \text{s.t.}\\ {} & \sum \limits_{j}{W_{j}}=1,\\ {} & {W_{j}}\geqslant 0,\hspace{1em}\text{for all}\hspace{2.5pt}j.\end{aligned}\]Model (1) is equivalent to Model (2) by using ξ to denote the maximal deviation.
(2)
\[\begin{aligned}{}& \min \xi \\ {} & \text{s.t.}\\ {} & \bigg|\frac{{W_{B}}}{{W_{j}}}-{a_{Bj}}\bigg|\leqslant \xi ,\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & \bigg|\frac{{W_{j}}}{{W_{w}}}-{a_{jw}}\bigg|\leqslant \xi ,\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & \sum \limits_{j}{W_{j}}=1{W_{j}}\geqslant 0,\hspace{1em}\text{for all}\hspace{2.5pt}j.\end{aligned}\]Solving Model (2), the value of the optimal weight is obtained for each criterion. Based on the obtained weights, the criteria would be prioritized. The criterion with the highest weight would have the highest priority.
In this method, the consistency ratio, calculated by Eq. (1), is a value ranging from 0 to 1. A value closer to 0 indicates higher consistency.
where the Consistency Index is determined based on Table 1 proportional to the number of criteria.
(1)
\[ \mathit{Consistency}\hspace{2.5pt}\mathit{Ratio}=\frac{\xi }{\mathit{Consistency}\hspace{2.5pt}\mathit{Index}},\]Table 1
Consistency Index (Rezaei, 2015).
| The maximal preference degree of the best over the worst | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 |
| Consistency Index | 0.00 | 0.44 | 1.00 | 1.63 | 2.30 | 3.00 | 3.73 | 4.47 | 5.23 |
Considering the advantages, the BWM has gained ever-increasing investigation in recent years (Mi et al., 2019).
2.1 Applications of the BWM
In this section, we are going to examine applications of BWM that have been conducted in different fields. Rezaei et al. (2015) presented an integrated approach with the BWM, which includes willingness and capabilities as two features for classifying and evaluating the suppliers. This integrated approach helps companies divide their managerial resources effectively. In another research, Rezaei et al. (2016) discussed that choosing suppliers is a strategic decision that significantly influences the competitive advantage of a company. Due to the importance of selecting the suppliers and the sensitivity of this issue, they presented a novel method which aimed to select the suppliers during three phases, namely pre-selection, selection, and aggregation. They used BWM for initial screening. As for the decision-making phase, they made use of an MCDM method and the final decision was made at the aggregation phase based on the material price and the required annual quantity.
Scholars have applied the original BWM into different fields. Gupta and Barua (2016) conducted a study to identify the factors affecting technological innovations in India. Making use of the BWM, they presented a procedure for selecting the best empowerers. The primary purpose of their study was to identify vital enablers in the field of technological innovation in India. Rezaei et al. (2017) considered three alternatives related to transportation: through unit load devices, mixed unit load devices, and loose freight in trucks. To find the optimal freight bundling configuration, they considered three key performance indicators including cost, quality, and loading time. They used the BWM for finding the most suitable type of bundling configuration in the surface transportation. Their proposed framework facilitated risk assessment process. Torabi et al. (2016) used the BWM to present an enhanced risk assessment method within the business continuity management system framework. Sadaghiani et al. (2017) used the BWM to assess the importance of different types of energy such as oil and gas industries by sending questionnaires to academic experts on sustainable supply chain management. Their research helped companies to develop strategies for identifying external forces. The treatment of urban sewage sludge for reducing the threats of environmental pollution and its negative influences on human health is of high significance. Accordingly, Ren et al. (2017) presented a general framework for selecting the proper technology for the treatment of urban sewage. In this framework, the BWM was used to determine the weights of the criteria. Salimi and Rezaei (2018) measured the R&D performance using the BWM. The SERVQUAL model is designed to evaluate the service quality of a baggage handling system. Rezaei et al. (2018) collected the list of criteria for the SERVQUAL model and utilized the BWM to determine the weights of the criteria.
2.2 BWM Extensions in Uncertainty Theories
For some situations, the uncertain information of experts is easy to access and the uncertain BWMs have attracted the attention of scholars (Mi et al., 2019). Making some changes to the steps of the BWM, Rezaei (2016) presented a linear model, by which the weight of each criterion was obtained in intervals and the criteria were prioritized after comparing the obtained intervals. Using the fuzzy approach and triangular fuzzy numbers, Guo and Zhao (2017) presented the fuzzy BWM (FBWM). They believed that the FBWM obtained the suitable ranking of alternatives and that it has more comparison consistency than the BWM. In addition, Hafezalkotob and Hafezalkotob (2017) presented a fuzzy approach which combined individual and group alternatives, and used the FBWM in their approach in which triangular fuzzy numbers were employed. Gupta et al. (2017) presented an accurate and comprehensive approach for identifying various factors of managing energy efficiency in India. They believed that identifying the barriers is not enough and other proper solutions would also be presented. To rank the existing barriers, they used the BWM in their research. Mou et al. (2016) used intuitionistic fuzzy multiplicative preference relations for ranking criteria or alternatives and it can be considered as a tool for combining with other MCDM methods. It should be noted that their method did not consider the case that the consistency degree is not acceptable. Aboutorab et al. (2018) proposed the ZBWM by utilizing the Z-number approach to consider the uncertainty of the input data in MCDM problems. Their method has lower inconsistency ratio compared with the BWM.
Mou et al. (2017) presented a graph-based group decision-making approach for intuitionistic fuzzy BWM. Moreover, they presented three numerical examples to show the efficiency of the proposed method. Mi and Liao (2019) enabled BWM to accept hesitant numbers as input and Liao et al. (2019) fused the hesitant linguistic information in BWM. Mi et al. (2019) conducted a survey of the BWM related publications between 2015 and 2019. Hafezalkotob et al. (2019) proposed the interval MULTIMOORA and group interval BWM.
Since there is no research item on the combination of the grey system theory and BWM, the current study attempts to present the GBWM. As mentioned, the grey theory has features such as reality and vague aspects, and can be employed with an incomplete data, therefore, it can be a useful approach to solve decision-making problems (Zavadskas et al., 2009). The result calculated by the GBWM is more reliable than the one calculated by the fuzzy BWM because the GBWM has a smaller inconsistency ratio compared with the fuzzy BWM. We transform the weight determination model in the GBWM into linear models which can be solved efficiently and own high reliability. It should be noted that the rankings derived by the FBWM and GBWM are the same in the comparative analysis. Moreover, after defuzzification, the weights obtained by the FBWM fall within the weights that are obtained from the GBWM. Hence, the proposed method has appropriate performance in calculating weights. The calculation of grey numbers is simple compared to other uncertainty approaches. In the next section, the basic concepts of grey system theory are examined.
3 Grey System Theory
The grey system theory, proposed in 1980s by Julong (1989), is a mathematical concept which has a widespread application in MCDM. It is considered as a highly effective method for encountering uncertainty problems associated with unknown and incomplete information (Liu and Lin, 2010). Generally, the information pertaining to the preferences of decision-makers for certain criteria and the reasons for such preferences are expressed based on the qualitative judgments of decision-makers. Also, in practice, the judgments of decision-makers are often uncertain and thus cannot be represented by accurate numerical values. The grey theory is one of the concepts used for studying uncertainty and incompleteness. This theory has been used in the mathematical analysis of incomplete information system (Chalekaee et al., 2019; Mahmoudi and Feylizadeh, 2018). The importance degrees of criteria in a decision-making process can be expressed by numerical intervals. These numerical intervals would include uncertain information. In other words, the accurate values of grey numbers are unknown, but the interval which covers a value is almost known (Liu et al., 2017). Since we will compare the GBWM with the FBWM in the current study, the grey system theory and fuzzy set theory (Zadeh, 1965) are compared with each other from different aspects in Table 2 (Mahmoudi et al., 2019).
Table 2
The comparison between grey system theory and fuzzy set theory.
| Uncertainty research | Grey system | Fuzzy math |
| Research object | Poor information | Cognitive |
| Basic set | Grey number set | Fuzzy set |
| Describe method | Possibility function | Membership function |
| Procedure | Sequence operator | Cut set |
| Data requirement | Any distribution | Known membership |
| Emphasis | Intension | Extension |
| Objective | Law of reality | Cognitive expression |
| Characteristic | Small data | Depend on experience |
3.1 Preliminaries of Grey Numbers
In this section, the related preliminaries of grey numbers are reviewed.
A grey number (Liu et al., 2012) is expressed as
To denote the central point of grey numbers, the “Kernel” of grey numbers (Guo et al., 2017) is proposed as:
For two grey numbers $\otimes A$ and $\otimes B$, the arithmetic operations can be defined (Mahmoudi and Feylizadeh, 2018; Turskis and Zavadskas, 2010) as:
(3)
\[\begin{aligned}{}& \otimes A+\otimes B=[\underline{A}+\underline{B},\overline{A}+\overline{B}],\end{aligned}\](4)
\[\begin{aligned}{}& \otimes A-\otimes B=[\underline{A}-\overline{B},\overline{A}-\underline{B}],\end{aligned}\](5)
\[\begin{aligned}{}& \otimes A.\otimes B=\big[\operatorname{Min}\{\underline{A}\underline{B},\overline{A}\overline{B},\overline{A}\underline{B},\underline{A}\overline{B}\},\operatorname{Max}\{\underline{A}\underline{B},\overline{A}\overline{B},\overline{A}\underline{B},\underline{A}\overline{B}\}\big],\end{aligned}\](6)
\[\begin{aligned}{}& \frac{\otimes A}{\otimes B}=\otimes A.\otimes {B^{-1}}=\bigg[\operatorname{Min}\bigg\{\frac{\underline{A}}{\underline{B}},\frac{\underline{A}}{\overline{B}},\frac{\overline{A}}{\underline{B}},\frac{\overline{A}}{\overline{B}}\bigg\},\operatorname{Max}\bigg\{\frac{\underline{A}}{\underline{B}},\frac{\underline{A}}{\overline{B}},\frac{\overline{A}}{\underline{B}},\frac{\overline{A}}{\overline{B}}\bigg\}\bigg].\end{aligned}\]To compare the grey numbers, the greyness degree of $\otimes A$ was calculated as [32]:
where Ω represents the background of grey numbers and μ is the length of the background of grey numbers. For two grey numbers $\otimes A=[\underline{A},\overline{A}]$ and $\otimes B=[\underline{B},\overline{B}]$, if $\otimes \hat{A}<\otimes \hat{B}$, then $\otimes A{<_{G}}\otimes B;$ if $\otimes \hat{A}=\otimes \hat{B}$, then if ${g^{o}}(\otimes A)={g^{o}}(\otimes B)$, then $\otimes A{=_{G}}\otimes B$; if ${g^{o}}(\otimes A)<{g^{o}}(\otimes B)$, then $\otimes A{>_{G}}\otimes B$.
3.2 Grey Linear Programming (GLP)
Different methods have been presented for solving GLP models. Huang et al. (1995) presented a method for solving grey mixed-integer linear programming. Their model was suitable for grey models with the same sign in lower and upper bounds of grey numbers. After that, Li (2007) proposed another method to solve GLP problems named “Covered Solution”. The disadvantages of the method were complex calculations and it sometimes fails to meet stop conditions. Hajiagha et al. (2012) proposed a method to solve the GLP problem by using a multi-objective concept, yet their method presented the wrong solution for GLP problems as proved by Mahmoudi et al. (2018a). Li et al. (2014) proposed a method based on the concept of Covered Solution method, yet it had some problems, similarly as Li (2007)’s method. Nasseri et al. (2016) presented a new method using the primal simplex algorithm to solve GLP problems, but their method could solve the GLP problems just with the grey objective function. Liu et al. (2009) presented a positioned programming for solving GLP models. This method truly enjoys simplicity and covers all uncertainties in grey numbers. Moreover, this method can present crisp values based on ρ, β and δ parameters that are determined by the decision maker.
The current study employs the positioned programming method to solve GLP problems in the form of Model (3) (Liu et al., 2009; Mahmoudi et al., 2018b).
Such that:
It should be noted that $A(\otimes )$ is the grey consumption matrix, $C(\otimes )$ is the grey price vector, $b(\otimes )$ is the grey constraint vector, and X is the problem decision vector. The parameters employed in Eqs. (10) to (12) are defined in Eqs. (13) to (15):
(3)
\[\begin{aligned}{}& \operatorname{Max}S=C(\otimes )X\\ {} & \text{s.t.}\\ {} & A(\otimes \hspace{2.5pt})X\leqslant b(\otimes ),\\ {} & X\geqslant 0.\end{aligned}\](10)
\[\begin{aligned}{}& C(\otimes )={\big[{C_{1}}(\otimes ),{C_{2}}(\otimes ),\dots ,{C_{n}}(\otimes )\big]^{T}},\end{aligned}\](11)
\[\begin{aligned}{}& b(\otimes )={\big[{b_{1}}(\otimes ),{b_{2}}(\otimes ),\dots ,{b_{m}}(\otimes )\big]^{T}},\end{aligned}\](12)
\[\begin{aligned}{}& A(\otimes )=\left[\begin{array}{c@{\hskip4.0pt}c@{\hskip4.0pt}c@{\hskip4.0pt}c}{a_{11}}(\otimes )\hspace{1em}& {a_{11}}(\otimes )\hspace{1em}& \dots \hspace{1em}& {a_{1n}}(\otimes )\\ {} {a_{21}}(\otimes )\hspace{1em}& {a_{11}}(\otimes )\hspace{1em}& \dots \hspace{1em}& {a_{2n}}(\otimes )\\ {} \dots \hspace{1em}& \dots \hspace{1em}& \dots \hspace{1em}& \dots \\ {} {a_{m1}}(\otimes )\hspace{1em}& {a_{11}}(\otimes )\hspace{1em}& \dots \hspace{1em}& {a_{mn}}(\otimes )\end{array}\right].\end{aligned}\](13)
\[\begin{aligned}{}& {c_{j}}(\otimes )\in [{\underline{c}_{j}},{\overline{c}_{j}}],\hspace{1em}{\underline{c}_{j}}\geqslant 0,\hspace{2.5pt}j=1,2,\dots ,n,\end{aligned}\]To solve the GLP model, it ought to be whitened first.
Definition 1.
If the values of ${\delta _{ij}}$, ${\beta _{j}}$ and ${\rho _{j}}$ for $i=1\dots m$ and $j=1\dots n$ fall within the closed interval $[0,1]$, the white values of the grey parameters are defined by Eqs. (16) to (18).
(16)
\[\begin{aligned}{}& {\tilde{c}_{j}}(\otimes )={\rho _{j}}{\overline{c}_{j}}+(1-{\rho _{j}}){\underline{c}_{j}};\hspace{1em}j=1,2,\dots ,n,\end{aligned}\]After the whitening stage, we are left with Model (4):
If $(\rho ,\beta ,\delta )=(0,0,1)$, we have the lowest value after solving model (4). It can be displayed by $\operatorname{Max}\underline{s}$. On the other hand, if $(\rho ,\beta ,\delta )=(1,1,0)$, we have the highest value after solving model (4), represented by $\operatorname{Max}\overline{s}$.
This theorem has been proven in Liu et al. (2009).
3.3 Grey Best Worst Method (GBWM)
In this section, the GBWM is presented. The eight steps of this method are outlined in detail as follows:
Step 1: Determine the criteria set $\{{c_{1}},{c_{2}},\hspace{2.5pt}{c_{3}},\dots ,{c_{n}}\}$ by decision-makers.
Step 2: Each decision-maker determines the best and the worst criteria. If there are k experts, then k best criteria and the k worst criteria would exist: $\{{B^{p1}},{B^{p2}},{B^{p3}},\dots ,{B^{pk}}\}$ and $\{{W^{p1}},{W^{p2}},{W^{p3}},\dots ,{W^{pk}}\}$.
Step 3: At this step, each decision-maker determines the degrees of preference of the best criterion to the other criteria using the linguistic variables presented in Table 3. Equation (20) expresses best-to-others (BO) vectors of experts.
(20)
\[ \begin{aligned}{}& \otimes {A_{B}^{p1}}=\big(\otimes {a_{B1}^{p1}},\otimes {a_{B2}^{p1}},\otimes {a_{B3}^{p1}}\cdots \otimes {a_{Bn}^{p1}}\big),\\ {} & \hspace{99.58464pt}\dots \dots \\ {} & \otimes {A_{B}^{pk}}=\big(\otimes {a_{B1}^{pk}},\otimes {a_{B2}^{pk}},\otimes {a_{B3}^{pk}}\cdots \otimes {a_{Bn}^{pk}}\big).\end{aligned}\]Table 3
Linguistic variables of decision-makers.
| Linguistic variable | Value |
| Equally Important (EI) | $[1,1]$ |
| Weakly Important (WI) | $\big[\frac{2}{3},\frac{3}{2}\big]$ |
| Fairly Important (FI) | $\big[\frac{3}{2},\frac{5}{2}\big]$ |
| Very Important (VI) | $\big[\frac{5}{2},\frac{7}{2}\big]$ |
| Absolutely Important (AI) | $\big[\frac{7}{2},\frac{9}{2}\big]$ |
In Eq. (20), $\otimes {A_{B}^{p1}}$ represents the opinion of the first decision-maker determining the degree of preference of the best criterion to criteria 1 to n.
Step 4: In Eq. (21), $\otimes {A_{W}^{pk}}$ represents the viewpoint of the kth decision-maker determining the degree of preference of criteria 1 to n to the worst criterion. Equation (21) expresses others-to-worst (OW) vectors.
(21)
\[ \begin{aligned}{}& \otimes {A_{W}^{p1}}={\big(\otimes {a_{1W}^{p1}},\otimes {a_{2W}^{p1}},\otimes {a_{3W}^{p1}}\cdots \otimes {a_{nW}^{p1}}\big)^{T}},\\ {} & \hspace{99.58464pt}\dots \dots \\ {} & \otimes {A_{W}^{pk}}={\big(\otimes {a_{1W}^{pk}},\otimes {a_{2W}^{pk}},\otimes {a_{3W}^{pk}}\cdots \otimes {a_{nW}^{pk}}\big)^{T}}.\end{aligned}\]
Step 5: At this stage, the degree of optimal weight for each criterion is determined. Since the inputs of the problem are considered in grey numbers, Model (2) is converted into a grey model. For this purpose, we consider Model (5):
(5)
\[\begin{aligned}{}& \min \xi \\ {} & \text{s.t.}\\ {} & \bigg|\frac{{W_{B}}}{{W_{j}}}-{a_{Bj}}\bigg|\leqslant \xi ,\hspace{1em}\text{for all}j,\\ {} & \bigg|\frac{{W_{j}}}{{W_{w}}}-{a_{jw}}\bigg|\leqslant \xi ,\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & \sum \limits_{j}{W_{j}}=1,\\ {} & {W_{j}}\geqslant 0,\hspace{1em}\text{for all}\hspace{2.5pt}j.\end{aligned}\]According to the features of absolute value, Model (5) is equivalent to Model (6):
(6)
\[\begin{aligned}{}& \min \xi \\ {} & \text{s.t.}\hspace{2.5pt}\\ {} & \frac{{W_{B}}}{{W_{j}}}-{a_{Bj}}\leqslant \xi ,\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & -\frac{{W_{B}}}{{W_{j}}}+{a_{Bj}}\leqslant \xi ,\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & \frac{{W_{j}}}{{W_{w}}}-{a_{jw}}\leqslant \xi ,\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & -\frac{{W_{j}}}{{W_{w}}}+{a_{jw}}\leqslant \xi ,\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & \sum \limits_{j}{W_{j}}=1,\\ {} & {W_{j}}\geqslant 0,\hspace{1em}\text{for all}\hspace{2.5pt}j.\end{aligned}\]To cross multiply the constraints of Model (6), Model (7) is ultimately obtained:
(7)
\[\begin{aligned}{}& \min \xi \\ {} & \text{s.t.}\hspace{2.5pt}\\ {} & {W_{B}}-{W_{j}}{a_{Bj}}\leqslant {W_{j}}\xi ,\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & -{W_{B}}+{W_{j}}{a_{Bj}}\leqslant {W_{j}}\xi ,\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & {W_{j}}-{W_{w}}{a_{jw}}\leqslant {W_{w}}\xi ,\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & -{W_{j}}+{W_{w}}{a_{jw}}\leqslant {W_{w}}\xi ,\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & \sum \limits_{j}{W_{j}}=1,\\ {} & {W_{j}}\geqslant 0,\hspace{1em}\text{for all}\hspace{2.5pt}j.\end{aligned}\]Due to the multiplication of variables ξ and ${W_{j}}$, Model (7) is a non-linear model. Since solving grey non-linear models involves high levels of complexity, the model should be converted into a linear model. Also, since two continuous variables cause the non-linearity in Model (7), the McCormick method (Hijazi et al., 2017; McCormick, 1976) can be used for linearization. The steps are as follows:
In Eq. (22), variables ${x_{1}}$ and ${x_{2}}$ are continuous and have specific upper and lower limits. Also, variable ${\varnothing _{1}}$ has been considered as the product of multiplying the variables ${x_{1}}$ and ${x_{2}}$. By considering Eq. (22) and adding four constraints mentioned in Eqs. (23) to (26), the linearization operation for variables ${x_{1}}$ and ${x_{2}}$ is undertaken.
(23)
\[\begin{aligned}{}& {\varnothing _{1}}\geqslant {x_{1}^{L}}{x_{2}}+{x_{2}^{L}}{x_{1}}-{x_{1}^{L}}{x_{2}^{L}},\end{aligned}\](24)
\[\begin{aligned}{}& {\varnothing _{1}}\geqslant {x_{1}^{U}}{x_{2}}+{x_{2}^{U}}{x_{1}}-{x_{1}^{U}}{x_{2}^{U}},\end{aligned}\]In this section, the nonlinear Model (7) is converted into a linear model. Using Eq. (22), we first have the assumptions mentioned in Eq. (27):
The assumptions mentioned in Eq. (27) are not sufficient for linearization and the range of variable ξ should be determined in this regard. Based on Eq. (1), for variable ξ, we have
The value of $\mathit{CI}$ is determined based on Table 1 and the value of $\mathit{CR}$ is always within the interval of $[0,1]$. In real contexts, the decision-maker may want the inconsistency rate not to be greater than the specified value of A. Therefore, the range of $\mathit{CR}$ is considered as $[0,A]$. Equation (29) is thus formed:
Based on Eqs. (27) and (29) and the constraints from Eq. (23) to Eq. (26), the nonlinear Model (7) is changed into the grey linear Model (8):
(8)
\[\begin{aligned}{}& \min \otimes \xi \\ {} & \text{s.t.}\\ {} & \otimes {W_{B}}-\otimes {W_{j}}\otimes {a_{Bj}}\leqslant \otimes {\varnothing _{1}},\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & -\otimes {W_{B}}+\otimes {W_{j}}\otimes {a_{Bj}}\leqslant \otimes {\varnothing _{1}},\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & \otimes {\varnothing _{1}}\geqslant 0,\\ {} & \otimes {\varnothing _{1}}\geqslant \otimes \xi +\mathit{CI}.A.\otimes {W_{j}}-\mathit{CI}.A,\\ {} & \otimes {\varnothing _{1}}\leqslant \mathit{CI}.A.\otimes {W_{j}},\\ {} & \otimes {\varnothing _{1}}\leqslant \otimes \xi ,\\ {} & \otimes {W_{j}}-\otimes {W_{w}}\otimes {a_{jw}}\leqslant \otimes {\varnothing _{2}},\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & -\otimes {W_{j}}+\otimes {W_{w}}\otimes {a_{jw}}\leqslant \otimes {\varnothing _{2}},\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & \otimes {\varnothing _{2}}\geqslant 0,\\ {} & \otimes {\varnothing _{2}}\geqslant \otimes \xi +\mathit{CI}.A.\otimes ,{W_{w}}-\mathit{CI}.A,\\ {} & \otimes {\varnothing _{2}}\leqslant \mathit{CI}.A.\otimes {W_{w}},\\ {} & \otimes {\varnothing _{2}}\leqslant \otimes \xi ,\\ {} & \overline{{W_{j}}}-\underline{{W_{j}}}\geqslant \varepsilon ,\hspace{1em}\text{for all}\hspace{2.5pt}j,\\ {} & \sum \limits_{j}\otimes {W_{j}}=[0.8,1.2],\\ {} & \otimes {W_{j}}\geqslant 0,\hspace{1em}\text{for all}\hspace{2.5pt}j,\hspace{2.5pt}0\leqslant A\leqslant 1,\hspace{2.5pt}\mathit{CI}\geqslant 0.\end{aligned}\]To find the optimal grey weights, the grey linear Model (8) can be formed and then solved using on positioned programming approach.
Step 6: Based on the opinion of each decision-maker, a weight is assigned to each criterion. To integrate the viewpoint of the experts with regard to each criterion, the grey geometric mean relation as defined in Eq. (30) is used. Since the opinions of decision-makers are of different significances, the weights (${W_{k}}$) for each expert are determined using the linguistic variables presented in Table 4.
Step 7. At this stage, the obtained weights are normalized by Eq. (31) (Dey and Chakraborty, 2016):
(31)
\[ \otimes {W_{j}^{\ast }}=\bigg(\frac{\underline{{W_{j}}}}{\frac{1}{2}[{\textstyle\textstyle\sum _{j=1}^{n}}\underline{{W_{j}}}+{\textstyle\textstyle\sum _{j=1}^{n}}\overline{{W_{j}}}]},\frac{\overline{{W_{j}}}}{\frac{1}{2}[{\textstyle\textstyle\sum _{j=1}^{n}}\underline{{W_{j}}}+{\textstyle\textstyle\sum _{j=1}^{n}}\overline{{W_{j}}}]}\bigg).\]Table 4
Grey linguistic variables for determining the significance of each expert.
| Very low | Low | Medium low | Medium | Medium-high | High | Very high |
| [$1,2$] | [$2,3$] | [$3,4$] | [$4,5$] | [$5,6$] | [$6,7$] | [$7,8$] |
Step 8. To sort the grey interval numbers obtained for the weights of criteria, the order relation can be used. Another method to compare the weights of the criteria is grey possibility degree as shown in Eq. (9). To do so, the matrix of grey possibility degree is formed as:
\[G{P_{ij}}=\substack{\begin{array}{c@{\hskip4.0pt}c@{\hskip4.0pt}c@{\hskip4.0pt}c}\hspace{2em}A\hspace{54.06006pt}& B\hspace{36.98866pt}& \cdots \hspace{36.98866pt}& N\end{array}\\ {} \substack{A\\ {} B\\ {} \vdots \\ {} N}\left[\begin{array}{c@{\hskip4.0pt}c@{\hskip4.0pt}c@{\hskip4.0pt}c}P(\otimes A\otimes \leqslant A)\hspace{1em}& P(\otimes B\leqslant \otimes A)\hspace{1em}& \cdots \hspace{1em}& P(\otimes N\leqslant A\otimes )\\ {} P(\otimes A\leqslant \otimes B)\hspace{1em}& P(\otimes B\leqslant \otimes B)\hspace{1em}& \dots \hspace{1em}& P(\otimes N\leqslant \otimes B)\\ {} \vdots \hspace{1em}& \vdots \hspace{1em}& \ddots \hspace{1em}& \vdots \\ {} P(\otimes A\leqslant \otimes N)\hspace{1em}& P(\otimes B\leqslant \otimes N)\hspace{1em}& \cdots \hspace{1em}& P(\otimes N\leqslant \otimes N)\end{array}\right].}\]
Ultimately, we have the following results:
\[{P_{ij}}=\begin{array}{l}\hspace{28.45274pt}\begin{array}{c@{\hskip4.0pt}c@{\hskip4.0pt}c@{\hskip4.0pt}c}A\hspace{2em}& B\hspace{1em}\hspace{2.5pt}\hspace{2.5pt}& \cdots \hspace{2em}& N\end{array}\\ {} \substack{A\\ {} B\\ {} \vdots \\ {} N}\left[\begin{array}{c@{\hskip4.0pt}c@{\hskip4.0pt}c@{\hskip4.0pt}c}{P_{AA}}\hspace{1em}& {P_{BA}}\hspace{1em}& \dots \hspace{1em}& {P_{NA}}\\ {} {P_{AB}}\hspace{1em}& {P_{BB}}\hspace{1em}& \dots \hspace{1em}& {P_{NB}}\\ {} \vdots \hspace{1em}& \vdots \hspace{1em}& \ddots \hspace{1em}& \vdots \\ {} {P_{AN}}\hspace{1em}& {P_{BN}}\hspace{1em}& \dots \hspace{1em}& {P_{NN}}\end{array}\right],\end{array}\]
where ${D_{ij}}=\left\{\begin{array}{l@{\hskip4.0pt}l}1\hspace{1em}& P(i\leqslant j)>0.5\hspace{2.5pt}i,j=A,\dots ,N,\\ {} 0\hspace{1em}& P(i\leqslant j)\leqslant 0.5\hspace{2.5pt}i,j=A,\dots ,N.\end{array}\right.$
Table 5
Consistency Index for linguistic grey numbers.
| Linguistic terms | Equally Important (EI) | Weakly Important (WI) | Fairly Important (FI) | Very Important (VI) | Absolutely Important (AI) |
| $\otimes {a_{BW}}$ | $[1,1]$ | $[\frac{2}{3},\frac{3}{2}]$ | $[\frac{3}{2},\frac{5}{2}]$ | $[\frac{5}{2},\frac{7}{2}]$ | $[\frac{7}{2},\frac{9}{2}]$ |
| CI-GBWM | 0.00 | 0.20 | 0.71 | 1.31 | 1.96 |
By the sum of the horizontal components of the matrix ${P_{ij}}$, the scores of criteria are obtained. Based on these scores, the criteria are prioritized. Finally, based on Eq. (1), the consistency ratio can be calculated while the consistency index has been shown in Table 5. For the first time, the consistency ratios for the GBWM are calculated in the current paper. To calculate CI-GBWM, we have employed the following equation (Rezaei, 2015).
4 Practical Example and Model Validation
In this section, we implement the GBWM to solve an MCDM problem with multiple experts and analyses the important parameters in the GBWM. Then, the comparative analysis with respect to the GBWM and the fuzzy BWM is performed to verify the validity in ranking results and the advantages in keeping high reliability.
4.1 Data Collection and Implementation by the GBWM
In this section, the collected data for an MCDM problem about purchasing a car is described. Different criteria may be considered for purchasing a car. Accordingly, three experts are consulted for a group decision making. The question is which criterion is the most important and how can we find optimal weights for the criteria. The solution is using GBWM for group decision making and under uncertainty conditions.
Step 1: Based on the opinions of the experts, quality, price, comfort, safety, and style are the most important criteria for purchasing a car.
Step 2: Each expert determines the best and the worst criteria as:
Step 3: The decision-makers determine the preference degrees of the best criterion to the other criteria using the linguistic variables presented in Table 6.
Step 4. The preference of each criterion to the worst criterion is determined by the experts as shown in Table 7.
Table 6
The degree of preferences of the best criterion to the other criteria.
| Experts | Best-to-others | Price | Quality | Comfort | Safety | Style |
| 1 | Best criterion: price | EI | WI | VI | VI | AI |
| 2 | Best criterion: price | EI | FI | AI | AI | VI |
| 3 | Best criterion: quality | WI | EI | AI | VI | AI |
Step 5. Based on the information collected through Steps 3 and 4 and using the data presented in Table 3, the linguistic variables are converted into grey numbers. Then, based on Model (8), Model (9) is formed for the first decision-maker as follows:
(9)
\[\begin{aligned}{}& \min \otimes \xi ,\\ {} & \text{s.t.}\hspace{2.5pt}\otimes {W_{1}}-\bigg[\frac{2}{3},\frac{3}{2}\bigg]\otimes {W_{2}}\leqslant \otimes {\varnothing _{1}},\\ {} & -\otimes {W_{1}}+\bigg[\frac{2}{3},\frac{3}{2}\bigg]\otimes {W_{2}}\leqslant \otimes {\varnothing _{1}},\\ {} & \otimes {\varnothing _{1}}\geqslant 0,\\ {} & \otimes {\varnothing _{1}}\geqslant \otimes \xi +5\times 1\times \otimes {W_{2}}-5\times 1,\\ {} & \otimes {\varnothing _{1}}\leqslant 5\times 1\times \otimes {W_{2}},\\ {} & \otimes {\varnothing _{1}}\leqslant \otimes \xi ,\\ {} & \otimes {W_{1}}-\bigg[\frac{5}{2},\frac{7}{2}\bigg]\otimes {W_{3}}\leqslant \otimes {\varnothing _{2}},\\ {} & -\otimes {W_{1}}+\bigg[\frac{5}{2},\frac{7}{2}\bigg]\otimes {W_{3}}\leqslant \otimes {\varnothing _{2}},\\ {} & \otimes {\varnothing _{2}}\geqslant 0,\\ {} & \otimes {\varnothing _{2}}\geqslant \otimes \xi +5\times 1\times \otimes {W_{3}}-5\times 1,\\ {} & \otimes {\varnothing _{2}}\leqslant 5\times 1\times \otimes {W_{3}},\\ {} & \otimes {\varnothing _{2}}\leqslant \otimes \xi ,\\ {} & \otimes {W_{1}}-\bigg[\frac{5}{2},\frac{7}{2}\bigg]\otimes {W_{4}}\leqslant \otimes {\varnothing _{3}},\\ {} & -\otimes {W_{1}}+\bigg[\frac{5}{2},\frac{7}{2}\bigg]\otimes {W_{4}}\leqslant \otimes {\varnothing _{3}},\\ {} & \otimes {\varnothing _{3}}\geqslant 0,\\ {} & \otimes {\varnothing _{3}}\geqslant \otimes \xi +5\times 1\times \otimes {W_{4}}-5\times 1,\\ {} & \otimes {\varnothing _{3}}\leqslant 5\times 1\times \otimes {W_{4}},\\ {} & \otimes {\varnothing _{3}}\leqslant \otimes \xi ,\\ {} & \otimes {W_{1}}-\bigg[\frac{7}{2},\frac{9}{2}\bigg]\otimes {W_{5}}\leqslant \otimes {\varnothing _{4}},\end{aligned}\]
\[\begin{aligned}{}& -\otimes {W_{1}}+\bigg[\frac{7}{2},\frac{9}{2}\bigg]\otimes {W_{5}}\leqslant \otimes {\varnothing _{4}},\\ {} & \otimes {\varnothing _{4}}\geqslant 0,\\ {} & \otimes {\varnothing _{4}}\geqslant \otimes \xi +5\times 1\times \otimes {W_{5}}-5\times 1,\\ {} & \otimes {\varnothing _{4}}\leqslant 5\times 1\times \otimes {W_{5}},\\ {} & \otimes {\varnothing _{4}}\leqslant \otimes \xi ,\\ {} & \otimes {W_{1}}-\bigg[\frac{7}{2},\frac{9}{2}\bigg]\otimes {W_{5}}\leqslant \otimes {\varnothing _{5}},\\ {} & -\otimes {W_{1}}+\bigg[\frac{7}{2},\frac{9}{2}\bigg]\otimes {W_{5}}\leqslant \otimes {\varnothing _{5}},\\ {} & \otimes {\varnothing _{5}}\geqslant 0,\\ {} & \otimes {\varnothing _{5}}\geqslant \otimes \xi +5\times 1\times \otimes {W_{5}}-5\times 1,\\ {} & \otimes {\varnothing _{5}}\leqslant 5\times 1\times \otimes {W_{5}},\\ {} & \otimes {\varnothing _{5}}\leqslant \otimes \xi ,\\ {} & \otimes {W_{2}}-\bigg[\frac{5}{2},\frac{7}{2}\bigg]\otimes {W_{5}}\leqslant \otimes {\varnothing _{5}},\\ {} & -\otimes {W_{2}}+\bigg[\frac{5}{2},\frac{7}{2}\bigg]\otimes {W_{5}}\leqslant \otimes {\varnothing _{5}},\\ {} & \otimes {W_{3}}-\bigg[\frac{3}{2},\frac{5}{2}\bigg]\otimes {W_{5}}\leqslant \otimes {\varnothing _{5}},\\ {} & -\otimes {W_{3}}+\bigg[\frac{3}{2},\frac{5}{2}\bigg]\otimes {W_{5}}\leqslant \otimes {\varnothing _{5}},\\ {} & \otimes {W_{4}}-\bigg[\frac{3}{2},\frac{5}{2}\bigg]\otimes {W_{5}}\leqslant \otimes {\varnothing _{5}},\\ {} & -\otimes {W_{4}}+\bigg[\frac{3}{2},\frac{5}{2}\bigg]\otimes {W_{5}}\leqslant \otimes {\varnothing _{5}},\\ {} & \overline{{W_{1}}}-\underline{{W_{1}}}\geqslant 0.001,\hspace{1em}\overline{{W_{2}}}-\underline{{W_{2}}}\geqslant 0.001,\hspace{1em}\overline{{W_{3}}}-\underline{{W_{3}}}\geqslant 0.001,\\ {} & \overline{{W_{4}}}-\underline{{W_{4}}}\geqslant 0.001,\hspace{2em}\overline{{W_{5}}}-\underline{{W_{5}}}\geqslant 0.001,\\ {} & \otimes {W_{1}}+\otimes {W_{2}}+\otimes {W_{3}}+\otimes {W_{4}}+\otimes {W_{5}}=[0.8,1.2],\\ {} & \otimes {W_{1}}\geqslant 0,\hspace{1em}\otimes {W_{2}}\geqslant 0,\hspace{1em}\otimes {W_{3}}\geqslant 0,\hspace{1em}\otimes {W_{4}}\geqslant 0,\hspace{1em}\otimes {W_{5}}\geqslant 0.\end{aligned}\]
Table 7
The degree of preference of each criterion to the worst criterion.
| Others-to-worst | Expert | ||
| 1 | 2 | 3 | |
| Worst criterion: Style | Worst criterion: Safety | Worst criterion: Comfort | |
| Price | AI | AI | VI |
| Quality | VI | AI | AI |
| Comfort | FI | WI | EI |
| Safety | FI | EI | FI |
| Style | EI | FI | WI |
After solving Model (9) for this problem using on positioned programming approach, the weight values of the criteria for the 1st decision-maker are obtained as shown by $P1$ in Table 8. After formulating the model for the second and third decision-makers and solving them, columns $P2$ and $P3$ pertaining to the opinions of the 2nd and 3rd decision-makers are also obtained.
Table 8
The weight of each criterion based on the opinions of decision-makers.
| Variable | P1 (Medium) | P2 (Medium-low) | P3 (Low) | |||
| Upper | Lower | Upper | Lower | Upper | Lower | |
| $\otimes {W_{1}}$ | 0.3107275 | 0.2516854 | 0.4569733 | 0.3334929 | 0.4000000 | 0.2647273 |
| $\otimes {W_{2}}$ | 0.4660912 | 0.2726592 | 0.3184965 | 0.1773898 | 0.3733333 | 0.2443636 |
| $\otimes {W_{3}}$ | 0.1035758 | 0.0838951 | 0.1364985 | 0.0985499 | 0.1600000 | 0.1047273 |
| $\otimes {W_{4}}$ | 0.1864365 | 0.1078652 | 0.0969337 | 0.0638603 | 0.1600000 | 0.1047273 |
| $\otimes {W_{5}}$ | 0.1331689 | 0.0838951 | 0.1910979 | 0.1267070 | 0.1066667 | 0.0814545 |
The opinions of the experts regarding each criterion have been compared in Fig. 1. According to Fig. 1, ‘price’ has the highest score by the first two decision-makers, which seems rational considering the fact that it has also been selected as the ‘best’ criterion. Also, the 3rd decision-maker has selected ‘quality’ as the best criterion. Similarly, the worst criterion for each decision-maker can be observed in Fig. 1. Overall, Fig. 1 demonstrates the differences between the viewpoints of experts with regard to each criterion.
Step 6. Using Eq. (30) and according to the data presented in Table 8, the weight of each criterion is presented in Table 9.
Step 7. The normalized weights are presented in Table 10. Also, Fig. 2 represents the final obtained weights.
Table 9
The aggregated weights based on the opinions of decision-makers.
| Variable | Lower | Upper |
| $\otimes {W_{1}}$ | 0.282710725 | 0.39256978 |
| $\otimes {W_{2}}$ | 0.22503049 | 0.37427681 |
| $\otimes {W_{3}}$ | 0.09769076 | 0.13611538 |
| $\otimes {W_{4}}$ | 0.08939217 | 0.14066120 |
| $\otimes {W_{5}}$ | 0.09500056 | 0.13694033 |
Table 10
Final normalized weights.
| Variable | Lower bound | Upper bound |
| $\otimes {W_{1}^{\ast }}$ | 0.286965 | 0.3984682 |
| $\otimes {W_{2}^{\ast }}$ | 0.228412 | 0.3799003 |
| $\otimes {W_{3}^{\ast }}$ | 0.099159 | 0.1381605 |
| $\otimes {W_{4}^{\ast }}$ | 0.090735 | 0.1427746 |
| $\otimes {W_{5}^{\ast }}$ | 0.096428 | 0.1389979 |
Step 8. Based on the order relation of grey numbers, the criteria are prioritized as:
Another method which yields a relatively similar result is the formation of a matrix of grey possibility degree which is concluded as Eq. (32).
(32)
\[\begin{aligned}{}& G{P_{ij}}=\left[\begin{array}{c@{\hskip4.0pt}c@{\hskip4.0pt}c@{\hskip4.0pt}c@{\hskip4.0pt}c}0.50\hspace{1em}& 0.64\hspace{1em}& 1.00\hspace{1em}& 1.00\hspace{1em}& 1.00\\ {} 0.35\hspace{1em}& 0.50\hspace{1em}& 1.00\hspace{1em}& 1.00\hspace{1em}& 1.00\\ {} 0.00\hspace{1em}& 0.00\hspace{1em}& 0.50\hspace{1em}& 0.52\hspace{1em}& 0.51\\ {} 0.00\hspace{1em}& 0.00\hspace{1em}& 0.47\hspace{1em}& 0.50\hspace{1em}& 0.48\\ {} 0.00\hspace{1em}& 0.00\hspace{1em}& 0.48\hspace{1em}& 0.51\hspace{1em}& 0.50\end{array}\right],\end{aligned}\](33)
\[\begin{aligned}{}& {P_{ij}}=\left[\begin{array}{c@{\hskip4.0pt}c@{\hskip4.0pt}c@{\hskip4.0pt}c@{\hskip4.0pt}c}0.00\hspace{1em}& 1.00\hspace{1em}& 1.00\hspace{1em}& 1.00\hspace{1em}& 1.00\\ {} 0.00\hspace{1em}& 0.00\hspace{1em}& 1.00\hspace{1em}& 1.00\hspace{1em}& 1.00\\ {} 0.00\hspace{1em}& 0.00\hspace{1em}& 0.00\hspace{1em}& 1.00\hspace{1em}& 1.00\\ {} 0.00\hspace{1em}& 0.00\hspace{1em}& 0.00\hspace{1em}& 0.00\hspace{1em}& 0.00\\ {} 0.00\hspace{1em}& 0.00\hspace{1em}& 0.00\hspace{1em}& 1.00\hspace{1em}& 0.00\end{array}\right]\substack{4\\ {} 3\\ {} 2\\ {} 0\\ {} 1}\end{aligned}\]Based on the sum of the horizontal components of the matrix in Eq. (33), the prioritization of the criteria is concluded as $\mathrm{Price}>\mathrm{Quality}>\mathrm{Comfort}>\mathrm{Style}>\mathrm{Safety}$.
4.2 Sensitivity Analysis
In this section, we aim to analyse the sensitivity of the example solved in the previous section. A primary reason for undertaking sensitivity analysis is to investigate the changes of output parameters resulting from changes in the input data. Since the current study makes use of the positioned programming approach, the parameter β is considered as constant in the sensitivity analysis. The sensitivity of the consistency ratio is calculated according to the changes of the parameters ρ and δ. In other words, we are going to examine the impact of parameters ρ and δ on consistency ratio.
As is depicted in Fig. 3, the parameter ρ in the positioned programming is indifferent to changes in the consistency ratio. The parameter δ has a direct impact on the consistency ratio, such that any increase in this parameter would result in an increase in the consistency ratio. Fig. 3 demonstrates the degree of sensitivity of the 1st decision-maker.
It can be seen from Fig. 4 that the consistency ratio has a direct relationship with the parameter δ in the positioned programming. The consistency ratio is indeterminate to the parameter ρ. Fig. 4 is based on the viewpoint of the 2nd decision-maker.
According to Fig. 5, the consistency ratio has a direct correlation with the parameter δ in the positioned programming. If a need arises, the definitive response could be achieved by determining the parameters ρ, β and δ based on the positioned programming method. The lower the parameter δ is and the higher the parameters ρ and β are, the lower the inconsistency response ratio will be achieved. Determining suitable parameters for positioned programming is a challenge as well. Finding the optimal values for these parameters is not an easy job and it depends on different factors. In future researches, Scholar can provide a methodology for finding optimal values. This is another benefit of GBWM that can present a crisp solution based on expert’s needs in different conditions.
4.3 Comparative Analyses and Discussions
In this section, we solve the numerical example by the FBWM (Guo and Zhao, 2017) for each expert to analyse the weights of criteria and the consistency index compared with the GBWM proposed in this paper. We tabulate the detailed computation results of the three experts in Tables 11–13. It should be noted that the results of GBWM have been obtained based on ideal and critical solutions using on position programming method while we have $(\rho ,\beta ,\delta )=(0,0,1)$ and $(\rho ,\beta ,\delta )=(1,1,0)$. These two situations include the highest and lowest consistency ratio for the problem.
Table 11
The deduced weights of criteria and consistency index by the FBWM and GBWM for expert 1.
| FBWM (Guo and Zhao, 2017) CR = 0.2180 | GBWM (This paper) CR = $[0,0.0519]$ | ||||
| Weights | Fuzzy weights | Defuzzified | Rank | Grey weights | Rank |
| W1 | $(0.3723,0.3723,0.3771)$ | 0.3731 | 1 | $[0.2647,0.4000]$ | 1 |
| W2 | $(0.1965,0.2609,0.3403)$ | 0.2634 | 2 | $[0.2444,0.3733]$ | 2 |
| W3 | $(0.1017,0.1363,0.1796)$ | 0.1377 | 3 | $[0.1047,0.1600]$ | 3 |
| W4 | $(0.1017,0.1363,0.1796)$ | 0.1377 | 3 | $[0.1047,0.1600]$ | 3 |
| W5 | $(0.0867,0.0867,0.0948)$ | 0.0880 | 4 | $[0.0815,0.1067]$ | 4 |
Table 12
The deduced weights of criteria and consistency index by the FBWM and GBWM for expert 2.
| FBWM (Guo and Zhao, 2017) CR = 0.3426 | GBWM (This paper) CR = $[0.0106,0.0561]$ | ||||
| Weights | Fuzzy weights | Defuzzified | Rank | Grey weights | Rank |
| W1 | $(0.3271,0.3663,0.4055)$ | 0.3663 | 1 | $[0.3335,0.4570]$ | 1 |
| W2 | $(0.2218,0.2757,0.3948)$ | 0.2866 | 2 | $[0.1774,0.3185]$ | 2 |
| W3 | $(0.1049,0.1100,0.1156)$ | 0.1101 | 4 | $[0.0985,0.1365]$ | 4 |
| W4 | $(0.0784,0.0784,0.0784)$ | 0.0784 | 5 | $[0.0639,0.0969]$ | 5 |
| W5 | $(0.1434,0.1573,0.1789)$ | 0.1586 | 3 | $[0.1267,0.1911]$ | 3 |
Table 13
The deduced weights of criteria and consistency index by the FBWM and GBWM for expert 3.
| FBWM (Guo and Zhao, 2017) CR = 0.2180 | GBWM (This paper) CR = $[0,0.0535]$ | ||||
| Weights | FBWM weights | Defuzzified | Rank | GBWM weights | Rank |
| W1 | $(0.2165,0.2681,0.3498)$ | 0.2731 | 2 | $[0.2517,0.3107]$ | 2 |
| W2 | $(0.3827,0.3827,0.4173)$ | 0.3884 | 1 | $[0.2727,0.4661]$ | 1 |
| W3 | $(0.0891,0.0891,0.0974)$ | 0.0905 | 5 | $[0.0839,0.1036]$ | 5 |
| W4 | $(0.1062,0.1429,0.1846)$ | 0.1438 | 3 | $[0.1864,0.1079]$ | 3 |
| W5 | $(0.0847,0.1043,0.1234)$ | 0.1042 | 4 | $[0.1332,0.0839]$ | 4 |
Based on the information given in Tables 11, 12, 13, we can obtain some characteristics of the GBWM compared with the FBWM.
-
(1) The ranking result computed by the GBWM is valid considering the same ranking of alternatives for different input values.
-
(2) The result calculated by the GBWM is more reliable than the one calculated by the FBWM because the GBWM has a smaller inconsistency ratio compared with that of the FBWM. Fig. 6 shows a comparison of the consistency ratio for the three experts in the numerical example.
-
(3) After defuzzifications of triangular fuzzy weights of criteria, the defuzzied weights belong to the grey weights deduced by the GBWM. The GBWM narrows the feasible space for potential weights of criteria by the positioned programming to obtain the reliability of weights. It can be shown in Fig. 7.
It is interesting to mention that the linguistic variables for both GBWM and FBWM are same, but GBWM employs a grey linear model as a core model and that is why the results are more reliable than of FBWM method. On the other hand, employing grey system theory can contribute to decrease the volume of calculations. Therefore, it is clear why current research suggests using GBMW method. Moreover, GBWM method can provide a crisp solution for decision-maker if there is a need. Decision-makers should just provide suitable values for the parameters in a grey linear model and this is another advantage of GBMW method. In conclusion, the GBWM is valid in deducing the weights of criteria and advantageous in keeping reliable results and requiring less computational complexity.
5 Conclusions
BWM has shown an acceptable performance, but the fact that it does not consider the uncertainties of the decision-making environment may reduce the performance of the BWM in real-world. In reality, most information is not clear and ambiguous. The grey system theory is a suitable approach for taking into account the uncertainties of the decision-making environment. In this regard, the current study presented the GBWM. The result calculated by the GBWM is more reliable than that calculated by the FBWM because the GBWM has a smaller inconsistency ratio compared with FBWM. Moreover, the GBWM uses a linear model that is able to present global-optimum weights for MCDM problems. It is interesting to mention that GBWM method can provide a crisp solution based on the experts’ needs. However, some limitations include determining suitable parameters for using on position programming method during solving grey linear programming, which it could be examined by scholars in future. Also, future studies may use grey-fuzzy hybrid approaches for the BWM and compare the results with those derived by the GBWM and FBWM methods. Furthermore, scholars can use different methods for linearization of the core model of the BWM and it may provide a better solution for MCDM problems.
