Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Multi-criteria group decision making based on Archimedean power partitioned Muirhead mean operators of q-rung orthopair fuzzy numbers

  • Yuchu Qin,

    Roles Conceptualization, Software, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation School of Computing and Engineering, University of Huddersfield, Huddersfield, West Yorkshire, England, United Kingdom

  • Qunfen Qi ,

    Roles Conceptualization, Funding acquisition, Writing – review & editing

    q.qi@hud.ac.uk

    Affiliation School of Computing and Engineering, University of Huddersfield, Huddersfield, West Yorkshire, England, United Kingdom

  • Paul J. Scott,

    Roles Methodology, Supervision, Writing – review & editing

    Affiliation School of Computing and Engineering, University of Huddersfield, Huddersfield, West Yorkshire, England, United Kingdom

  • Xiangqian Jiang

    Roles Supervision, Writing – review & editing

    Affiliation School of Computing and Engineering, University of Huddersfield, Huddersfield, West Yorkshire, England, United Kingdom

Abstract

Two critical tasks in multi-criteria group decision making (MCGDM) are to describe criterion values and to aggregate the described information to generate a ranking of alternatives. A flexible and superior tool for the first task is q-rung orthopair fuzzy number (qROFN) and an effective tool for the second task is aggregation operator. So far, nearly thirty different aggregation operators of qROFNs have been presented. Each operator has its distinctive characteristics and can work well for specific purpose. However, there is not yet an operator which can provide desirable generality and flexibility in aggregating criterion values, dealing with the heterogeneous interrelationships among criteria, and reducing the influence of extreme criterion values. To provide such an aggregation operator, Muirhead mean operator, power average operator, partitioned average operator, and Archimedean T-norm and T-conorm operations are concurrently introduced into q-rung orthopair fuzzy sets, and an Archimedean power partitioned Muirhead mean operator of qROFNs and its weighted form are presented and a MCGDM method based on the weighted operator is proposed in this paper. The generalised expressions of the two operators are firstly defined. Their properties are explored and proved and their specific expressions are constructed. On the basis of the specific expressions, a method for solving the MCGDM problems based on qROFNs is then designed. Finally, the feasibility and effectiveness of the method is demonstrated via a numerical example, a set of experiments, and qualitative and quantitative comparisons.

1. Introduction

Multi-criteria group decision making (MCGDM) refers to the process of finding optimal alternatives in complex scenarios via synthetically evaluating the values of multiple criteria of all alternatives provided by a group of domain experts [1]. In this process, there are two critical tasks. One critical task is to describe the values of different criteria accurately and effectively. For such description, there are many different kinds of available tools, where fuzzy set is a well-known kind [212]. To date, nearly thirty different types of fuzzy sets have been presented [13]. Yager’s generalised orthopair fuzzy set [14], commonly known as q-rung orthopair fuzzy set (qROFS), is one of the most important and popular types among them.

A qROFS consists of an element and a q-rung orthopair membership grade (commonly known as q-rung orthopair fuzzy number (qROFN), which is used to quantify the degrees of membership and non-membership of the element to the qROFS. In a qROFS, both the degrees of membership and non-membership and the sum of the q-th (q = 1, 2, 3, …) power of the degree of membership and the q-th power of the degree of non-membership are restricted to [0, 1]. In other words, the rung q in a qROFS is adjustable under the premise of satisfying this condition. Because of this characteristic, qROFS can be regarded as the generalisation of Zadeh’s fuzzy set (FS) [15], Atanassov’s intuitionistic fuzzy set (IFS) [16], and Yager’s Pythagorean fuzzy set (PFS) [17]. This is because qROFS will reduce to FS when q = 1 and the sum of the degrees of membership and non-membership is equal to 1, will reduce to IFS when q = 1, and will reduce to PFS when q = 2. In addition, the expressiveness of a qROFS will continue to increase as q increases, which provides enough freedom for the description of fuzzy information. Due to such advantage, qROFSs have received extensive attention in the field of MCGDM during the past few years. Various research topics regarding qROFSs for MCGDM, such as correlation and correlation coefficient of qROFSs [18], distance measures of qROFSs [19], similarity measures of qROFSs [20], application of qROFSs in practical MCGDM problems [21, 22], operational rules of qROFNs [2325], and aggregation operators of qROFNs [26−38], are gaining importance and popularity within academia.

The other critical task in MCGDM is to fuse the described criterion information to generate a ranking of all alternatives. For such fusion, aggregation operator is regarded as an effective tool. So far, over twenty different aggregation operators of qROFNs have been presented. They are the weighted averaging (WA) operator and the weighted geometric (WG) operator presented by Liu and Wang [26], the weighted Bonferroni mean (WBM) operator and the weighted geometric Bonferroni mean (WGBM) operator presented by Liu and Liu [27], the weighted Archimedean Bonferroni mean (WABM) operators presented by Liu and Wang [28], the weighted partitioned Bonferroni mean (WPBM) operator and the weighted partitioned geometric Bonferroni mean (WPGBM) operator presented by Yang and Pang [29], the weighted Heronian mean (WHM) operator and the weighted geometric Heronian mean (WGHM) operator presented by Wei et al. [30], the WHM* operator (This operator is different from the WHM operator presented by Wei et al. [30], though they have the same names) and the weighted partitioned Heronian mean (WPHM) operator presented by Liu et al. [31], the weighted Maclaurin symmetric mean (WMSM) operator and the weighted geometric Maclaurin symmetric mean (WGMSM) operator presented by Wei et al. [32], the weighted power Maclaurin symmetric mean (WPMSM) operator presented by Liu et al. [33], the weighted power partitioned Maclaurin symmetric mean (WPPMSM) operator presented by Bai et al. [34], the weighted Muirhead mean (WMM) operator and the weighted geometric Muirhead mean (WGMM) operator presented by Wang et al. [35], the weighted extended Bonferroni mean (WEBM) operator presented by Liu et al. [36], the weighted exponential (WE) operator presented by Peng et al. [37], and the weighted point (WP) operators presented by Xing et al. [38]. Each operator has its own characteristics and can work well for its specific purpose. But there is not yet an operator that has the following three characteristics at the same time: (1) Provide satisfying generality and flexibility in the aggregation of qROFNs; (2) Deal with the situation in which the criteria are divided into several parts and there are interrelationships among different criteria in each part whereas the criteria in different parts are independent of each other; (3) Reduce the negative effect of the unduly high or unduly low criterion values on the aggregation results.

In practical MCGDM problems, aggregation of criterion values is a complex process, in which the preferences of decision makers may change frequently. An ideal aggregation operator should be general and flexible enough to adapt to such change. Moreover, there are usually complex relationships among the different criteria considered in the problems. It is also of importance for an aggregation operator to capture the complex interrelationships of different criteria to generate more reasonable aggregation results [28]. Further, the values of criteria are generally assessed by domain experts. It is often difficult to ensure the absolute objectivity, which means that a few biased experts will give biased assessment values [33]. To obtain reasonable aggregation results, it is of necessity to reduce the negative influence of biased criterion values in the aggregation. Based on these considerations, the motivations of the present paper are explained as follows:

  1. To develop an aggregation operator of qROFNs which can capture the complex interrelationships among criteria, the Muirhead mean (MM) operator [39] and partitioned average operator are introduced. The MM operator, which is a generalisation of the generalised arithmetic average operator, Bonferroni mean (BM) operator [40, 41], Maclaurin symmetric mean (MSM) operator [42, 43], and generalised geometric average operator, is an all-in-one aggregation operator for capturing the interrelationships of criteria. It is applicable in the cases where all criteria are independent of each other, where there are interrelationships between any two criteria, and where there are interrelationships between any multiple (three or more) criteria [4449]. The partitioned average operator is an aggregation operator that has the capability of aggregating the parameters in different partitions using the same aggregation operator and aggregating the aggregation results of different partitions using the arithmetic average operator [5055].
  2. To enable the aggregation operator to reduce the negative effect of extreme criterion values on the aggregation results, the power average (PA) operator [56] is combined into the partitioned MM operator. The PA operator is an aggregation operator that can assign weights to the aggregated parameters by calculating the support degrees between these parameters, which makes it capable of reducing the negative influence of unreasonable parameter values on the aggregation results [5762].
  3. To improve the generality and flexibility of the combined aggregation operator, the operational rules of qROFNs based on the Archimedean T-norm and T-conorm (ATT) are leveraged to perform the operations in the operator. The ATT are operations for generalising the logical conjunction and disjunction to fuzzy logic. They are important tools that can generate versatile and flexible operational rules for fuzzy numbers and the aggregation operators based on them are rather general and flexible for aggregating fuzzy information [28, 6368].

As can be summarised from the motivations above, this paper aims to present a set of Archimedean power partitioned MM operators of qROFNs and propose a MCGDM method based on them. This aim is achieved through combining the MM operator, the partitioned average operator, the PA operator, and the ATT. As a result, the presented aggregation operators combine all of their characteristics.

The remainder of the paper is organised as follows. A brief introduction of some basic concepts is provided Section 2. Sections 3 explains the details of the presented Archimedean power partitioned MM operators. The specific process of the proposed MCGDM method is described in Section 4. Section 5 demonstrates and evaluates the presented operators and proposed method via example, experiments, and comparisons. Section 6 ends the paper with a conclusion.

2. Preliminaries

In this section, some prerequisites in qROFS theory, operational rules of qROFNs based on ATT, PA operator, and MM operator are briefly introduced to facilitate the understanding of the present paper.

2.1. qROFS theory

qROFS [14] is the generalisation of FS [15], IFS [16], and PFS [17]. Its formal definition is as follow:

Definition 1 [14]. A qROFS S in a finite universe of discourse X is S = {<x, μS(x), νS(x)> | xX}, where μS: X → [0, 1] is the degree of membership of xX to S, and νS: X → [0, 1] is the degree of non-membership of xX to S, such that 0 ≤ (μS(x))q + (νS(x))q ≤ 1 (q = 1, 2, 3,…). The degree of hesitancy of xX to S is πS(x) = (1 − (μS(x))q − (νS(x))q)1/q.

For convenience, a pair <μS(x), νS(x)> is called a qROFN, which is commonly simplified as Q = <μ, ν>. To compare two qROFNs, their scores and accuracies are required, which can be calculated according to the following definitions:

Definition 2 [26]. Let Q = <μ, ν> be a qROFN. Then its score is S(Q) = μqνq. Obviously, −1 ≤ S(Q) ≤ 1.

Definition 3 [26]. Let Q = <μ, ν> be a qROFN. Then its accuracy is A(Q) = μq + νq. Obviously, 0 ≤ A(Q) ≤ 1.

Using S(Q) and A(Q), two qROFNs can be compared via the following definition:

Definition 4 [26]. Let Q1 = <μ1, ν1> and Q2 = <μ2, ν2> be any two qROFNs, S(Q1) and S(Q2) be respectively the scores of Q1 and Q2, and A(Q1) and A(Q2) be respectively the accuracies of Q1 and Q2. Then: (1) If S(Q1) > S(Q2), then Q1 > Q2; (2) If S(Q1) = S(Q2) and A(Q1) > A(Q2), then Q1 > Q2; (3) If S(Q1) = S(Q2) and A(Q1) = A(Q2), then Q1 = Q2.

To calculate the distance between two qROFNs, a distance measure of qROFNs is required. The following definition provides the Minkowski-type distance measure of qROFNs:

Definition 5 [19]. Let Q1 = <μ1, ν1> and Q2 = <μ2, ν2> be any two qROFNs. Then the Minkowski-type distance between Q1 and Q2 is D(Q1, Q2) = (0.5|μ1μ2|b + 0.5|ν1ν2|b)1/b (b = 1, 2, …).

If b = 1, the Minkowski-type distance between Q1 and Q2 will reduce to the Hamming distance between Q1 and Q2: D(Q1, Q2) = 0.5(|μ1μ2| + |ν1ν2|); If b = 2, the Minkowski-type distance between Q1 and Q2 will reduce to the Euclidean distance between Q1 and Q2: (0.5|μ1μ2|2 + 0.5|ν1ν2|2)1/2; If b = ∞, the Minkowski-type distance between Q1 and Q2 will reduce to the Chebyshev distance between Q1 and Q2: D(Q1, Q2) = max{|μ1μ2|, |ν1ν2|}.

2.2. Operational rules of qROFNs based on ATT

Based on ATT, a set of general and versatile operational rules of qROFNs were established by Liu and Wang [28]. The formal definition of the rules is as follow:

Definition 6 [28]. Let Q = <μ, ν>, Q1 = <μ1, ν1>, and Q2 = <μ2, ν2> be any three qROFNs, and a and b be any two real numbers and a, b > 0. Then the sum, product, multiplication, and power operations of qROFNs based on the Archimedean T-norm T(x, y) = f−1(f(x) + f(y)) and its T-conorm T*(x, y) = g−1(g(x) + g(y)) can be respectively defined as follows: (1) (2) (3) (4)

2.3. PA operator

The PA operator, introduced by Yager [56], can assign weights to the aggregated arguments via computing the degrees of support between these arguments. This makes it possible to reduce the negative influence of the unduly high or unduly low argument values on the aggregation results. The formal definition of the PA operator is as follow:

Definition 7 [56]. Let (a1, a2, …, an) be a collection of crisp numbers, Sup(ai, aj) = 1 − D(ai, aj) (i, j = 1, 2, …, n and ji; D(ai, aj) is the distance between ai and aj) be the degree of support for ai from aj which has the following properties: (1) 0 ≤ Sup(ai, aj) ≤ 1; (2) Sup(ai, aj) = Sup(aj, ai); (3) Sup(ai, aj) ≥ Sup(ar, as) if |aiaj| ≤ |aras|, and (5) Then the aggregation function (6) is called the PA operator.

2.4. Partitioned average operator

The partitioned average operator can aggregate the arguments in different partitions using the same aggregation operator and aggregate the aggregation results of different partitions using the arithmetic average operator [51]. Its formal definition is as follow:

Definition 8 [51]. Let (a1, a2, …, an) be a collection of crisp numbers, S = {a1, a2, …, an} be a set of a1, a2, …, an, Sk = {a1, a2, …, a|Sk|} (k = 1, 2, …, N) be N partitions of S (i.e. S1S2 ∪ … ∪ SN = S and S1S2 ∩ … ∩ SN = Ø), and AO be a specific aggregation operator. Then the aggregation function (7) is called the partitioned average operator.

2.5. MM operator

The MM operator was firstly introduced to aggregate crisp numbers by Muirhead [39]. It has prominent characteristics in capturing the interrelationships among multiple aggregated arguments and providing a general form of a number of other aggregation operators. The formal definition of the MM operator is as follow:

Definition 9 [39]. Let (a1, a2, …, an) be a collection of crisp numbers, Δ = (δ1, δ2, …, δn) (where δ1, δ2, …, δn ≥ 0 but not at the same time δ1 = δ2 = … = δn = 0) be a collection of n real numbers, p(i) be a permutation of (1, 2,…, n), and Pn be the set of all permutations of (1, 2,…, n). Then the aggregation function (8) is called the MM operator.

In this operator, whether the interrelationships are considered depends on the values of δi (i = 1, 2, …, n): (1) If δ1 > 0 and δ2 = δ3 = … = δn = 0, then the interrelationships are not considered; (2) If δ1, δ2 > 0 and δ3 = δ4 = … = δn = 0, then the interrelationships between two crisp numbers are considered; (3) If δ1, δ2, …, δk > 0 (k = 3, 4, …, n) and δk+1 = δk+2 = … = δn = 0, then the interrelationships among k crisp numbers are considered.

3. Aggregation operators

In this section, a q-rung orthopair fuzzy Archimedean power partitioned MM (qROFAPPMM) operator and a q-rung orthopair fuzzy weighted Archimedean power partitioned MM (qROFWAPPMM) operator are presented. The properties of the two operators are explored and their specific cases are discussed.

3.1. qROFAPPMM operator

A qROFAPPMM operator is a power partitioned MM operator of qROFNs, in which the sum, product, multiplication, and power operations are performed using the operational rules of qROFNs based on ATT. Its formal definition is as follow:

Definition 10. Let Q1, Q2, …, Qn (Qi = <μi, νi>, i = 1, 2, …, n) be n qROFNs (q = 1, 2, 3, …), (Q1, Q2, …, Qn) be a collection of Q1, Q2, …, Qn, S = {Q1, Q2, …, Qn} be an ordered set of Q1, Q2, …, Qn, Sk = {Q1, Q2, …, Q|Sk|} (k = 1, 2, …, N) be N partitions of S (i.e. S1S2 ∪ … ∪ SN = S and S1S2 ∩ … ∩ SN = Ø), δ1, δ2, …, δ|Sk| (k = 1, 2, …, N and δ1, δ2, …, δ|Sk| ≥ 0 but not at the same time δ1 = δ2 = … = δ|Sk| = 0) be |Sk| real numbers that respectively correspond to Q1, Q2, …, Q|Sk|, Δk = (δ1, δ2, …, δ|Sk|) be a collection of δ1, δ2, …, δ|Sk|, Δ = (Δ1, Δ2, …, ΔN) be a collection of Δ1, Δ2, …, ΔN, p(ik) be a permutation of (1, 2,…, |Sk|), P|Sk| be a set of all permutations of (1, 2,…, |Sk|), QiQj and QiQj (i, j = 1, 2, …, n) be respectively the sum and product operations of Qi and Qj based on ATT, aQr and Qsb (r, s = 1, 2, …, n; a, b > 0) be respectively the multiplication operation of Qr and the power operation of Qs based on ATT, Sup(Qr, Qs) = 1 − D(Qr, Qs) (r, s = 1, 2, …, n and sr; D(Qr, Qs) is the distance between Qr and Qs) be the degree of support for Qr from Qs which satisfy 0 ≤ Sup(Qr, Qs) ≤ 1, Sup(Qr, Qs) = Sup(Qs, Qr), and Sup(Qr, Qs) ≥ Sup(Qu, Qv) if |QrQs| ≤ |QuQv|, and (9) Then the aggregation function (10) is called the qROFAPPMM operator. In this operator, the values of δik (ik = 1, 2, …, |Sk|) are used to capture the interrelationships among the aggregated qROFNs in each of the N partitions S: (1) If δ1 > 0 and δ2 = δ3 = … = δ|Sk| = 0, then the qROFNs in the k-th partition Sk are independent of each other; (2) If δ1, δ2 > 0 and δ3 = δ4 = … = δ|Sk| = 0, then the interrelationships between two qROFNs in Sk are considered; (3) If δ1, δ2, …, δr > 0 (r = 3, 4, …, n) and δr+1 = δr+2 = … = δ|Sk| = 0, then the interrelationships among r qROFNs in Sk are considered.

According to Eq (1)‒(4) and (10), the following theorem is obtained:

Theorem 1. Let Q1, Q2, …, Qn (Qi = <μi, νi>, i = 1, 2, …, n) be n qROFNs (q = 1, 2, 3, …) and (Q1, Q2, …, Qn) be a collection of Q1, Q2, …, Qn. Then (11) and it is still a qROFN, where (12) (13) and ξp(ik) is a PA factor which can be calculated via the following Equation: (14)

For the details regarding the proof of this theorem, please refer to Appendix A in S1 File. The following two theorems respectively state the idempotency and boundedness of the qROFAPPMM operator:

Theorem 2. Let Q1, Q2, …, Qn (Qi = <μi, νi>, i = 1, 2, …, n) be n qROFNs (q = 1, 2, 3, …) and (Q1, Q2, …, Qn) be a collection of Q1, Q2, …, Qn. If Qi = Q = <μQ, νQ> for all i = 1, 2, …, n, then qROFAPPMMΔ(Q1, Q2, …, Qn) = Q = <μQ, νQ>.

Theorem 3. Let Q1, Q2, …, Qn (Qi = <μi, νi>, i = 1, 2, …, n) be n qROFNs (q = 1, 2, 3, …), (Q1, Q2, …, Qn) be a collection of Q1, Q2, …, Qn, QUB = <max{μi}, min{νi}>, and QLB = <min{μi}, max{νi}>. Then QLBqROFAPPMMΔ(Q1, Q2, …, Qn) ≤ QUB.

For the details regarding the proofs of these two theorems, please refer to Appendixes B and C in S1 File, respectively.

From Eq (11), it is not difficult to see that the qROFAPPMM operator has the capabilities to capture the interrelationships of the aggregated qROFNs and to reduce the negative influence of unreasonable aggregated qROFNs, because the operator is constructed via combining the MM, partitioned average, and PA operators. In addition, it can also be seen that the qROFAPPMM operator has desirable generality, since different specific qROFAPPMM operators will be obtained if different if specific functions are assigned to f. For example, if the additive generators of Algebraic, Einstein, Hamacher, and Frank T-norms and T-conorms are assigned to f, then four specific operators can be respectively constructed as follows:

  1. If f(t) = −Intq, then g(t) = −In(1−tq), f−1(t) = (et)1/q, and g−1(t) = (1−et)1/q. A q-rung orthopair fuzzy Archimedean Algebraic power partitioned MM (qROFAAPPMM) operator is constructed according to Eq (11): (15) where ξp(ik) is a PA factor which can be calculated via Eq (14).
  2. If f(t) = In[(2−tq)/tq], then g(t) = In[(1+tq)/(1−tq)], f−1(t) = [2/(et+1)]1/q, and g−1(t) = [(et−1)/(et+1)]1/q. A q-rung orthopair fuzzy Archimedean Einstein power partitioned MM (qROFAEPPMM) operator is constructed according to Eq (11): (16) where (17) (18) and ξp(ik) is a PA factor which can be calculated via Eq (14).
  3. If f(t) = In{[λ+(1−λ)tq]/tq} (λ > 0), then g(t) = In{[λ+(1−λ)(1−tq)]/(1−tq)}, f−1(t) = [λ/(et+λ−1)]1/q, and g−1(t) = [(et−1)/ (et+λ−1)]1/q. A q-rung orthopair fuzzy Archimedean Hamacher power partitioned MM (qROFAHPPMM) operator is constructed according to Eq (11): (19) where (20) (21) and ξp(ik) is a PA factor which can be calculated via Eq (14).
  4. If f(t) = −In[(ε−1)/(εy−1)] (y = tq; ε > 1), then g(t) = −In[(ε−1)/(ε1−y−1)], f−1(t) = {logε[(ε−1+et)/et]}1/q, and g−1(t) = {1−logε[(ε−1+et)/et]}1/q. A q-rung orthopair fuzzy Archimedean Frank power partitioned MM (qROFAFPPMM) operator is constructed according to Eq (11): (22) where (23) (24) (25) (26) (27) (28) (29) (30) and ξp(ik) is a PA factor which can be calculated via Eq (14).

3.2. qROFWAPPMM operator

The qROFAPPMM operator has advantages in having satisfying generality and flexibility, capturing the complex interrelationships among qROFNs, and reducing the negative effect of unreasonable qROFNs on the aggregation results. But it does not consider the relative importance of each aggregated qROFN. To overcome this limitation, weights are introduced and a qROFWAPPMM operator is presented. The formal definition of this operator is as follow:

Definition 11. On the basis of Definition 10, let w1, w2, …, wn be respectively the weights of Q1, Q2, …, Qn such that 0 ≤ w1, w2, …, wn ≤ 1 and w1+w2+…+wn = 1. Then the aggregation function (31) is called the qROFWAPPMM operator. In this operator, the functions of δik (ik = 1, 2, …, |Sk|) are the same as the functions of δik in the qROFAPPMM operator (see Eq (11)).

According to Eqs (1)‒(4) and (31), the following theorem is obtained:

Theorem 4. Let Q1, Q2, …, Qn (Qi = <μi, νi>, i = 1, 2, …, n) be n qROFNs (q = 1, 2, 3, …) and (Q1, Q2, …, Qn) be a collection of Q1, Q2, …, Qn. Then (32) and it is still a qROFN, where (33) (34) and ξp(ik) and ξt are two PA factors which can be calculated via Eq (14).

The proof of Theorem 4 is similar to the proof of Theorem 1 (see Appendix A in S1 File) and is omitted here. It is worth nothing that the qROFWAPPMM operator no longer has the properties of idempotency and boundedness due to the introduce of weights.

For Eq (32), if the additive generators of Algebraic, Einstein, Hamacher, and Frank T-norms and T-conorms are assigned to f, then four specific operators can be respectively constructed as follows:

  1. If f(t) = −Intq, then g(t) = −In(1−tq), f−1(t) = (et)1/q, and g−1(t) = (1−et)1/q. A q-rung orthopair fuzzy weighted Archimedean Algebraic power partitioned MM (qROFWAAPPMM) operator is constructed according to Eq (32): (35) where ξp(ik) and ξt are two PA factors which can be calculated via Eq (14).
  2. If f(t) = In[(2−tq)/tq], then g(t) = In[(1+tq)/(1−tq)], f−1(t) = [2/(et+1)]1/q, and g−1(t) = [(et−1)/(et+1)]1/q. A q-rung orthopair fuzzy weighted Archimedean Einstein power partitioned MM (qROFWAEPPMM) operator is constructed according to Eq (32): (36) where (37) (38) and ξp(ik) and ξt are two PA factors which can be calculated via Eq (14).
  3. If f(t) = In{[λ+(1−λ)tq]/tq} (λ > 0), then g(t) = In{[λ+(1−λ)(1−tq)]/(1−tq)}, f−1(t) = [λ/(et+λ−1)]1/q, and g−1(t) = [(et−1)/ (et+λ−1)]1/q. A q-rung orthopair fuzzy weighted Archimedean Hamacher power partitioned MM (qROFWAHPPMM) operator is constructed according to Eq (32): (39) where (40) (41) and ξp(ik) and ξt are two PA factors which can be calculated via Eq (14).
  4. If f(t) = −In[(ε−1)/(εy−1)] (y = tq; ε > 1), then g(t) = −In[(ε−1)/(ε1−y−1)], f−1(t) = {logε[(ε−1+et)/et]}1/q, and g−1(t) = {1− logε[(ε−1+et)/et]}1/q. A q-rung orthopair fuzzy weighted Archimedean Frank power partitioned MM (qROFWAFPPMM) operator is constructed according to Eq (32): (42) where (43) (44) (45) (46) (47) (48) (49) (50) and ξp(ik) and ξt are two PA factors which can be calculated via Eq (14).

4. MCGDM method

In this section, a MCGDM method based on the qROFWAPPMM operator is proposed to solve the MCGDM problems based on qROFNs.

In general, a MCGDM problem based on qROFNs can be described via a set of alternatives A = {A1, A2, …, Am}, a set of criteria C = {C1, C2, …, Cn} such that C is divided into N partitions Ck = {C1, C2, …, C|Ck|} (k = 1, 2, …, N) and there are interrelationships among different criteria in each Ck whereas the criteria in different Ck are independent of each other, a vector of weights of criteria w = [w1, w2, …, wn] such that 0 ≤ w1, w2, …, wn ≤ 1, w1+w2+…+wn = 1, and each element respectively denotes the relative importance of C1, C2, …, Cn, a set of experts E = {E1, E2, …, EM}, a vector of weights of experts ϖ = [ϖ1, ϖ2, …, ϖM] such that 0 ≤ ϖ1, ϖ2, …, ϖM ≤ 1, ϖ1+ϖ2+…+ϖM = 1, and each element respectively denotes the relative importance of E1, E2, …, EM, and M q-rung orthopair fuzzy decision matrices Mh = [Qh,i,j]m×n (h = 1, 2,…, M; i = 1, 2,…, m; j = 1, 2,…, n) such that Qh,i,j = <μh,i,j, νh,i,j> is a qROFN which denotes the evaluation value of Cj with respect to Ai provided by Ek. Based on these components, the MCGDM problem can be described as: Determining the optimal alternative with the help of a ranking of the elements of A based on C, Mh, w, and ϖ. Using the qROFWAPPMM operator, the problem is solved according to the following steps:

  1. Normalise the q-rung orthopair fuzzy decision matrices Mh. Generally, a MCGDM problem may contain two types of criteria, i.e. benefit and cost criteria, which affect the results of decision making positively and negatively, respectively. To eliminate the negative effect, the q-rung orthopair fuzzy decision matrices Mh = [Qh,i,j]m×n = <μh,i,j, νh,i,j> are normalised as (51)
  2. Calculate the power weights of Qh,i,j. The power weights of Qh,i,j are computed using (52) where D(Qh,i,j, Qx,i,j) (D(Qz,i,j, Qy,i,j)) is the Minkowski-type distance between Qh,i,j and Qx,i,j (Qz,i,j and Qy,i,j) that can be computed via the Equation in Definition 5.
  3. Calculate the collective values of Qh,i,j. Taking the normalised q-rung orthopair fuzzy decision matrices MN,h and the expert weight vector ϖ as input, the collective values of Qh,i,j are computed using (53) where qROFWAPPMM is an arbitrary specific qROFWAPPMM operator, such as the qROFWAAPPMM operator in Eq (35), the qROFWAEPPMM operator in Eq (36), the qROFWAHPPMM operator in Eq (39), and the qROFWAFPPMM operator in Eq (42), and the values of the elements in Δ = (δ1, δ2, …, δM) (It is worth nothing that there is only one partition) are determined via identifying the interrelationships among different experts. In general, all experts should be mutually independent. Thus δ1 > 0 and δ2 = δ3 = … = δM = 0.
  4. Calculate the power weights of Qi,j. The power weights of Qi,j are computed using (54) where D(Qi,j, Qi,r) (D(Qi,r, Qi,s)) is the Minkowski-type distance between Qi,j and Qi,r (Qi,t and Qi,s) that can be computed via the Equation in Definition 5.
  5. Calculate the collective values of Qi,j. On the basis of the N partitions Ck = {C1, C2, …, C|Ck|}, suppose Si = {Qi,1, Qi,2, …, Qi,n} is an ordered set of Qi,1, Qi,2, …, Qi,n, Si,k = {Qi,1, Qi,2, …, Qi,|Sk|} are N partitions of S that correspond to Ck, δ1, δ2, …, δ|Si,k| (k = 1, 2, …, N and δ1, δ2, …, δ|Si,k| ≥ 0 but not at the same time δ1 = δ2 = … = δ|Si,k| = 0) are |Si,k| real numbers that respectively correspond to Qi,1, Qi,2, …, Qi,|Sk|, Δk = (δ1, δ2, …, δ|Si,k|) is a collection of δ1, δ2, …, δ|Si,k|, Δ = (Δ1, Δ2, …, ΔN) is a collection of Δ1, Δ2, …, ΔN, p(ik) is a permutation of (1, 2,…, |Si,k|), and P|Si,k| is a set of all permutations of (1, 2,…, |Si,k|). Then the collective values of Qi,j are computed using (55) where qROFWAPPMM is the same specific qROFWAPPMM operator used in Eq (53), and the values of the elements in Δ = (Δ1, Δ2, …, ΔN) (Δk = (δ1, δ2, …, δ|Si,k|)) are determined via identifying the interrelationships among different criteria in Ck. When all criteria in Ck are mutually independent, then δ1 > 0 and δ2 = δ3 = … = δ|Si,k| = 0; When there are interrelationships between any two criteria in Ck, δ1, δ2 > 0 and δ3 = δ4 = … = δ|Si,k| = 0; When there are interrelationships among any d (d = 3, 4, …, n) criteria in Ck, δ1, δ2, …, δd > 0 and δd+1 = δd+2 = … = δ|Si,k| = 0.
  6. Calculate the scores and accuracies of Qi. The scores and accuracies of Qi can be respectively computed using the Equations in Definitions 2 and 3.
  7. Generate a ranking of Ai. On the basis of the scores and accuracies of Qi, a ranking of Ai can be generated according to the comparison rules in Definition 4.
  8. Determine the optimal alternative. The optimal alternative is determined with the help of the ranking.

5. Example, experiments, and comparisons

In this section, a numerical example is firstly used to illustrate the working process of the proposed MCGDM method. Then a set of experiments are carried out to explore the influence of different specific operators and parameter values on the aggregation results. Finally, qualitative and quantitative comparisons to the existing methods are made to demonstrate the feasibility and effectiveness of the proposed method.

5.1. Example

A numerical example about the determination of the best industry for investment from five possible industries (adapted on the basis of Reference [31]) is used to demonstrate the proposed MCGDM method.

To make full use of idle capital, the board of directors of a company decided to invest in a new industry. Five industries were identified as possible industries for investment after preliminary research. The five alternative industries are medical industry (A1), real estate development industry (A2), Internet industry (A3), education and training industry (A4), and manufacturing industry (A5). To select the best industry for investment, the board of directors appointed an expert panel, which consists of four different experts E1, E2, E3, and E4. The relative importance of these experts is quantified by the weight vector ϖ = [0.30, 0.22, 0.28, 0.20]. The four experts were asked to evaluate the five alternative industries on the basis of five criteria, which are the amount of capital profit (C1), the market potential (C2), the risk of capital loss (C3), the growth potential (C4), and the stability of policy (C5). The relative importance of these criteria is measured by the weight vector w = [0.20, 0.20, 0.15, 0.25, 0.20]. According to the structure of interrelationships, the five criteria are divided into two partitions C1 = {C1, C3, C5} and C2 = {C2, C4} and there are interrelationships among the three criteria in C1, so do the two criteria in C2, and C1 and C2 are independent of each other. To provide enough freedom in the evaluation of the values of the five criteria of each alternative industry, experts were allowed to use qROFNs. The evaluation results of the four experts are respectively listed in the following four matrices:

With the known conditions above, the determination of the best industry for investment can be carried out leveraging the proposed MCGDM method. Its process consists of the following eight steps:

  1. Normalise the q-rung orthopair fuzzy decision matrices Mh (h = 1, 2, 3, 4). Among the five criteria, the amount of capital profit (C1), the market potential (C2), the growth potential (C4), and the stability of policy (C5) are four benefit criteria and the risk of capital loss (C3) is a cost criterion. Based on this, the four matrices Mh are normalised according to Eq (51) and four normalised matrices are obtained as follows:
  2. Calculate the power weights of Qh,i,j. According to Eq (52), the power weights of Qh,i,j are computed and the computed results are listed in the following matrices:
  3. Calculate the collective values of Qh,i,j. Taking the normalised q-rung orthopair fuzzy decision matrices MN,h and the power weight matrices [Wh,i,j]5×5 as input, the collective values of Qh,i,j can be computed according to Eq (53). Without loss of generality (It is worth nothing that the qROFWAHPPMM operator (see Eq (39)) will reduce to the qROFWAAPPMM operator (see Eq (35)) when λ = 1 and will reduce to the qROFWAEPPMM operator (see Eq (36)) when λ = 2), the qROFWAHPPMM operator in Eq (39) (When adapting this operator, q = 3, λ = 3, and Δ = (δ1, δ2, δ3, δ4) = (1, 0, 0, 0)) is used in Eq (53) to complete the computation. The computed results are listed in the following matrix:
  4. Calculate the power weights of Qi,j. According to Eq (54), the power weights of Qi,j are computed and the computed results are listed in the following matrix:
  5. Calculate the collective values of Qi,j. On the basis of the two partitions C1 = {C1, C3, C5} and C2 = {C2, C4}, the matrix MN is divided into Si,1 = {Qi,1, Qi,3, Qi,5} and Si,2 = {Qi,2, Qi,4} (i = 1, 2, 3, 4, 5). Taking them and the power weight matrix [Wi,j]5×5 as input, the collective values of Qi,j can be computed according to Eq (55). Since the qROFWAHPPMM operator in Eq (39) has been used in the third step, this operator (When adapting the operator, q = 3, λ = 3, and Δ = (Δ1, Δ2) = ((δ1, δ2, δ3), (δ1, δ2)) = ((1, 2, 3), (1, 2))) is also leveraged in this step to complete the computation. The computed results are listed as follows:
    Q1 = <0.5453, 0.5668>, Q2 = <0.5159, 0.7220>, Q3 = <0.6024, 0.6563>, Q4 = <0.5303, 0.6991>, Q5 = <0.5278, 0.6953>
  6. Calculate the scores and accuracies of Qi. The scores and accuracies of Qi are respectively computed using the Equations in Definitions 2 and 3. The computed results are respectively listed as follows:
    S(Q1) = −0.0199, S(Q2) = −0.2390, S(Q3) = −0.0641, S(Q4) = −0.1925, S(Q5) = −0.1891
    A(Q1) = 0.3443, A(Q2) = 0.5137, A(Q3) = 0.5013, A(Q4) = 0.4908, A(Q5) = 0.4832
  7. Generate a ranking of Ai. On the basis of the S(Qi) and A(Qi), a ranking of Ai is generated according to the comparison rules in Definition 4. The ranking is as follow: A1A3A5A4A2
  8. Determine the optimal alternative. The best industry for investment is determined as medical industry (A1).

5.2. Experiments

To explore the effect of using different specific operators and assigning different parameter values on the aggregation results, the following three experiments were carried out:

(1) Experiment 1 aims to show the influence of using different specific operators on the aggregation results. In this experiment, the presented qROFWAAPPMM (see Eq (35)), qROFWAEPPMM (see Eq (36)), qROFWAHPPMM (see Eq (39)), and qROFWAFPPMM (see Eq (42)) operators are respectively used to calculate the collective values of Qh,i,j and the collective values of Qi,j in the numerical example (When calculating the collective values of Qh,i,j, q = 3, λ = 3, ε = 2, and Δ = (δ1, δ2, δ3, δ4) = (1, 0, 0, 0); When calculating the collective values of Qi,j, q = 3, λ = 3, ε = 2, and Δ = (Δ1, Δ2) = ((δ1, δ2, δ3), (δ1, δ2)) = ((1, 2, 3), (1, 2))). The results of the experiment are the calculated scores of Qi and the generated rankings of Ai, which are listed in Table 1. As can be seen from the table, there is slight difference among the scores of the same Qi calculated by the four pairs of specific operators, and the rankings of Ai also indicate small difference with respect to the four pairs of specific operators. These indicate that the use of different specific operators has no obvious influence on the aggregation results.

(2) Experiment 2 aims to show the influence of assigning different q values on the aggregation results. In this experiment, the presented qROFWAHPPMM (see Eq (39) and qROFWAFPPMM (see Eq (42)) operators (Since the qROFWAHPPMM operator is the generalisation of the qROFWAAPPMM and qROFWAEPPMM operators, they are not included in this experiment and will not be included in the subsequent experiments and comparisons for the sake of simplicity) are respectively used to calculate the collective values of Qh,i,j and the collective values of Qi,j in the numerical example (When calculating the collective values of Qh,i,j, λ = 3, ε = 2, and Δ = (δ1, δ2, δ3, δ4) = (1, 0, 0, 0); When calculating the collective values of Qi,j, λ = 3, ε = 2, and Δ = (Δ1, Δ2) = ((δ1, δ2, δ3), (δ1, δ2)) = ((1, 2, 3), (1, 2))). The results of the experiment are the calculated scores of Qi and the generated rankings of Ai, which are depicted in Fig 1. From the figure, it can be seen that the ranking will change as q changes. For the pair of qROFWAHPPMMs, the ranking is A1A3A5A4A2 when q = 3, changes to A1A3A4A5A2 when q = 4, 5, 6, and becomes A3A1A4A5A2 when q = 7, 8, 9, 10; The best alternative changes from A1 to A3 from q = 7. For the pair of qROFWAFPPMMs, the ranking is A1A3A4A5A2 when q = 3, changes to A3A1A4A5A2 when q = 4, 5, 6, 7, 8, 9, and becomes A3A4A1A5A2 when q = 10; The best alternative changes from A1 to A3 from q = 4. Based on these results, it is recommended that the smallest q which can meet 0 ≤ μq + νq ≤ 1 is assigned in practical applications. For example, since MN,3 contains <0.9, 0.6>, q is assigned 3 since 0.92 + 0.62 > 1 and 0.93 + 0.63 < 1.

(3) Experiment 3 aims to show the influence of assigning different λ (ε) values on the aggregation results. In this experiment, the presented qROFWAHPPMM (see Eq (39) and qROFWAFPPMM (see Eq (42)) operators are respectively used to calculate the collective values of Qh,i,j and the collective values of Qi,j in the numerical example (When calculating the collective values of Qh,i,j, q = 3 and Δ = (δ1, δ2, δ3, δ4) = (1, 0, 0, 0); When calculating the collective values of Qi,j, q = 3 and Δ = (Δ1, Δ2) = ((δ1, δ2, δ3), (δ1, δ2)) = ((1, 2, 3), (1, 2))). The results of the experiment are the calculated scores of Qi and the generated rankings of Ai, which are depicted in Fig 2. It can be seen from the figure that the scores computed by the pair of qROFWAHPPMMs (qROFWAFPPMMs) gradually decrease as λ (ε) gradually increases. Therefore, the parameter λ (ε) can be seen as a pessimistic factor for MCGDM problems. Generally, if the attitude of a decision maker is neutral, a small λ (ε) (e.g. λ = 1, 2, 3; ε = 2, 3, 4) is recommended. If the attitude is pessimistic enough, a bigger λ (ε) can be assigned when the pair of qROFWAHPPMMs (qROFWAFPPMMs) is used. Otherwise, a smaller λ (ε) is recommended.

5.3. Comparisons

As mentioned in the introduction, more than twenty different aggregation operators of qROFNs have been presented within academia. Representative examples are the WA and WG [26], WBM and WGBM [27], WABM [28], WPBM and WPGBM [29], WHM and WGHM [30], WHM* and WPHM [31], WMSM and WGMSM [32], WPMSM [33], WPPMSM [34], WMM and WGMM [35], WEBM [36], WE [37], and WP [38] operators. In this subsection, qualitative and quantitative comparisons between the MCGDM methods based on these operators and the proposed MCGDM method are carried out to demonstrate its feasibility and effectiveness.

5.3.1. Qualitative comparison.

Generally, a qualitative comparison among different MCGDM methods can be carried out by comparing their characteristics. For the twenty existing methods and the proposed method, the generality and flexibility in the aggregation of qROFNs, the capability to deal with the interrelationships among different criteria, and the capability to reduce the negative influence of the unduly high or unduly low criterion values on the aggregation results are selected as the comparison characteristics. The results of the comparison are shown in Table 2. The details of the comparison are explained as follows:

thumbnail
Table 2. The results of the qualitative comparison.

Note: Heterogeneous interrelationships refer to the situation in which the criteria are divided into several parts and there are interrelationships among different criteria in each part whereas the criteria in different parts are independent of each other.

https://doi.org/10.1371/journal.pone.0221759.t002

  1. Generality and flexibility: For the WP method, any one of the twenty different WP operators can be used in the aggregation. Therefore, its generality and flexibility can be seen as moderate. The generality and flexibility of the WABM method and the proposed (WAPPMM) method are desirable since the aggregations are based on the operations of any family of ATTs. The aggregations in the remaining methods are based on the operation of a specific family of ATT. Relatively, they have limited generality and flexibility.
  2. When all criteria are independent of each other: It is no doubt that all of the listed methods can deal with this case.
  3. When there are interrelationships between any two criteria: The WE, WP, WA, and WG methods are only suitable for the independent case. All other methods have the capability of dealing with the case in which there are interrelationships between any two criteria.
  4. When there are interrelationships among any multiple criteria: The WMSM, WGMSM, WPMSM, WPPMSM, WMM, and WGMM methods and the proposed method can handle this situation because of the use of the MSM or MM operator.
  5. When there are heterogeneous interrelationships among different criteria: The WPBM, WPGBM, WPHM, WPPMSM, and WEBM methods and the proposed method can deal with this case due to the combination of the partitioned average operator.
  6. Capability to reduce the negative effect: The WPMSM and WPPMSM methods and the proposed method have this capability because of the combination of the PA operator.

As can be summarised from the qualitative comparison above, the proposed method has desirable generality and flexibility at both aggregating the q-rung orthopair fuzzy information and dealing with the interrelationships of criteria, and has the capability to reduce the negative influence caused by the deviation of some criterion values.

5.3.2. Quantitative comparison.

In general, a quantitative comparison among different MCGDM methods can be carried out using the same numerical example. Here the numerical example in subsection 5.1 is used to quantitatively compare the proposed method to the WA, WG, WBM, WGBM, WABM, WPBM, WPGBM, WHM, WGHM, WHM*, WPHM, WMSM, WGMSM, WPMSM, WPPMSM, WMM, and WGMM methods (Please note that the WEBM, WE, and WP operators were not included in the quantitative comparison because the WEBM method was pre-sented to handle the situation with unknown criterion weights, the WE method was presented to deal with the case where the values of criteria are expressed by a fuzzy number in Zadeh’s FS and the values of weights are described by qROFNs, the WP method was presented to control the uncertainty of the evaluation results of some experts, and the proposed method does not consider these characteristics). In this comparison, the qROFNs Qh,1,4 and Qh,1,5 (h = 1, 2, 3, 4) in Mh were constantly adjusted from high score to low score according to Table 3. It is easy to guess that such adjustments will lower the ranking of A1, because the score becomes lower and lower. To confirm this conjecture, the methods were implemented to generate the changes of the places of A1 in their respective rankings. The specific generation process is as follows:

  1. Calculate the collective values of Qh,i,j. A specific operator presented in each method, as listed in Table 4, was leveraged to compute the collective values of Qh,i,j in the numerical example on the basis of certain Qh,1,4 and Qh,1,5 and certain parameters (whose values are also listed in Table 4).
  2. Calculate the collective values of Qi,j. For each method, the same specific operator used in the calculation of the collective values of Qh,i,j, as listed in Table 4, was leveraged to compute the collective values of Qi,j on the basis of the computed collective values of Qh,i,j of each group of Qh,1,4 and Qh,1,5 and certain parameters (whose values are also listed in Table 4).
  3. Calculate the scores of Qi. For each method, the scores of Qi of each group of Qh,1,4 and Qh,1,5 were computed according to the calculated collective values of Qi,j of this group of Qh,1,4 and Qh,1,5 and the Equation in Definition 2.
  4. Generate the changes of the places of A1. According to the calculated scores of Qi of each group of Qh,1,4 and Qh,1,5 of each method (The calculated scores of Q1 of the nine groups of Qh,1,4 and Qh,1,5 of each method are listed in Table 5), the changes of the places of A1 in the rankings of all comparison methods were generated and are depicted in Figs 312.
thumbnail
Fig 3. The changes of the places of A1 in the rankings of the WA and WG methods.

https://doi.org/10.1371/journal.pone.0221759.g003

thumbnail
Fig 4. The changes of the places of A1 in the rankings of the WBM and WGBM methods.

https://doi.org/10.1371/journal.pone.0221759.g004

thumbnail
Fig 5. The changes of the places of A1 in the rankings of the WABM method.

https://doi.org/10.1371/journal.pone.0221759.g005

thumbnail
Fig 6. The changes of the places of A1 in the rankings of the WPBM and WPGBM methods.

https://doi.org/10.1371/journal.pone.0221759.g006

thumbnail
Fig 7. The changes of the places of A1 in the rankings of the WHM and WGHM methods.

https://doi.org/10.1371/journal.pone.0221759.g007

thumbnail
Fig 8. The changes of the places of A1 in the rankings of the WHM* and WPHM methods.

https://doi.org/10.1371/journal.pone.0221759.g008

thumbnail
Fig 9. The changes of the places of A1 in the rankings of the WMSM and WGMSM methods.

https://doi.org/10.1371/journal.pone.0221759.g009

thumbnail
Fig 10. The changes of the places of A1 in the rankings of the WPMSM and WPPMSM methods.

https://doi.org/10.1371/journal.pone.0221759.g010

thumbnail
Fig 11. The changes of the places of A1 in the rankings of the WMM and WGMM methods.

https://doi.org/10.1371/journal.pone.0221759.g011

thumbnail
Fig 12. The changes of the places of A1 in the rankings of the WAPPMM method.

https://doi.org/10.1371/journal.pone.0221759.g012

thumbnail
Table 3. The values of the adjusted qROFNs in the quantitative comparison.

https://doi.org/10.1371/journal.pone.0221759.t003

thumbnail
Table 4. The specific operators and their parameter values for calculating the collective values of Qh,i,j and Qi,j.

https://doi.org/10.1371/journal.pone.0221759.t004

thumbnail
Table 5. The calculated scores of Q1 of the nine groups of Qh,1,4 and Qh,1,5 of all comparison methods.

https://doi.org/10.1371/journal.pone.0221759.t005

From the ten figures, it can be seen that the results of all comparison methods are consistent with the conjecture. This verifies the feasibility and effectiveness of the proposed method. In addition, as can be found from the figures of WA, WBM, WAHBM, WAFBM, WPBM, WHM, WHM*, WPHM, WMSM, WPMSM, WPPMSM, WMM, WAHPPMM, and WAFPPMM (In addition to these operators, the remaining comparison operators belong to geometric operators. Because the presented operators do not be-long to geometric operators, the geometric operators are not included in the discussion of the quantitative comparison), the rankings of A1 generated by the WA, WAFBM, WMSM, WPMSM, WPPMSM, and WMM methods do not start from the first, which is obviously different from the WBM, WAHBM, WPBM, WHM, WHM*, WPHM, WAHPPMM, and WAFPPMM methods. From the given data in the numerical example it is not difficult to determine that A1 is the best alternative. Therefore, the results of the latter group of methods are more reasonable than that of the former group of methods.

Among the methods in the latter group, the WBM, WAHBM, WHM, and WHM* methods can deal with the case in which all criteria are independent of each other or there are interrelationships between any two criteria, and the WPBM and WPHM methods are applicable for the situation where the criteria are divided into several partitions and there are interrelationships between any two criteria in each partition whereas the criteria in different partitions are independent of each other. These six methods are different with the proposed (WAHPPMM, WAFPPMM) method in characteristics.

Except the six methods, it is also of necessity to compare the results of the WMM method and the proposed (WAHPPMM, WAFPPMM) method. The difference in characteristics between them is that the proposed method can provide desirable generality and flexibility, deal with the heterogeneous relationships among criteria, and reduce the negative effect of the deviation of criterion values. The intuitive manifestation of such difference is that their results are significantly different. As can be seen from Fig 11(A) and Fig 12(A), the ranking generated by the WMM method is A3A1A4A5A2 at 0, changes to A3A4A5A1A2 at 2, and becomes A3A4A5A2A1 at 5; while the ranking generated by the WAHPPMM method is A1A3A5A4A2 at 0, changes to A3A1A5A4A2 at 2 and A3A5A4A1A2 at 5, and becomes A3A5A4A2A1 at 7. Hence, the rankings of the two methods are different, which is caused by the capability to handle the interrelationships among criteria. In addition, the place of A1 in the rankings generated by the WMM method descends faster than that in the rankings of the WAHPPMM method. This is because the WAHPPMM method has the capability to reduce the influence of the distortion of criterion values and the WMM method does not have such capability.

Finally, it should be pointed out that the quantitative comparison does not aim to find out the best method, but to demonstrate the feasibility and effectiveness of the proposed method and illustrate the difference of different methods. Generally, it is difficult to conclude that one MCGDM method is better than the others because each method have its specific characteristics, which determine its specific application scenario. A wider range of methods offer decision makers a greater flexibility when selecting a proper method for their specific application scenario.

6. Conclusion

In this paper, an Archimedean power partitioned MM operator and a weighted Archimedean power partitioned MM operator have been presented to solve the MCGDM problems based on qROFNs. The idempotency and boundedness of the Archimedean power partitioned MM operator have been proved and the four specific expressions of the two operators have been constructed leveraging the operational rules of qROFNs based on the Algebraic, Einstein, Hamacher, and Frank families of ATTs and their additive generators. On the basis of the presented weighted Archimedean power partitioned MM operator, a method for solving the MCGDM problems based on qROFNs has been proposed. The paper has also provided a numerical example coupled with a set of experiments to illustrate the working process of the proposed method and reported qualitative and quantitative comparisons to demonstrate its feasibility and effectiveness. The results of the comparisons suggest that the proposed method is general and flexible at both aggregation of criterion values and capture of criterion interrelationships, and concurrently has the capability to handle the heterogeneous interrelationships of criteria and reduce the negative effect of the biased criterion values.

Future work will focus especially on extending the presented operators from the aspect of dealing with more complex interrelationships of criteria and risk attitudes of decision makers. Further, the application of the proposed method in solving practical decision making problems, such as manufacturing process selection, part build orientation determination, medical diagnosis, and resource evaluation, will also be studied.

Supporting information

S1 File.

Appendixes A, B, and C: The proofs of Theorems 1, 2, and 3.

https://doi.org/10.1371/journal.pone.0221759.s001

(PDF)

Acknowledgments

The authors are very grateful to the five anonymous reviewers for their insightful comments for the improvement of the paper and to Dr. Peizhi Shi at the School of Computer Science, The University of Manchester, UK for his kind help in writing the Java implementation code of the WMM and WGMM operators. The authors also would like to acknowledge the financial supports by the EPSRC UKRI Innovation Fellowship (Ref. EP/S001328/1).

References

  1. 1. Greco S, Figueira J, Ehrgott M. Multiple Criteria Decision Analysis: State of the Art Surveys. New York: Springer-Verlag New York; 2016.
  2. 2. Liu P. Some Hamacher aggregation operators based on the interval-valued intuitionistic fuzzy numbers and their application to group decision making. IEEE Transactions on Fuzzy systems 2013; 22(1): 83–97.
  3. 3. Wan S, Zhu Y. Triangular intuitionistic fuzzy triple Bonferroni harmonic mean operators and application to multi-attribute group decision making. Iranian Journal of Fuzzy Systems 2016; 13(5): 117–145.
  4. 4. Jamkhaneh EB, Garg H. Some new operations over the generalized intuitionistic fuzzy sets and their application to decision-making process. Granular Computing 2018; 3(2): 111–122.
  5. 5. Liu P, Tang G. Some intuitionistic fuzzy prioritized interactive Einstein Choquet operators and their application in decision making. IEEE Access 2018; 6: 72357–72371.
  6. 6. Dong J, Chen Y, Wan S. A cosine similarity based QUALIFLEX approach with hesitant fuzzy linguistic term sets for financial performance evaluation. Applied Soft Computing 2018; 69: 316–329.
  7. 7. Wan S, Jin Z, Dong J. Pythagorean fuzzy mathematical programming method for multi-attribute group decision making with Pythagorean fuzzy truth degrees. Knowledge and Information Systems 2018; 55(2): 437–466.
  8. 8. Wang J, Gao H, Wei G, Wei Y. Methods for multiple-attribute group decision making with q-rung interval-valued orthopair fuzzy information and their applications to the selection of green suppliers. Symmetry 2019; 11(1): 56.
  9. 9. Yang W, Wang C, Liu Y, Sun Y. Hesitant Pythagorean fuzzy interaction aggregation operators and their application in multiple attribute decision-making. Complex & Intelligent Systems 2019; 5(2): 199–216.
  10. 10. Garg H, Kaur G. TOPSIS based on nonlinear-programming methodology for solving decision-making problems under cubic intuitionistic fuzzy set environment. Computational and Applied Mathematics 2019; 38(3): 114.
  11. 11. Garg H, Kumar K. Linguistic interval-valued Atanassov intuitionistic fuzzy sets and their applications to group decision-making problems. IEEE Transactions on Fuzzy Systems 2019;
  12. 12. Wan S, Jin Z, Dong J. A new order relation for Pythagorean fuzzy numbers and application to multi-attribute group decision making. Knowledge and Information Systems 2019; https://doi.org/10.1007/s10115-019-01369-8.
  13. 13. Bustince H, Barrenechea E, Pagola M, Fernandez J, Xu Z, Bedregal B, Montero J, Hagras H, Herrera F, De Baets B. A Historical Account of Types of Fuzzy Sets and Their Relationships. IEEE Transactions on Fuzzy Systems 2016; 24(1): 179–194.
  14. 14. Yager RR. Generalized Orthopair Fuzzy Sets. IEEE Transactions on Fuzzy Systems 2017; 25(5): 1222–1230.
  15. 15. Zadeh LA. Fuzzy Sets. Information and Control 1965; 8(3): 338–353.
  16. 16. Atanassov KT. Intuitionistic Fuzzy Sets. Fuzzy Sets and Systems 1986; 20(1): 87–96.
  17. 17. Yager RR. Pythagorean Membership Grades in Multicriteria Decision Making. IEEE Transactions on Fuzzy Systems 2014; 22(4): 958–965.
  18. 18. Du WS. Correlation and correlation coefficient of generalized orthopair fuzzy sets. International Journal of Intelligent Systems 2019; 34(4): 564–583.
  19. 19. Du WS. Minkowski-type distance measures for generalized orthopair fuzzy sets. International Journal of Intelligent Systems 2018; 33(4): 802–817.
  20. 20. Wang P, Wang J, Wei G, Wei C. Similarity measures of q-rung orthopair fuzzy sets based on cosine function and their applications. Mathematics 2019; 7(4): 340.
  21. 21. Wang R, Li Y. A Novel Approach for Green Supplier Selection under a q-Rung Orthopair Fuzzy Environment. Symmetry, 2018; 10(12): 687.
  22. 22. Peng X, Dai J. Research on the assessment of classroom teaching quality with q-rung orthopair fuzzy information based on multiparametric similarity measure and combinative distance-based assessment. International Journal of Intelligent Systems 2019; https://doi.org/10.1002/int.22109.
  23. 23. Du WS. Research on arithmetic operations over generalized orthopair fuzzy sets. International Journal of Intelligent Systems 2019; 34(5): 709–732.
  24. 24. Gao J, Liang Z, Shang J, Xu Z. Continuities, Derivatives and Differentials of q-Rung Orthopair Fuzzy Functions, IEEE Transactions on Fuzzy Systems 2018;
  25. 25. Ye J, Ai Z, Xu Z. Single variable differential calculus under q-rung orthopair fuzzy environment: Limit, derivative, chain rules, and its application. International Journal of Intelligent Systems 2019; https://doi.org/10.1002/int.22100.
  26. 26. Liu P, Wang P. Some q-rung orthopair fuzzy aggregation operators and their applications to multiple-attribute decision making. International Journal of Intelligent Systems 2018; 33(2): 259–280.
  27. 27. Liu P, Liu J. Some q-rung orthopair fuzzy Bonferroni mean operators and their application to multi-attribute group decision making. International Journal of Intelligent Systems 2018; 33(2): 315–347.
  28. 28. Liu P, Wang P. Multiple-Attribute Decision Making based on Archimedean Bonferroni Operators of q-Rung Orthopair Fuzzy Numbers. IEEE Transactions on Fuzzy Systems 2019; 27(5): 834–848.
  29. 29. Yang W, Pang Y. New q-rung orthopair fuzzy partitioned Bonferroni mean operators and their application in multiple attribute decision making. International Journal of Intelligent Systems 2019; 34(3): 439–476.
  30. 30. Wei G, Gao H, Wei Y. Some q-rung orthopair fuzzy Heronian mean operators in multiple attribute decision making. International Journal of Intelligent Systems 2018; 33(7): 1426–1458.
  31. 31. Liu Z, Wang S, Liu P. Multiple attribute group decision making based on q-rung orthopair fuzzy Heronian mean operators. International Journal of Intelligent Systems 2018; 33(12): 2341–2363.
  32. 32. Wei G, Wei C, Wang J, Gao H, Wei Y. Some q-rung orthopair fuzzy Maclaurin symmetric mean operators and their applications to potential evaluation of emerging technology commercialization. International Journal of Intelligent Systems 2019; 34(1): 50–81.
  33. 33. Liu P, Chen SM, Wang P. Multiple-Attribute Group Decision-Making Based on q-Rung Orthopair Fuzzy Power Maclaurin Symmetric Mean Operators. IEEE Transactions on Systems, Man, and Cybernetics: Systems 2018;
  34. 34. Bai K, Zhu X, Wang J, Zhang R. Some partitioned Maclaurin symmetric mean based on q-rung orthopair fuzzy information for dealing with multi-attribute group decision making. Symmetry 2018; 10(9): 383.
  35. 35. Wang J, Zhang R, Zhu X, Zhou Z, Shang X, Li W. Some q-rung orthopair fuzzy Muirhead means with their application to multi- attribute group decision making. Journal of Intelligent & Fuzzy Systems 2019; 36(2): 1599–1614.
  36. 36. Liu Z, Liu P, Liang X. Multiple attribute decision-making method for dealing with heterogeneous relationship among attributes and unknown attribute weight information under q-rung orthopair fuzzy environment. International Journal of Intelligent Systems 2018; 33(9): 1900–1928.
  37. 37. Peng X, Dai J, Garg H. Exponential operation and aggregation operator for q-rung orthopair fuzzy set and their decision-making method with a new score function. International Journal of Intelligent Systems 2018; 33(11): 2255–2282.
  38. 38. Xing Y, Zhang R, Zhou Z, Wang J. Some q-rung orthopair fuzzy point weighted aggregation operators for multi-attribute decision making. Soft Computing 2019; https://doi.org/10.1007/s00500-018-03712-7.
  39. 39. Muirhead RF. Some methods applicable to identities and inequalities of symmetric algebraic functions of n letters. Proceedings of the Edinburgh Mathematical Society 1902; 21: 144–162.
  40. 40. Garg H, Arora R. Bonferroni mean aggregation operators under intuitionistic fuzzy soft set environment and their applications to decision-making. Journal of the Operational Research Society 2018; 69(11): 1711–1724.
  41. 41. Kaur G, Garg H. Multi-attribute decision-making based on Bonferroni mean operators under cubic intuitionistic fuzzy set environment. Entropy 2018; 20(1): 65.
  42. 42. Yang W, Pang Y. New Pythagorean fuzzy interaction Maclaurin symmetric mean operators and their application in multiple attribute decision making. IEEE Access 2018; 6: 39241–39260.
  43. 43. Wei G, Garg H, Gao H, Wei C. Interval-valued Pythagorean fuzzy Maclaurin symmetric mean operators in multiple attribute decision making. IEEE Access 2018; 6: 67866–67884.
  44. 44. Liu P, Li D. Some Muirhead mean operators for intuitionistic fuzzy numbers and their applications to group decision making. PLOS ONE. 2017; 12(1): e0168767. pmid:28103244
  45. 45. Garg H. Multi-criteria decision-making method based on prioritized Muirhead mean aggregation operator under neutrosophic set environment. Symmetry 2018; 10(7): 280.
  46. 46. Li L, Zhang R, Wang J, Zhu X, Xing Y. Pythagorean fuzzy power Muirhead mean operators with their application to multi-attribute decision making. Journal of Intelligent & Fuzzy Systems 2018; 35(2): 2035–2050.
  47. 47. Liu P, Teng F. Some Muirhead mean operators for probabilistic linguistic term sets and their applications to multiple attribute decision-making. Applied Soft Computing. 2018; 68: 396–431.
  48. 48. Xu Y, Shang X, Wang J, Zhao H, Zhang R, Bai K. Some Interval-Valued q-Rung Dual Hesitant Fuzzy Muirhead Mean Operators with Their Application to Multi-Attribute Decision-Making. IEEE Access 2019; 7: 54724–54745.
  49. 49. Wang J, Wei G, Wei C, Wei Y. Dual Hesitant q-Rung Orthopair Fuzzy Muirhead Mean Operators in Multiple Attribute Decision Making. IEEE Access 2019; 67139–67166.
  50. 50. Liu P, Chen SM, Liu J. Multiple attribute group decision making based on intuitionistic fuzzy interaction partitioned Bonferroni mean operators. Information Sciences 2017; 411: 98–121.
  51. 51. Yang W, Shi J, Liu Y, Pang Y, Lin R. Pythagorean fuzzy interaction partitioned Bonferroni mean operators and their application in multiple-attribute decision-making. Complexity 2018; 2018: 3606245.
  52. 52. Liu P, Liu J, Merigó JM. Partitioned Heronian means based on linguistic intuitionistic fuzzy numbers for dealing with multi-attribute group decision making. Applied Soft Computing 2018; 62: 395–422.
  53. 53. Liu P, Liu J. Partitioned Bonferroni mean based on two-dimensional uncertain linguistic variables for multiattribute group decision making. International Journal of Intelligent Systems 2019; 34(2): 155–187.
  54. 54. Liu P, Gao H, Ma J. Novel green supplier selection method by combining quality function deployment with partitioned Bonferroni mean operator in interval type-2 fuzzy environment. Information Sciences 2019; 490: 292–316.
  55. 55. Zhu X, Bai K, Wang J, Zhang R, Xing Y. Pythagorean fuzzy interaction power partitioned Bonferroni means with applications to multi-attribute group decision making. Journal of Intelligent & Fuzzy Systems 2019; 36(4): 3423–3438.
  56. 56. Yager RR. The Power Average Operator. IEEE Transactions on Systems, Man, and Cybernetics-Part A: Systems and Humans. 2001; 31(6): 724–731.
  57. 57. Wan S, Dong J. Power geometric operators of trapezoidal intuitionistic fuzzy numbers and application to multi-attribute group decision making. Applied Soft Computing 2015; 29: 153–168.
  58. 58. Wan S, Yi Z. Power average of trapezoidal intuitionistic fuzzy numbers using strict t-norms and t-conorms. IEEE Transactions on Fuzzy Systems 2015; 24(5): 1035–1047.
  59. 59. Wei G, Lu M. Pythagorean fuzzy power aggregation operators in multiple attribute decision making. International Journal of Intelligent Systems 2018; 33(1): 169–186.
  60. 60. Rani D, Garg H. Complex intuitionistic fuzzy power aggregation operators and their applications in multicriteria decision-making. Expert Systems 2018; 35(6): e12325.
  61. 61. Liu P, Khan Q, Mahmood T, Hassan N. T-Spherical Fuzzy Power Muirhead Mean Operator Based on Novel Operational Laws and Their Application in Multi-Attribute Group Decision Making. IEEE Access 2019; 7: 22613–22632.
  62. 62. Garg H, Arora R. Generalized intuitionistic fuzzy soft power aggregation operator based on t-norm and their application in multicriteria decision-making. International Journal of Intelligent Systems 2019; 34(2): 215–246.
  63. 63. Garg H. Generalized intuitionistic fuzzy interactive geometric interaction operators using Einstein t-norm and t-conorm and their application to decision making. Computers & Industrial Engineering 2016; 101: 53–69.
  64. 64. Lei Q, Xu Z, Bustince H, Fernandez J. Intuitionistic fuzzy integrals based on Archimedean t-conorms and t-norms. Information Sciences. 2016; 327: 57–70.
  65. 65. Liu P. The aggregation operators based on Archimedean t-conorm and t-norm for single-valued neutrosophic numbers and their application to decision making. International Journal of Fuzzy Systems 2016; 18(5): 849–863.
  66. 66. Zhang N, Yao Z, Zhou Y, Wei G. Some new dual hesitant fuzzy linguistic operators based on Archimedean t-norm and t-conorm. Neural Computing and Applications 2018; https://doi.org/10.1007/s00521-018-3534-x.
  67. 67. Kaur G, Garg H. Generalized cubic intuitionistic fuzzy aggregation operators using t-norm operations and their applications to group decision-making process. Arabian Journal for Science and Engineering 2019; 44(3): 2775–2794.
  68. 68. Garg H, Rani D. New generalised Bonferroni mean aggregation operators of complex intuitionistic fuzzy information based on Archimedean t-norm and t-conorm. Journal of Experimental & Theoretical Artificial Intelligence 2019; https://doi.org/10.1080/0952813X.2019.1620871.