Next Article in Journal
Comparing the Normalized Difference Vegetation Index with the Google Street View Measure of Vegetation to Assess Associations between Greenness, Walkability, Recreational Physical Activity, and Health in Ottawa, Canada
Next Article in Special Issue
The Oakville Oil Refinery Closure and Its Influence on Local Hospitalizations: A Natural Experiment on Sulfur Dioxide
Previous Article in Journal
Effect of Bioaugmentation on Biogas Yields and Kinetics in Anaerobic Digestion of Sewage Sludge
Previous Article in Special Issue
Predicting Infectious Disease Using Deep Learning and Big Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel Multiattribute Decision-Making Method Based on Point–Choquet Aggregation Operators and Its Application in Supporting the Hierarchical Medical Treatment System in China

1
School of Economics and Management, Beijing Jiaotong University, Beijing 100044, China
2
School of Mechanical, Electronic and Control Engineering, Beijing Jiaotong University, Beijing 100044, China
*
Authors to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2018, 15(8), 1718; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15081718
Submission received: 29 June 2018 / Revised: 27 July 2018 / Accepted: 27 July 2018 / Published: 10 August 2018

Abstract

:
The hierarchical medical treatment system is an efficient way to solve the problem of insufficient and unbalanced medical resources in China. Essentially, classifying the different degrees of diseases according to the doctor’s diagnosis is a key step in pushing forward the hierarchical medical treatment system. This paper proposes a framework to solve the problem where diagnosis values are given as picture fuzzy numbers (PFNs). Point operators can reduce the uncertainty of doctor’s diagnosis and get intensive information in the process of decision making, and the Choquet integral operator can consider correlations among symptoms. In order to take full advantage of these two kinds of operators, in this paper, we firstly define some point operators under the picture fuzzy environment, and further propose a new class of picture fuzzy point–Choquet integral aggregation operators. Moreover, some desirable properties of these operators are also investigated in detail. Then, a novel approach based on these operators for multiattribute decision-making problems in the picture fuzzy context is introduced. Finally, we give an example to illustrate the applicability of the new approach in assisting hierarchical medical treatment system. This is of great significance for integrating the medical resources of the whole society and improving the service efficiency of the medical service system.

1. Introduction

With increasing environmental issues, lung diseases are becoming a serious health problem in China. As the medical facilities in grade III, class A hospitals are much better than those of in other small hospitals, people prefer to go to those relatively high-level hospitals for treatment. As a result, overcrowding in large hospitals is common, far exceeding the coping capacity. At the same time, however, small hospitals or clinics waste medical resources. Under such circumstances, how to better allocate limited medical resources and improve the input and output efficiency of the health care system are new challenges for the medical system in China.
Developing a hierarchical medical treatment system is regarded as key an effective way to solve the problem of insufficient and unbalanced medical resources, in which medical institutions at various levels receive patients according to the degree and urgency of the diseases they have. In such a system, common illnesses are treated at basic clinics, with patients transferred to more specialized facilities if their condition demands it. Serious illnesses should be treated in higher-grade hospitals. At the same time, higher-grade hospitals can also transfer patients down to lower-grade ones as their condition stabilizes. Thus, determining the severity of the illness is a key action in this system. At present, with the increase in the number of patients with lung diseases, establishing an appropriate approach to divide patients under different conditions into different levels of hospitals is an effective way to make full use of limited medical resources and cure more patients with lung diseases. However, the diagnosis for patient’s condition often involves multiple correlative criteria and thus can be described as multiattribute decision making (MADM) problems. This paper proposes a general framework in order to solve the MADM problem, which can be applied in the above scenario.
The essence of MADM is the process of ranking the alternatives and selecting an optimal scheme among a set of alternatives with respect to a list of attribute value. Recently, MADM has received much attention from scholars and has been widely applied to economic management and daily life. For example, Tang et al. [1] proposed an algorithm for group decision making with incomplete hesitant fuzzy linguistic preference relations and applied it to flood disaster risk evaluation. Qi [2] developed two effective multicriteria decision making (MCDM) approaches based on defined prioritized average aggregation operators and applied them to tackle complex emergency response solutions evaluation problems. Lin [3] proposed a linear program and a procedure for solving linguistic MADM problems with risk preferences and incomplete weight information, and further applied it to low-carbon tourism destination selection. Due to the increased complexity of real decision-making problems, we usually have to face the difficulty of representing attribute values appropriately. Chatterjee et al. [4] proposed a novel hybrid method encompassing factor relationship and multi-attributive border approximation area comparison methods for selection and evaluation of non-traditional machining process. Roya et al. [5] proposed a rough group analytic hierarchy process approach to the evaluation supplier criteria in the company for producing metal washers for the automotive industry. Vasiljević et al. [6] developed rough strength relational decision making and trial evaluation laboratory model to analyze the individual priorities of key success factors of hospital’s performance measures. As a generalization of the intuitionistic fuzzy set (IFS) [7], the picture fuzzy set (PFS) introduced by Cuong [8] is a very effective tool to express the complex fuzzy information because it is characterized by three functions expressing the degree of positive, neutral, and negative memberships at the same time. Because of this advantage, the PFS has been widely investigated and quite a few achievements have been made [9,10,11,12,13]. Among them, an important research topic in the research fields of MADM is aggregation operator theory, that can aggregate a collection of individual evaluated values into one. Abbas et al. [14] presented a comprehensive review on aggregation operator theory and decision-making approaches between 1986 and 2017. Among these aggregation operators, traditional aggregation operators, such as arithmetic and geometric operators for the IFS and neutral averaging operators [15] are based on the assumption that the attributes are independent of one another. However, the attributes of the problem are often correlative in the real decision-making process, especially in medical diagnosis. For example, to evaluate patients based on the following symptoms of lung diseases: (vital signs, body temperature, cough and hemoptysis), we want to place more emphasis on hemoptysis than on body temperature. However, on the other hand, we also want to pay more attention to patients who have severe hemoptysis and high body temperature, because hemoptysis and hyperthermia are two classical symptoms of pneumonia. Therefore, we need to find some new ways to deal with these situations where the decision data are correlative. The Choquet integral [16] introduced by Choquet is a useful tool to address the problem. Many scholars have made quite a few achievements in this field and applied the Choquet integral in MADM problems. By using Choquet integral and quasi-arithmetic means, Zhou and Chen [17] proposed a combined continuous quasi-arithmetic Choquet integral operator and a combined continuous generalized Choquet integral operator. In order to globally reflect the interactions between elements, Meng and Zhang [18] further defined the probabilistic generalized semivalue-induced continuous Choquet weighted averaging operator and the induced continuous Choquet geometric mean operator. Xu [19] used the Choquet integral to propose some operators for aggregating intuitionistic fuzzy values with correlative weights and further extended those operators to interval-valued intuitionistic fuzzy sets. Yager [20] proposed an approximation to the Choquet integral criteria aggregation that did not require ordering. By extending Marichal’s concept of entropy for fuzzy measures, Liu et al. [21] proposed a new method for determining fuzzy measures of the Choquet integral. Wen et al. [22] introduced Choquet integral-based linguistic operators under fuzzy heterogeneous environments for supplier selection in supply chain management. Some scholars also extended the Choquet integral to other fuzzy environments, such as in interval intuitionistic fuzzy information [23], the dual hesitant fuzzy environment [24], the interval-valued intuitionistic hesitant fuzzy environment [19] and the Pythagorean fuzzy environment [25]. Point operators are another aggregation tool to reduce the uncertainty of the aggregated arguments and thus obtain intensive information in the process of decision making. Since the point operator was proposed [26], it has been applied to many fields and has attracted increasing attention. Liu and Wang [27] proposed some point operators to translate IFS into another one. Xia and Xu [28] used the point operators to propose some operators for aggregating intuitionistic fuzzy values, and further extended those operators to intuitionistic multiplicative sets [29]. Peng [30,31], and Xing [32] also extended point operators to Pythagorean fuzzy sets, interval-valued Pythagorean fuzzy sets, and dual hesitant fuzzy sets, respectively.
However, the medical diagnosis problem in the real world is complex than many other applications. For instance: (1) We need to exactly express fuzzy information, and picture fuzzy numbers (PFNs) can depict doctors’ diagnoses for patients with respect to the symptoms; (2) We need to consider correlations among symptoms, and then the Choquet integral operator can be utilized to solve this problem; and (3) We need to reduce the uncertainty of doctor’s diagnosis data and get intensive information when diagnosing diseases. We can select point operators to achieve this function by adjusting the degree of doctor’s diagnosis data with some parameters. In order to solve above problems simultaneously, it is necessary to combine point operator with Choquet integral operator under picture fuzzy environment. Thus, the goal of this paper is to establish a new decision-making method that can not only control the certainty of doctor’s diagnosis data, but also deal with these situations where the diagnosis data are correlative. Then we apply new decision-making method to judge patient condition, and patients with different conditions are divided into different levels of hospitals instead of all patients rushing to large hospitals.
The rest of this paper is organized as follows. In the following section, we review some basic concepts related to PFS and the Choquet integral. In Section 3, we define some picture fuzzy point operators. In Section 4, by combining the point operators with Choquet integral operator, we propose the picture fuzzy point–Choquet averaging (PFPCA) operator, the picture fuzzy point–Choquet geometric (PFPCG) operator, the generalized picture fuzzy point–Choquet averaging (GPFPCA) operator and the generalized picture fuzzy point–Choquet geometric (GPFPCG) operator. Some prominent properties and special cases of these proposed operators are also studied. In Section 5, we introduce a novel method for solving MADM with picture fuzzy information based on the proposed operators. In Section 6, we provide an application example about assisting the hierarchical medical system to show the performance of new method.

2. Preliminaries

In the section, we briefly review some basic notions including PFS and the Choquet integral.

2.1. Picture Fuzzy Sets

Definition 1
[8]. Let X be an ordinary fixed set; then a picture fuzzy set P defined on X is given by
P = { x , μ p ( x ) , η p ( x ) , v p ( x ) | x X } ,
where μ p ( x ) is the positive degree of x X , and η p ( x ) and v p ( x ) are the neutral degree and negative degree, respectively, satisfying
μ p ( x ) + η p ( x ) + v p ( x ) 1 .
The uncertainty associated with PFS π P ( x ) = 1 μ p ( x ) η p ( x ) v p ( x ) is also defined. In the case η p ( x ) = 0 , PFS is reduced to the IFS, and when both μ p ( x ) , v p ( x ) = 0 , PFS is reduced to the fuzzy set.
For simplicity, we use the pair ( μ ( x ) , η ( x ) , v ( x ) ) to denote a general PFN that can be denoted by p = ( μ , η , v ) .
Given three PFNs p = ( μ , η , v ) , p 1 = ( μ 1 , η 1 , v 1 ) , p 1 = ( μ 2 , η 2 , v 2 ) , Cuong [8] defined the operations of intersection, union, complement and inclusion for them, which can be described as below:
p 1 p 2 = ( m i n ( μ 1 , μ 2 ) , m a x ( η 1 , η 2 ) , m a x ( v 1 , v 2 ) ) ,
p 1 p 2 = ( m a x ( μ 1 , μ 2 ) , m i n ( η 1 , η 2 ) , m i n ( v 1 , v 2 ) ) ,
p c = ( v , η , μ ) ,
p 1 p 2 ,   if   μ 1 μ 2 , η 1 η 2   and   v 1 v 2 .
Wei [9] further defines some operational laws for PFNs as shown below:
p 1 p 2 = ( ( μ 1 + μ 2 μ 1 μ 2 ) , η 1 η 2 , v 1 v 2 ) ,
p 1 p 2 = ( μ 1 μ 2 , η 1 + η 2 η 1 η 2 , v 1 + v 2 v 1 v 2 ) ,
λ p = ( 1 ( 1 μ ) λ , η λ , v λ ) ,
p λ = ( μ λ , 1 ( 1 η ) λ , 1 ( 1 v ) λ ) .
Definition 2
[13]. For two PFNs p 1 = ( μ 1 , η 1 , v 1 ) , p 2 = ( μ 2 , η 2 , v 2 ) , their relations are defined as follows:
p 1 p 2   iff   x X , μ 1 μ 2 , v 1 v 2 ,
p 1 = p 2   iff   x X , μ 1 = μ 2 , v 1 = v 2 .
In order to rank the PFNs, Garg [13] gave the score function and accuracy function of PFNs.
Definition 3
[13]. Suppose that p = ( μ , η , v ) is a PFN; then the score function of p is shown as follows:
S p = μ p v p .
Definition 4
[13]. Suppose that p = ( μ , η , v ) is a PFN; then the accuracy function of p is shown as follows:
H p = μ p + η p + v p .
Based on the score and accuracy function of PFN, Garg further defines the following ranking rules to compare two PFNs.
Definition 5.
For two PFNs:
if   S p 1 > S p 2 ,   then   p 1 > p 2 ,
if   S p 1 = S p 2 ,   then
if   H p 1 > H p 2 ,   t h e n   p 1 > p 2 ,
if   H p 1 = H p 2 ,   t h e n   p 1 = p 2 .

2.2. Choquet Integral Operator

The fuzzy measure can be used to define a weight on each combination of criteria in the Choquet integral model. In this subsection, we introduce the definitions of fuzzy measure and Choquet integral.
Definition 6
[33]. A fuzzy measure on X is a set function ρ : Γ ( x ) [ 0 , 1 ] , with the following conditions:
(1) 
ρ ( ϕ ) = 0 , ρ ( X ) = 1 (boundary conditions),
(2) 
A , B X and A B , then ρ ( A ) ρ ( B ) (monotonicity).
However, we generally need to determine 2 n 2 values for n criteria, which is quite complex, and thus it is not easy to give such fuzzy measure according to Definition 6. Therefore, the following σ -fuzzy measure ρ is further defined:
ρ ( A B )   =   ρ ( A ) + ρ ( B ) + σ ρ ( A ) ρ ( B ) ,
where A B = ϕ , and the parameter σ [ 1 , + ) denotes the interaction between attributes. In Equation (17):
(1)
If σ = 0, then σ -fuzzy measure ρ reduces to ρ ( A B ) = ρ ( A ) + ρ ( B ) , A B = ϕ , which is defined as an additive measure.
In this situation, if all the elements in X are independent, we get
ρ ( A ) = x i A ρ ( x i ) .
(2)
If all the elements in X are finite, then
ρ ( A ) = ρ ( i = 1 n x i ) = { 1 σ [ i = 1 n ( 1 + σ ρ ( x i ) ) 1 ] , σ 0 x i A ρ ( x i ) , σ = 0 ,
where x i x j = Φ , for i , j = 1 , 2 n , and i j .
(3)
If ρ 0 , then σ -fuzzy measure ρ reduces to ρ ( A B ) ρ ( A ) + ρ ( B ) , which is defined as a super-additive measure.
(4)
If 1 ρ 0 , then σ -fuzzy measure ρ reduces to ρ ( A B ) ρ ( A ) + ρ ( B ) , which is defined as a sub-additive measure.
When using a fuzzy measure to model the importance of decision criteria set S, a well-known aggregation function is the Choquet integral [16].
Definition 7.
Let f be a positive real-valued function on X and ρ be a fuzzy measure on X . The discrete Choquet integral of f with respect to ρ is defined as
( C ) f d ρ = i = 1 n [ ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ] f σ ( i ) ,
where σ ( i ) denotes a permutation of ( 1 , 2 n ) such that f σ ( 1 ) f σ ( 2 ) f σ ( n ) , and A σ ( 0 ) = ϕ , A σ ( i ) = { x σ ( 1 ) , x σ ( i ) } .

3. Some Point Operations for Picture Fuzzy Numbers and Their Properties

Motivated by the idea of intuitionistic fuzzy point operators [28] and dual hesitant fuzzy point operators [32], we will define a series of picture fuzzy point operations to obtain more intensive information and further analyze some desirable properties of these operations, which are very useful in the remainder of this paper.
Definition 8.
For a PFN p = ( μ , η , v ) , let α , β , γ [ 0 , 1 ] , we define some PF point operators: PFNPFN as follows:
D α , β ( p ) = { μ p + α π p , η p + β π p , v p + ( 1 α β ) π p } ,
F α , β , γ ( p ) = { μ p + α π p , η p + β π p , v p + γ π p } ,
where α + β + γ 1
G α , β , γ ( p ) = { α μ p , β η p , γ v p } .
It is obvious that the above PF point operators transform a PFN into another one. From Equations (21) and (22), we know that D α ( p ) assigns all the uncertainty into the other three parts of a PFS, while F α , β , γ ( p ) only assigns part of the uncertainty. Meanwhile, we can get π D α ( p ) = 1 π p , and π F α , β , γ ( p ) = π p ( 1 α β γ ) , which means that F α , β , γ ( p ) and D α ( p ) can reduce the uncertainty of PFS, and increase the positive degree, neutral degree, and positive degree. Similarly, From Equation (23), we know that G α , β , γ ( p ) can reduce the positive degree, neutral degree, and positive degree, and π G α , β , γ = ( 1 α μ p β η p γ v p ) , which means that G α , β , γ ( p ) increases the uncertainty of PFS.
Then, we discuss some properties of the operator F α , β , γ ( p ) in detail.
Theorem 1.
Let p = ( μ , η , v ) be a PFN and taking α , β , γ [ 0 , 1 ] , then
( F α , β , γ ( p c ) ) c = F γ , β , α ( p ) ,
( G α , β , γ ( p c ) ) c = G γ , β , α ( p ) .
If
α = μ p μ p + η p + v p ,   β = η p μ p + η p + v p   and   γ = v p μ p + η p + v p   then   F α , β , γ ( p ) = ( α , β , γ ) .
Proof. 
We prove the Equation (24) holds, and (25), (26) can be proved analogously.
(1)
From p c = ( v , η , μ ) , we get
( F α , β , γ ( p c ) ) c = ( v p + α π p , η p + β π p , μ p + γ π p ) c = F γ , β , α ( p ) .
(2)
Then
F α , β , γ ( p ) = ( μ p + μ p μ p + η p + v p π p , η p + η p μ p + η p + v p π p , v p + v p μ p + η p + v p π p ) = ( μ p μ p + η p + v p , η p μ p + η p + v p , v p μ p + η p + v p ) = ( α , β , γ ) .
Based on the operations of the PFNs, let D α 0 ( p ) = F α , β , γ 0 ( p ) = G α , β , γ 0 ( p ) = H α , β , γ 0 ( p ) = p ; we then get the following Theorem 2. ☐
Theorem 2.
Let p = ( μ , η , v ) be a PFN and taking α , β , γ [ 0 , 1 ] , and α + β + γ 0 , then
D α n ( γ ) = { μ p + α π p , η p + β π p , v p + ( 1 α β ) π p } ,
F α , β , γ n ( p ) = ( μ p + α π p τ , η p + β π p τ , v p + γ π p τ ) ,
where τ = 1 ( 1 α β γ ) n α + β + γ ,
G α , β , γ n ( p ) = ( μ p α n , η p β n , v p γ n ) .
The proof of this theorem is provided in Appendix A.
In the following, a numeric example is forwarded to illustrate Theorems 1 and 2.
Example 1.
Let p = ( 0.15 ,   0.35 ,   0.25 ) be a PFN, then the point operators of p can be calculated according to Definition 8 (Suppose α = 0.4 , β = 0.3 , γ = 0.2 ). Firstly, we can obtain π p = 1 ( 0.15 + 0.35 + 0.25 ) = 0.25 , and τ = 1 ( 1 α β γ ) n α + β + γ = 1 0.1 n 0.9 , then we have
D α , β ( p ) = { 0.15 + 0.25 α , 0.35 + 0.25 β , 0.25 + 0.25 ( 1 α β ) } = ( 0.25 ,   0.425 ,   0.325 ) ,
F α , β , γ ( p ) = ( 0.15 + 0.25 α , 0.35 + 0.25 β , 0.25 + 0.25 γ ) = ( 0.25 ,   0.425 ,   0.3 ) ,
G α , β , γ ( p ) = ( 0.15 α , 0.35 β , 0.25 γ ) = ( 0.06 ,   0.105 ,   0.05 ) .
Similarly,
D α n ( γ ) = ( 0.25 ,   0.425 ,   0.325 ) ,
F α , β , γ n ( p ) = ( 0.15 + 1 0.1 n 9 , 0.35 + 0.75 × ( 1 0.1 n ) 9 , 0.25 + 5 × ( 1 0.1 n ) 9 ) ,
G α , β , γ n ( p ) = { 0.15 × 0.4 n , 0.35 × 0.3 n , 0.25 × 0.2 n } .
From Theorem 2, we can easily obtain the following properties.
Theorem 3.
Let p = ( μ , η , v ) be a PFS, and n be a positive integer. Taking α , β , γ [ 0 , 1 ] , then
( F α , β , γ n ( p c ) ) c = F γ , β , α n ( p ) ,
( G α , β , γ n ( p c ) ) c = G γ , β , α n ( p ) .
Theorem 4.
Let p = ( μ , η , v ) be a PFS, and n be a positive integer. Taking α , β , γ [ 0 , 1 ] , the relation is defined as A B if and only if μ F α , β , γ n ( p ) μ F α , β , γ n 1 ( p ) , and v F α , β , γ n ( p ) v F α , β , γ n 1 ( p ) , and then
( F α , β , γ n ( p c ) ) c = F γ , β , α n ( p ) ,
π F α , β , γ n π F α , β , γ n 1 .
If   α = μ p μ p + η p + v p , β = η p μ p + η p + v p , γ = v p μ p + η p + v p ,   then   F α , β , γ n ( p ) = F α , β , γ ( γ ) .
Definition 9.
Let α , β , γ [ 0 , 1 ] , and α + β + γ ≤ 1. We define the following limit:
lim n F α , β , γ n ( p ) = lim n { μ F α , β , γ n ( p ) , η F α , β , γ n ( p ) , v F α , β , γ n ( p ) } .
Theorem 5.
Let α , β , γ [ 0 , 1 ] , and α + β + γ 1 ; then we have
lim n F α , β , γ n ( p ) = = D α α + β + γ , β α + β + γ n ( p ) .
Proof of Theorem 5.
According to Theorem 7, we get
lim n μ F ξ ζ ( γ ) n = lim n ( μ p + α π p 1 ( 1 α β γ ) n α + β + γ ) = μ p + α α + β + γ π p ,
lim n η F ξ ζ ( γ ) n = lim n ( η p + β π p 1 ( 1 α β γ ) n α + β + γ ) = η p + β α + β + γ π p ,
lim n v F ξ ζ ( γ ) n = = v p + γ α + β + γ π p .
So we have
lim n F α , β , γ n ( p ) = lim n { μ F α , β , γ n ( p ) , η F α , β , γ n ( p ) , v F α , β , γ n ( p ) } = { μ p + α α + β + γ π p , η p + β α + β + γ π p , v p + γ α + β + γ π p } = D α α + β + γ , β α + β + γ n ( p )
. ☐

4. Picture Fuzzy Point–Choquet Integral Aggregation Operators and Their Properties

In order to get more intensive information from PFS and efficiently deal with correlations among arguments at the same time, we combine picture fuzzy point operators with the Choquet integral operator to propose some new class of aggregation operators for aggregating picture fuzzy information in this section. Some desirable properties of proposed aggregation operators are also discussed in detail.

4.1. Picture Fuzzy Point–Choquet Averaging Operator

Definition 10.
Let Ω be the set of all PFNs, and p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs, taking α i , β i , γ i [ 0 , 1 ] . Then we define the series of PFPCA operators): Ω m Ω , if
F ( C 1 ) p d ρ = P F P C A D α , β n ( p 1 , p 2 p n ) = i = 1 m ( ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ) D α σ ( i ) , β σ ( i ) n ( p σ ( i ) ) ,
F ( C 2 ) p d ρ = P F P C A F α , β , γ n ( p 1 , p 2 p n ) = i = 1 m ( ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ) F α σ ( i ) , β σ ( i ) , γ σ ( i ) n ( p σ ( i ) ) ,
F ( C 3 ) p d ρ = P F P C A G α , β , γ n ( p 1 , p 2 p n ) = i = 1 m ( ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ) G α σ ( i ) , β σ ( i ) , γ σ ( i ) n ( p σ ( i ) ) ,
where σ ( i ) denotes a permutation of ( 1 , 2 m ) such that p σ ( 1 ) p σ ( 2 ) p σ ( m ) , and G σ ( i ) is the attribute corresponding to p σ ( i ) .
By operational laws defined in Section 2.1, we can obtain the following theorem.
Theorem 6.
Let p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs, and σ ( i ) be a permutation of ( 1 , 2 m ) such that p σ ( 1 ) p σ ( 2 ) p σ ( m ) , G σ ( i ) is the attribute corresponding to p σ ( i ) , and A σ ( 0 ) = ϕ A σ ( i ) = { G σ ( 1 ) , G σ ( i ) } , taking ω ˜ i = ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) , α i , β i , γ i [ 0 , 1 ] , α i + β i + γ i 1 . Then, the aggregated values by the series of PFPCA operators are also PFNs:
P F P C A D α , β n ( p 1 , p 2 p n ) = { 1 i = 1 m ( 1 ( μ p σ ( i ) + α i π p σ ( i ) ) ) ω ˜ i , 1 i = 1 m ( 1 ( η p σ ( i ) + β i π p σ ( i ) ) ) ω ˜ i i = 1 m ( v p σ ( i ) + ( 1 α i β i ) π p σ ( i ) ) ω ˜ i } ,
P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) = ( ( 1 i = 1 m ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i ) , i = 1 m η F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i , i = 1 m v F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i )
where
μ F α i , β i , γ i n ( p i ) = μ p σ ( i ) + α i π p σ ( i ) 1 ( 1 α i β i γ i ) n α i + β i + γ i ,
η F α i , β i , γ i n ( p i ) = η p σ ( i ) + β i π p σ ( i ) 1 ( 1 α i β i γ i ) n α i + β i + γ i ,
v F α i , β i , γ i n ( p σ ( i ) ) = v p σ ( i ) + γ i π p σ ( i ) 1 ( 1 α i β i γ i ) n α i + β i + γ i ,
P F P C A G α , β , γ n ( p 1 , p 2 , , p m ) = ( ( 1 i = 1 m ( 1 μ p σ ( i ) α i n ) ω ˜ i ) , i = 1 m ( η p σ ( i ) β i n ) ω ˜ i , i = 1 m ( v p σ ( i ) γ i n ) ω ˜ i ) .
The proof of this theorem is provided in Appendix B.
In the following, a numeric example is forwarded to illustrate Theorem 6.
Example 2.
Let p 1 = ( 0.25 ,   0.35 ,   0.15 ) , p 2 = ( 0.42 ,   0.18 ,   0.37 ) , p 3 = ( 0.34 ,   0.27 ,   0.16 ) be PFN. Then we aggregate the three PFNs by the following steps:
Step 1. Identify the fuzzy measure of the n attributes of G according to Equations (17) and (19). Suppose that the fuzzy measures of attributes of G are given as follows:
ρ ( G 1 ) = 0.38 , ρ ( G 2 ) = 0.27 , ρ ( G 3 ) = 0.36 .
Firstly, according to Equation (19), the value of σ is obtained: σ = 0.029 , and then the fuzzy measures of attribute sets of G = { G 1 , G 2 , G 3 , G 4 } can be calculated by Equation (13), shown as follows:
ρ ( G 1 , G 2 ) = 0.65 , ρ ( G 1 , G 3 ) = 0.74 ,   ρ ( G 2 , G 3 ) = 0.63 ,   ρ ( G 1 , G 2 , G 3 ) = 1 .
Step 2. By score functions, we rearrange the three PFNs in descending order, shown as follows:
s ( p 1 ) = 0.1 , s ( p 2 ) = 0.05 , s ( p 3 ) = 0.18 ,
p σ ( 1 ) = ( 0.34 ,   0.27 ,   0.16 ) , p σ ( 2 ) = ( 0.25 ,   0.35 ,   0.15 ) , p σ ( 3 ) = ( 0.42 ,   0.18 ,   0.37 ) .
Then we can get
A σ ( 1 ) = { G 3 } , A 1 σ ( 2 ) = { G 1 , G 3 } , A 1 σ ( 4 ) = { G 1 , G 2 , G 3 } ,
ρ A σ ( 1 ) = ρ G 3 = 0.36 , ρ A σ ( 2 ) ρ A σ ( 1 ) = ρ G 1 G 3 ρ G 3 = 0.38 , ρ A σ ( 3 ) ρ A σ ( 2 ) = ρ G 2 G 3 G 4 ρ G 1 G 3 = 0.26 .
.
Step 3. Calculate the point operators of pi according to Definition 8 (Suppose α = 0.3 , β = 0.4 , γ = 0.1 , n = 3 ). Firstly, we can obtain π p σ ( 1 ) = 1 ( 0.25 + 0.35 + 0.15 ) = 0.25 , π p σ ( 2 ) = 1 ( 0.42 + 0.18 + 0.37 ) = 0.03 , π p σ ( 3 ) = 1 ( 0.34 + 0.27 + 0.16 ) = 0.23 . Then we have
F α , β , γ 3 ( p σ ( 1 ) ) = ( 0.25 + 0.25 × 1 0.2 3 0.8 , 0.35 + 0.25 × 1 0.2 3 0.8 , 0.15 + 0.25 × 1 0.2 3 0.8 ) = ( 0.56 , 0.66 , 0.46 ) ,
F α , β , γ 3 ( p σ ( 2 ) ) = ( 0.42 + 0.03 × 1 0.2 3 0.8 , 0.18 + 0.03 × 1 0.2 3 0.8 , 0.37 + 0.03 × 1 0.2 3 0.8 ) = ( 0.46 , 0.22 , 0.41 ) , F α , β , γ 3 ( p σ ( 3 ) ) = ( 0.34 + 0.23 × 1 0.2 3 0.8 , 0.27 + 0.23 × 1 0.2 3 0.8 , 0.16 + 0.23 × 1 0.2 3 0.8 ) = ( 0.63 , 0.56 , 0.45 )
Step 4. Utilize the P F P C A F α , β , γ n ( p 1 , p 2 , p 3 ) operator to aggregate the three PFNs and get the aggregated p as follows:
p = i = 1 m ( ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ) F α i , β i , γ i n ( p σ ( i ) ) = ( ( 1 i = 1 m ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ) , i = 1 m ( η F α i , β i , γ i n ( p σ ( i ) ) ) ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) , i = 1 m ( v F α i , β i , γ i n ( p σ ( i ) ) ) ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ) = ( ( 1 ( 1 0.56 ) 0.36 × ( 1 0.46 ) 0.38 × ( 1 0.63 ) 0.26 ) , 0.66 0.36 × 0.22 0.38 × 0.56 0.26 , 0.46 0.36 × 0.41 0.38 × 0.45 0.26 ) .
Example 2 gives a detailed portrait of the P F P C A F α , β , γ n ( p 1 , p 2 , p 3 ) operator. It should be pointed out that the P F P C A F α , β , γ n ( p 1 , p 2 , p 3 ) operator includes a reorder step and it is similar to the famous ordered weighted averaging (OWA) operator. In the following, we discuss some properties of the above PFPCA operators.
Theorem 7.
Let p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs. Taking k 0 , then
P F P C A D α , β n ( k p 1 , k p 2 , , k p m ) = k P F P C A D α , β n ( p 1 , p 2 , , p m ) ,
P F P C A F α , β , γ n ( k p 1 , k p 2 , , k p m ) = k P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) ,
P F P C A G α , β , γ n ( k p 1 , k p 2 , , k p m ) = k P F P C A G α , β , γ n ( p 1 , p 2 , , p m ) .
Proof. 
We prove the Equation (50) holds for all m, and the others can be proved analogously.
By the operational law in Section 2.2, we have
k p i = ( 1 ( 1 μ i ) k , η i k , v i k )  
and
P F P C A F α , β , γ n ( k p 1 , k p 2 , , k p m ) = ( 1 i = 1 m ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) k ω ˜ i , i = 1 m η F α i , β i , γ i n ( p i ) k ω ˜ i , i = 1 m v F α i , β i , γ i n ( p σ ( i ) ) k ω ˜ i ) ,
and hence
k P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) = k ( 1 i = 1 m ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i , i = 1 m η F α i , β i , γ i n ( p σ ( i ) ) k ω ˜ i , i = 1 m v F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i ) ( 1 i = 1 m ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) k ω ˜ i , i = 1 m η F α i , β i , γ i n ( p σ ( i ) ) k ω ˜ i , i = 1 m v F α i , β i , γ i n ( p σ ( i ) ) k ω ˜ i ) = P F P C A F α , β , γ n ( k p 1 , k p 2 , , k p m ) .
Therefore, Equation (50) holds, which completes the proof. ☐
Theorem 8.
Let pi and qi be two collections of PFNs, then
P F P C A D α , β n ( p 1 q 1 , p 2 q 2 , , p n q m ) = P F P C A D α , β n ( p 1 , p 2 , , p m ) P F P C A D α , β n ( q 1 , q 2 , , q m )
P F P C A F α , β , γ n ( p 1 q 1 , p 2 q 2 , , p n q m ) = P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) P F P C A F α , β , γ n ( q 1 , q 2 , , q m )
P F P C A G α , β , γ n ( p 1 q 1 , p 2 q 2 , , p n q m ) = P F P C A G α , β , γ n ( p 1 , p 2 , , p m ) P F P C A G α , β , γ n ( q 1 , q 2 , , q m )
Proof. 
We prove the Equation (53) holds for all m, and the others can be proved analogously.
By the operational law in Section 2.2, we have
p i q i = ( μ p i + μ q i μ p i μ q i , η p i η q i , v p i v q i ) ,
P F P C A F α , β , γ n ( p 1 q 1 , p 2 q 2 , , p n q m ) = ( ( 1 i = 1 m ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i ( 1 μ F α i , β i , γ i n ( q σ ( i ) ) ) ω ˜ i ) , i = 1 m η F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i η F α i , β i , γ i n ( q σ ( i ) ) ω ˜ i , i = 1 m v F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i v F α i , β i , γ i n ( q σ ( i ) ) ω ˜ i ) ,
P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) P F P C A F α , β , γ n ( q 1 , q 2 , , q m ) = ( 1 i = 1 m ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i ( 1 μ F α i , β i , γ i n ( q σ ( i ) ) ) ω ˜ i , i = 1 m η F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i η F α i , β i , γ i n ( q σ ( i ) ) ω ˜ i , i = 1 m ν F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i ν F α i , β i , γ i n ( q σ ( i ) ) ω ˜ i ) = P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) P F P C A F α , β , γ n ( q 1 , q 2 , , q m ) .
Therefore, Equation (53) holds, which completes the proof. ☐
Theorem 9.
(Idempotency). If p i = ( μ i , η i , v i ) are equal, i.e., p i = p = ( μ , η , v ) for all i, then
P F P C A D α , β n ( p 1 , p 2 , , p m ) = D α , β n ,
P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) = F α , β , γ n ,
P F P C A G α , β , γ n ( p 1 , p 2 , , p m ) = G α , β , γ n .
Proof. 
We prove the Equation (56) holds for all m, and the others can be proved analogously.
Since p i = p = ( μ , η , v ) for all i, then
P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) = ( ( 1 i = 1 m ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i ) , i = 1 m η F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i , i = 1 m v F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i ) = ( 1 ( 1 μ F α i , β i , γ i n ( p ) ) i = 1 m ω i , ( η F α i , β i , γ i n ( p ) ) i = 1 m ω i , ( v F α i , β i , γ i n ( p ) ) i = 1 m ω i ) = ( 1 ( 1 μ F α i , β i , γ i n ( p ) ) , η F α i , β i , γ i n ( p ) , v F α i , β i , γ i n ( p ) ) = ( μ F α i , β i , γ i n ( p ) , η F α i , β i , γ i n ( p ) , v F α i , β i , γ i n ( p ) ) = F α , β , γ n .
 ☐
Theorem 10.
(Monotonicity) Let p i = ( μ p i , η p i , v p i ) and q i = ( μ q i , η q i , v q i ) ( i = 1 , 2 , , m ) be two collections of PFN. If μ p i μ q i , η p i η q i and v p i v q i holds for all i ( i = 1 , 2 , , m ) , then
P F P C A D α , β n ( p 1 , p 2 , , p m ) P F P C A D α , β n ( q 1 , q 2 , , q m ) ,
P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) P F P C A F α , β , γ n ( q 1 , q 2 , , q m ) ,
P F P C A G α , β , γ n ( p 1 , p 2 , , p m ) P F P C A G α , β , γ n ( q 1 , q 2 , , q m ) .
Proof. 
We prove the Equation (59) holds for all m, and the others can be proved analogously.
By Theorem 6, we get
P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) = ( ( 1 i = 1 m ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i ) , i = 1 m η F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i , i = 1 m v F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i ) , P F P C A F α , β , γ n ( q 1 , q 2 , , q m ) = ( ( 1 i = 1 m ( 1 μ F α i , β i , γ i n ( q σ ( i ) ) ) ω ˜ i ) , i = 1 m η F α i , β i , γ i n ( q σ ( i ) ) ω ˜ i , i = 1 m v F α i , β i , γ i n ( q σ ( i ) ) ω ˜ i ) .
Since μ p i μ q i and v p i v q i , we can get
( 1 i = 1 m ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i ) ( 1 i = 1 m ( 1 μ F α i , β i , γ i n ( q σ ( i ) ) ) ω ˜ i )  
and
i = 1 m η F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i i = 1 m η F α i , β i , γ i n ( q σ ( i ) ) ω ˜ i , i = 1 m v F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i i = 1 m v F α i , β i , γ i n ( q σ ( i ) ) ω ˜ i .
By Definition 6, we get P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) P F P C A F α , β , γ n ( q 1 , q 2 , , q m ) . ☐
Theorem 11.
(Boundedness) Let p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs, then
d D α , β n P F P C A D α , β n d D α , β n + ,
d F α , β , γ n P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) d F α , β , γ n + ,
d G α , β , γ n P F P C A G α , β , γ n ( p 1 , p 2 , , p m ) d G α , β , γ n + ,
where d Δ + = ( m a x i ( μ Δ ) , m i n i ( v Δ ) ) and d Δ = ( m i n i ( μ Δ ) , m a x i ( v Δ ) ) and Δ denotes D α , β n , F α , β , γ n , G α , β , γ n .
Proof. 
We prove the Equation (62) holds for all m, and the others can be proved analogously.
From Theorem 6, we can get
P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) = ( ( 1 i = 1 m ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i ) , i = 1 m η F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i , i = 1 m v F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i ) .
By the definition of d F ξ , ζ n + , d F ξ , ζ n we can get
1 i = 1 m ( 1 m i n ( μ F α i , β i , γ i n ( p σ ( i ) ) ) ) ω ˜ i 1 i = 1 m ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i 1 i = 1 m ( 1 m a x ( μ F α i , β i , γ i n ( p σ ( i ) ) ) ) ω ˜ i F P C A F α , β , γ n ( p 1 , p 2 , , p m ) P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) F P C A F α , β , γ n ( p 1 + , p 2 + , , p m + ) .
By Definition 7, we get d F α , β , γ n P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) d F α , β , γ n + . ☐
By giving different values of the parameters, we get the following special cases.
Theorem 12.
Let p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs, then
(1) 
If ω i = ω ˜ i = ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) , then the series of PFPCA operators are all reduced to the series of picture fuzzy point averaging (PFPWA) operators. In particular, if m i = 1 m , ( i = 1 , 2 , , m ) , then PFPCA operators is reduced to a picture fuzzy averaging (PFA) operator, which is defined as:
P F A = ( 1 ( i = 1 m ( 1 μ i 2 ) ) 1 / m , ( i = 1 m η i ) 1 / m , ( i = 1 m v i ) 1 / m ) .  
(2) 
If n = 0 , ω ˜ i = ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) , and ρ ( A ) = i = 1 | A | ω ˜ i for all A X , where | A | is the number of the elements in set A, ω ˜ = ( ω ˜ 1 , ω ˜ 2 , ω ˜ m ) T , ω ˜ i [ 0 , 1 ] , i = 1 m ω ˜ i = 1 , then the PFPCA operator is reduced to a picture fuzzy order-weighted averaging (PFOWA) operator defined by Garg [13].
(3) 
If n = 0 , ω ˜ i = ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) , then the series of PFPCA operators are all reduced to the series of picture fuzzy weighted averaging (PFWA) operators defined by Garg [13].

4.2. Picture Fuzzy Point–Choquet Geometric Operator

Definition 11.
Let Ω be the set of all PFNs, and p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs, taking α i , β i , γ i [ 0 , 1 ] . Then we define the series of PFPCG operators: Ω m Ω , if
F ( C 4 ) p d ρ = P F P C G D α , β n ( p 1 , p 2 p m ) = i = 1 m ( ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ) D α σ ( i ) , β σ ( i ) n ( p σ ( i ) ) ,
F ( C 5 ) p d ρ = P F P C G F α , β , γ n ( p 1 , p 2 p m ) = i = 1 m ( ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ) F α σ ( i ) , β σ ( i ) , γ σ ( i ) n ( p σ ( i ) ) ,
F ( C 6 ) p d ρ = P F P C G G α , β , γ n ( p 1 , p 2 p m ) = i = 1 m ( ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ) G α σ ( i ) , β σ ( i ) , γ σ ( i ) n ( p σ ( i ) ) ,
where σ ( i ) denotes a permutation of ( 1 , 2 m ) such that p σ ( 1 ) p σ ( 2 ) p σ ( m ) , and G σ ( i ) is the attribute corresponding to p σ ( i ) , A σ ( i ) = { G σ ( 1 ) , G σ ( i ) } , A σ ( 0 ) = ϕ .
By operational laws defined in Section 2.1, we can obtain the following theorem.
Theorem 13.
Let p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs, and σ ( i ) be a permutation of ( 1 , 2 m ) such that p σ ( 1 ) p σ ( 2 ) p σ ( m ) , G σ ( i ) is the attribute corresponding to p σ ( i ) , A σ ( i ) = { G σ ( 1 ) , G σ ( i ) } , A σ ( 0 ) = ϕ . Taking ω ˜ i = ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) , α i , β i , γ i [ 0 , 1 ] , α i + β i + γ i 1 , then the aggregated values by the series of PFPCG operators are also PFNs, and
P F P C G D α , β n ( p 1 , p 2 p m ) = { i = 1 m ( μ p σ ( i ) + ( 1 γ i β i ) π p σ ( i ) ) ω ˜ i , 1 i = 1 m ( 1 ( η p σ ( i ) + β i π p σ ( i ) ) ) ω ˜ i , i = 1 m ( 1 ( v p σ ( i ) + γ i π p σ ( i ) ) ) ω ˜ i }
P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) = ( i = 1 m μ F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i , 1 i = 1 m ( 1 η F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i , 1 i = 1 m ( 1 v F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i )  
where
μ F α i , β i , γ i n ( p i ) = μ p σ ( i ) + α i π p σ ( i ) 1 ( 1 α i β i γ i ) n α i + β i + γ i ,
η F α i , β i , γ i n ( p i ) = η p σ ( i ) + β i π p σ ( i ) 1 ( 1 α i β i γ i ) n α i + β i + γ i ,
v F α i , β i , γ i n ( p σ ( i ) ) = v p σ ( i ) + γ i π p σ ( i ) 1 ( 1 α i β i γ i ) n α i + β i + γ i
P F P C G G α , β , γ n ( p 1 , p 2 , , p m ) = ( ( 1 i = 1 m ( 1 μ p σ ( i ) α i n ) ω ˜ i ) , i = 1 m ( η p σ ( i ) β i n ) ω ˜ i , i = 1 m ( v p σ ( i ) γ i n ) ω ˜ i )
Theorem 14.
Let p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs. Taking k 0 , then
P F P C G D α , β n ( p 1 k , p 2 k , , p m k ) = ( P F P C G D α , β n ( p 1 , p 2 , , p m ) ) k ,
P F P C G F α , β , γ n ( p 1 k , p 2 k , , p m k ) = ( P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) ) k ,
P F P C G G α , β , γ n ( p 1 k , p 2 k , , p m k ) = ( P F P C G G α , β , γ n ( p 1 , p 2 , , p m ) ) k .
Theorem 15.
Let p i = ( μ p i , η p i , v p i ) and q i = ( μ q i , η q i , v q i ) ( i = 1 , 2 , , m ) be two collections of PFNs, then
P F P C G D α , β n ( p 1 q 1 , p 2 q 2 , , p m q m ) = P F P C G D α , β n ( p 1 , p 2 , , p m ) P F P C G D α , β n ( q 1 , q 2 , , q m ) ,
P F P C G F α , β , γ n ( p 1 q 1 , p 2 q 2 , , p m q m ) = P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) P F P C G F α , β , γ n ( q 1 , q 2 , , q m ) ,
P F P C G G α , β , γ n ( p 1 q 1 , p 2 q 2 , , p m q m ) = P F P C G G α , β , γ n ( p 1 , p 2 , , p m ) P F P C G G α , β , γ n ( q 1 , q 2 , , q m ) .
Parallel to Theorems 9–11, the series of PFPCG operators have properties similar to PFPCA operators such as idempotency, monotonicity, and boundedness under some conditions, which are omitted in order to save space.
ω i = ω ˜ i = ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) .  

4.3. Generalized Picture Fuzzy Point–Choquet Averaging Operator

Definition 12.
Let p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs, taking α i , β i , γ i [ 0 , 1 ] , λ 0 , and α i + β i + γ i 1 . Then we define a series of GPFPCA operators: Ωm → Ω, if
F ( C 7 ) p d ρ = G P F P C A D α , β n ( p 1 , p 2 p n ) = ( i = 1 m ( ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ) ( D α σ ( i ) , β σ ( i ) n ( p σ ( i ) ) ) λ ) 1 / λ ,
F ( C 8 ) p d ρ = G P F P C A F α , β , γ n ( p 1 , p 2 p n ) = ( i = 1 m ( ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ) ( F α σ ( i ) , β σ ( i ) , γ σ ( i ) n ( p σ ( i ) ) ) λ ) 1 / λ ,
F ( C 9 ) p d ρ = G P F P C A G α , β , γ n ( p 1 , p 2 p n ) = ( i = 1 m ( ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ) ( G α σ ( i ) , β σ ( i ) , γ σ ( i ) n ( p σ ( i ) ) ) λ ) 1 / λ .
where G σ ( i ) is the attribute corresponding to p σ ( i ) , A σ ( i ) = { G σ ( 1 ) , G σ ( i ) } , A σ ( 0 ) = ϕ , and σ ( i ) denotes a permutation of ( 1 , 2 , , m ) such that p σ ( 1 ) p σ ( 2 ) p σ ( m ) .
By operational laws defined in Section 2.1, we can obtain the following theorem.
Theorem 16.
Let p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs, and σ ( i ) be a permutation of ( 1 , 2 m ) such that p σ ( 1 ) p σ ( 2 ) p σ ( m ) , G σ ( i ) is the attribute corresponding to p σ ( i ) , A σ ( 0 ) = ϕ , A σ ( i ) = { G σ ( 1 ) , G σ ( i ) } , and taking ω ˜ i = ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) , α i , β i , γ i [ 0 , 1 ] , α i + β i + γ i 1 , then the aggregated values by the series of GPFPCA operators are also PFNs.
(1)
G P F P C A D α , β n ( p 1 , p 2 p n ) = { 1 i = 1 m ( 1 ( μ p σ ( i ) + α i π p σ ( i ) ) λ ) ω ˜ i , 1 ( 1 i = 1 m ( 1 ( 1 η p σ ( i ) ( 1 β i ) π p σ ( i ) ) λ ) ω ˜ j ) 1 / λ ,   1 ( 1 i = 1 m ( 1 ( 1 v p σ ( i ) ( 1 γ i ) π p σ ( i ) ) λ ) ω ˜ j ) 1 / λ } ;
(2)
G P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) = ( ( 1 i = 1 m ( 1 μ F α i , β i , γ i n ( p σ i ) λ ) ω ˜ i ) 1 λ , 1 [ 1 i = 1 m ( 1 ( 1 η F α i , β i , γ i n ( p σ i ) ) λ ) ω ˜ i ] 1 λ , 1 [ 1 i = 1 m ( 1 ( 1 v F α i , β i , γ i n ( p σ i ) ) λ ) ω ˜ i ] 1 λ )
where
μ F α i , β i , γ i n ( p i ) = μ p σ ( i ) + α i π p σ ( i ) 1 ( 1 α i β i γ i ) n α i + β i + γ i ,
η F α i , β i , γ i n ( p i ) = η p σ ( i ) + β i π p σ ( i ) 1 ( 1 α i β i γ i ) n α i + β i + γ i ,
v F α i , β i , γ i n ( p σ ( i ) ) = v p σ ( i ) + γ i π p σ ( i ) 1 ( 1 α i β i γ i ) n α i + β i + γ i ;
(3) G P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) =
( ( 1 i = 1 m ( 1 μ p σ ( i ) λ α i n λ ) ω ˜ i ) 1 λ , 1 [ 1 i = 1 m ( 1 ( 1 η p σ ( i ) β i n ) λ ) ω ˜ i ] 1 λ , 1 [ 1 i = 1 m ( 1 ( 1 v p σ ( i ) γ i n ) λ ) ω ˜ i ] 1 λ ) .
Parallel to Theorems 9–11, the series of GPFPCA operators have properties similar to PFPCA operators such as idempotency, monotonicity, and boundedness under some conditions, which are omitted in order to save space.

4.4. Generalized Picture Fuzzy Point–Choquet Geometric Ooperator

Definition 13.
Let p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs, taking α i , β i , γ i [ 0 , 1 ] , λ 0 , and α i + β i + γ i 1 . Then we define a series of GPFPCG operators: Ωm → Ω, if
F ( C 10 ) p d ρ = G P F P C G D α , β n ( p 1 , p 2 p m ) = 1 λ i = 1 m ( λ D α σ ( i ) , β σ ( i ) n ( p σ ( i ) ) ) ( ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ) ,
F ( C 11 ) p d ρ = G P F P C G F α , β , γ n ( p 1 , p 2 p m ) = 1 λ i = 1 m ( λ F α σ ( i ) , β σ ( i ) γ σ ( i ) n ( p σ ( i ) ) ) ( ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ) ,
F ( C 12 ) p d ρ = G P F P C G G α , β , γ n ( p 1 , p 2 p m ) = 1 λ i = 1 m ( λ G α σ ( i ) , β σ ( i ) γ σ ( i ) n ( p σ ( i ) ) ) ( ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) ) .
Similarly, we can obtain the following theorem:
Theorem 17.
Let p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs, and σ ( i ) be a permutation of ( 1 , 2 m ) such that p σ ( 1 ) p σ ( 2 ) p σ ( m ) , G σ ( i ) is the attribute corresponding to p σ ( i ) , A σ ( 0 ) = ϕ , and A σ ( i ) = { G σ ( 1 ) , G σ ( i ) } . Taking ω ˜ i = ρ ( A σ ( i ) ) ρ ( A σ ( i 1 ) ) , α i , β i , γ i [ 0 , 1 ] , α i + β i + γ i 1 , then the aggregated values by the series of GPFPCG operators are also PFNs, and
(1)
G P F P C G D α , β n ( p 1 , p 2 p n ) = { 1 ( 1 i = 1 m ( 1 ( 1 μ p σ ( i ) ( 1 α i ) π p σ ( i ) ) λ ) ω ˜ i ) 1 / λ , 1 i = 1 m ( 1 ( η p σ ( i ) + β i π p σ ( i ) ) λ ) ω ˜ i , 1 i = 1 m ( 1 ( v p σ ( i ) + γ i π p σ ( i ) ) λ ) ω ˜ i } ;
(2)
G P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) = { 1 ( 1 i = 1 m ( 1 ( 1 μ p σ ( i ) ( 1 α σ ( i ) ) π p σ ( i ) ) λ ) ω ˜ i ) 1 / λ , 1 i = 1 m ( 1 ( η p σ ( i ) + β σ ( i ) π p σ ( i ) ) λ ) ω ˜ i , 1 i = 1 m ( 1 ( v p σ ( i ) + γ σ ( i ) π p σ ( i ) ) λ ) ω ˜ i } ,
where
μ F α i , β i , γ i n ( p i ) = μ p σ ( i ) + α i π p σ ( i ) 1 ( 1 α i β i γ i ) n α i + β i + γ i ,
η F α i , β i , γ i n ( p i ) = η p σ ( i ) + β i π p σ ( i ) 1 ( 1 α i β i γ i ) n α i + β i + γ i , v F α i , β i , γ i n ( p σ ( i ) ) = v p σ ( i ) + γ i π p σ ( i ) 1 ( 1 α i β i γ i ) n α i + β i + γ i ;
(3)
G P F P C G G α , β , γ n ( p 1 , p 2 , , p m ) = ( 1 [ 1 i = 1 m ( 1 ( 1 μ p σ ( i ) α i n ) λ ) ω i ] 1 λ , ( 1 i = 1 m ( 1 η p σ ( i ) λ β i n λ ) ω ˜ i ) 1 λ , ( 1 i = 1 m ( 1 v p σ ( i ) λ γ i n λ ) ω ˜ i ) 1 λ ) .
Parallel to Theorems 13–15, the series of GPFPCG operators have properties such as idempotency, monotonicity, and boundedness under some conditions, which are omitted in order to save space.
In fact, the correlations of these proposed aggregation operators can be further studied. Here, we take PFPCAF α , β , γ n as an example.
Theorem 18.
Let p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs, Then the operation of complement on them is as follows:
P F P C A F α , β , γ n ( p 1 c , p 2 c , , p m c ) = P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) c ,
P F P C G F α , β , γ n ( p 1 c , p 2 c , , p m c ) = P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) c ,
G P F P C A F α , β , γ n ( p 1 c , p 2 c , , p m c ) = G P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) c ,
G P F P C G F α , β , γ n ( p 1 c , p 2 c , , p m c ) = G P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) c .
By Theorems 3–5, we can easily obtain the following theorems.
Theorem 19.
Let p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs, then the operation of the complement to aggregation operators is as follows:
[ P F P C A F α , β , γ n ( p 1 c , p 2 c , , p m c ) ] c = P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) ,
[ P F P C G F α , β , γ n ( p 1 c , p 2 c , , p m c ) ] c = P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) ,
[ G P F P C G F α , β , γ n ( p 1 c , p 2 c , , p m c ) ] c = G P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) ,
[ G P F P C A F α , β , γ n ( p 1 c , p 2 c , , p m c ) ] c = G P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) .
Theorem 20.
Let p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs, then
lim n P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) = P F P C A D α α + β + γ , β α + β + γ , n ( p 1 , p 2 , , p m ) ,
lim n P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) = P F P C G D α α + β + γ , β α + β + γ , n ( p 1 , p 2 , , p m ) ,
lim n G P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) = G P F P C A D α α + β + γ , β α + β + γ , n ( p 1 , p 2 , , p m ) ,
lim n G P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) = G P F P C G D α α + β + γ , β α + β + γ , n ( p 1 , p 2 , , p m ) .
Theorem 21.
Let p i = ( μ i , η i , v i ) be a collection of PFNs, If α i = μ p i μ p i + η p i + v p i , β i = η p i μ p i + η p i + v p i , γ i = v p i μ p i + η p i + v p i , ( i = 1 , 2 , , m ) then
P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) = P F P C A F α , β , γ ( p 1 , p 2 , , p m ) ,
P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) = P F P C G F α , β , γ ( p 1 , p 2 , , p m ) ,
G P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) = G P F P C A F α , β , γ ( p 1 , p 2 , , p m ) ,
G P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) = G P F P C G F α , β , γ ( p 1 , p 2 , , p m ) .
In following, we discuss the differences and relationships between PFPCA, PFPCG, GPFPCA, and GPFPCG operators in detail.
In the case where λ = 1 , the GPFPCA operator reduces to the PFPCA operator in Definition 10, and the GPFPCG operator reduces to the PFPCG operator in Definition 11. On the other hand, the PFPCA operator is an arithmetic aggregation operator, and thus the PFPCG operator can be treated its geometric form. Similarly, GPFPCG operator is geometric form of GPFPCA operator. Since m x i λ i i = 1 m λ i x i when x i 0 , λ i 0 , i = 1 m λ i = 1 , the relationships between the aggregated values obtained by the PFPCA, PFPCG, GPFPCA, and GPFPCG operators are shown as follows:
Theorem 22.
Let p i = ( μ i , η i , v i ) ( i = 1 , 2 , , m ) be a collection of PFNs, then
P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) ,
G P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) G P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) ,
P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) G P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) ,
G P F P C G F α , β , γ n ( p 1 , p 2 , , p m ) P F P C A F α , β , γ n ( p 1 , p 2 , , p m ) .
Thus, we can conclude that the values obtained by the PFPCG operator are not bigger than the ones obtained by the PFPCA and GPFPCA. The values obtained by the GPFPCG operator are not bigger than the ones obtained by the PFPCA and GPFPCA operators for any value of λ i 0 . Therefore, decision makers can select the four different operators according to their preferences and actual needs.

5. A New Method to Multiattribute Decision-Making with Picture Fuzzy Information

In the present section, we introduce a novel approach to MADM under the picture fuzzy environment. A typical MADM problem with picture fuzzy information can be described as: let X = { x 1 , x 2 , , x m } be a set of alternatives, and G = { G 1 , G 2 , , G s } be a set of attributes. Decision makers are organized to make decisions over alternatives. For attribute G j ( j = 1 , 2 , , s ) of alternative x i ( i = 1 , 2 , , m ) , decision makers are required to use a PFN to express their preference information, which can be denoted as p i j = ( μ i j , η i j , v i j ) ( i = 1 , 2 , , m ; j = 1 , 2 , , s ) . Therefore, a picture fuzzy decision matrix can be obtained P = ( p i j ) m × s . In the following, based on the picture fuzzy aggregation operators, a novel approach to solve this problem is introduced.
Step 1. Generally, there are two kinds of attributes: benefit attributes and cost attributes. Therefore, the decision matrix should be normalized in the decision matrix by
p i j = { ( μ i j , η i j , v i j ) q G j I 1 ( v i j , η i j , μ i j ) q G j I 2 ,
where I1 represents benefit attributes and I2 represents cost attributes. Then a normalized decision matrix can be obtained.
Step 2. Identify the fuzzy measure of the attributes of G.
Step 3. Rearrange the PFNs in a descending order based on the score function S p by Definition 9 or accuracy function H p by Definition 10.
Step 4. For alternative x i ( i = 1 , 2 , , m ) , utilize the series of PFPCA operators, or the series of PFPCG operators, or the series of GPFPCA operators, or the series of GPFPCG operators to aggregate all the attributes values. Therefore, we can get overall values p i ( i = 1 , 2 , , m ) of alternatives.
Step 5. Calculate scores of p i ( i = 1 , 2 , , m ) .
Step 6. Rank alternatives x i ( i = 1 , 2 , , m ) according to the rank of the corresponding overall values.

6. Applications in Supporting the Hierarchical Medical Treatment System with the Proposed Approach

Air pollution is currently the principal issue in the field of environmental health and PM2.5 (fine particulate matter with a aerodynamic diameter of less than 2.5 μm) has become the most important air contaminant in most cities of China, increasing health risks to the Chinese population with respect to respiratory and lung system diseases. Abundant evidence has revealed that exposure to particulate matter air pollution increases the risk of lung cancer since particulate matter with hazardous substances can enter the human body through the respiratory system and is deposited into the lung, giving rise to the damage of pulmonary function. As the research results published in The Lancet [34,35] show, among the risk factors affecting the burden of disease, ambient air pollution rank fourth as risk factor that contributed most to disability adjusted life-years. Air pollution has caused long-lasting adverse effects on respiratory health, and the adverse effect represents a substantial burden with regard to disease prevention and management.
Today, the number of patients with lung disease is soaring due to the above air pollution. In addition, China’s grassroots medical services still leave much to be desired, which to some extent force patients with lung disease to rush to large hospitals even if they only have simple lung health issues. However, the number of patients exceeds the load capacity of the large hospitals, causing great pressure. Under such circumstances, the concept of a hierarchical medical treatment system in accordance with China’s actual conditions by 2020 was introduced into the 13th Five Year Plan (2016–2020). Through the hierarchical medical treatment system, patients with different conditions can choose to go different levels of hospitals instead of all patients rushing to grade III, class A hospitals. Essentially, classifying the different degrees of diseases is a key step in pushing forward the hierarchical system. Therefore, in the present case analysis, we focus on classifying the different degrees of lung diseases to support the hierarchical medical system in China. The specific statement about the medical diagnosis problem is described as follows:
Suppose four patients, denoted by x i (i = 1, 2, 3, 4), who are possibly infected with lung diseases, need to be diagnosed and distributed according to hierarchical medical treatment system. The four patients are diagnosed from the following four symptoms (attributes) of the lung diseases: G1: vital signs, including heart rate, blood pressure, and so on; G2: body temperature (shivering and hyperthermia are two classical symptoms of pneumonia); G3: the frequency of cough; and G4: the frequency of hemoptysis. We invited a doctor who is an expert in lung diseases from a large central hospital. Then, the judgments provided by the doctor for the four patients with respect to the symptoms were represented by PFNs and the decision matrix is shown in Table 1.
With the above four diagnostic criteria for the lung diseases, the patient’s condition can be judged by the doctor. According to the degree and urgency of lung diseases, patients can be distributed to different levels and types of hospitals. Patients with severe conditions should be treated in grade III, class A hospitals, and patients with less severe symptoms should be treated in grade II hospitals. Other common illnesses can be treated in local hospitals. As mentioned in Section 2, the proposed new decision-making method does not only control the certainty of doctor’s decision data, but also deals with these situations where the decision data are correlative. Thus, the new decision-making method is suitable to be employed here.

6.1. Decision-Making Process

(1) The decision-making steps based on the series of PFPCA operators
Step 1. As all the attributes (symptoms) are benefit attributes, the decision matrix does not need to be normalized.
Step 2. Identify the fuzzy measure of the n attributes of G. Suppose that the fuzzy measures of attributes of G are given as follows:
ρ ( G 1 ) = 0.2 , ρ ( G 2 ) = 0.3 , ρ ( G 3 ) = 0.2 , ρ ( G 4 ) = 0.4
The ρ -fuzzy measure is used to calculate the fuzzy measure of attribute sets. Firstly, according to Equation (19), the value of σ is obtained: σ = 0.237 , and then the fuzzy measures of attribute sets of G = { G 1 , G 2 , G 3 , G 4 } can be calculated by Equation (13), shown as follows:
ρ ( G 1 , G 2 ) = 0.486 , ρ ( G 1 , G 3 ) = 0.400 , ρ ( G 1 , G 4 ) = 0.580 ρ ( G 2 , G 4 ) = 0.680 , ρ ( G 3 , G 4 ) = 0.581 , ρ ( G 2 , G 3 ) = 0.486 ρ ( G 1 , G 2 , G 3 ) = 0.663 , ρ ( G 1 , G 2 , G 4 ) = 0.840 , ρ ( G 1 , G 3 , G 4 ) = 0.754 ρ ( G 2 , G 3 , G 4 ) = 0.840 , ρ ( G 1 , G 2 , G 3 , G 4 ) = 1 .
Step 3. According to Table 1, by score functions, rearrange the PFNs in descending order, shown as follows:
p 1 σ ( 1 ) = ( 0.5 ,   0.3 , 0.1 ) , p 1 σ ( 2 ) = ( 0.5 ,   0.1 ,   0.3 ) , p 1 σ ( 3 ) = ( 0.2 ,   0.3 ,   0.4 ) , p 1 σ ( 4 ) = ( 0.6 ,   0.1 , 0.2 ) ,
p 2 σ ( 1 ) = ( 0.6 ,   0.3 ,   0.1 ) , p 2 σ ( 2 ) = ( 0.5 ,   0.2 ,   0.2 ) , p 2 σ ( 3 ) = ( 0.4 , 0.4 ,   0.1 ) , p 2 σ ( 4 ) = ( 0.7 ,   0.1 ,   0.2 ) ,
p 3 σ ( 1 ) = ( 0.6 ,   0.2 ,   0.1 ) , p 3 σ ( 2 ) = ( 0.4 ,   0.1 ,   0.3 ) , p 3 σ ( 3 ) = ( 0.2 ,   0.2 , 0.3 ) , p 3 σ ( 4 ) = ( 0.4 ,   0.3 ,   0.3 ) ,
p 4 σ ( 1 ) = ( 0.2 ,   0.3 ,   0.2 ) , p 4 σ ( 2 ) = ( 0.6 ,   0.1 ,   0.3 ) , p 4 σ ( 3 ) = ( 0.1 ,   0.3 ,   0.5 ) , p 4 σ ( 4 ) = ( 0.1 ,   0.2 , 0.6 ) .
Then we can get
A 1 σ ( 1 ) = { G 2 } , A 1 σ ( 2 ) = { G 2 , G 3 } , A 1 σ ( 3 ) = { G 2 , G 3 , G 4 } , A 1 σ ( 4 ) = { G 1 , G 2 , G 3 , G 4 } ,
A 2 σ ( 1 ) = { G 2 } , A 2 σ ( 2 ) = { G 2 , G 3 } , A 2 σ ( 3 ) = { G 1 , G 2 , G 3 } , A 2 σ ( 4 ) = { G 1 , G 2 , G 3 , G 4 } ,
A 3 σ ( 1 ) = { G 2 } , A 2 σ ( 2 ) = { G 2 , G 3 } , A 2 σ ( 3 ) = { G 1 , G 2 , G 3 } , A 3 σ ( 4 ) = { G 1 , G 2 , G 3 , G 4 } ,
A 4 σ ( 1 ) = { G 4 } , A 4 σ ( 2 ) = { G 1 , G 4 } , A 4 σ ( 3 ) = { G 1 , G 3 , G 4 } , A 4 σ ( 4 ) = { G 1 , G 2 , G 3 , G 4 } .
Taking patientas x1 an example,
ρ 1 A σ ( 1 ) = ρ G 2 = 0.3 , ρ 1 A σ ( 2 ) ρ 1 A σ ( 1 ) = ρ G 2 G 3 ρ G 2 = 0.186 ,
ρ 1 A σ ( 3 ) ρ 1 A σ ( 2 ) = ρ G 2 G 3 G 4 ρ G 2 G 3 = 0.354 , ρ 1 A σ ( 4 ) ρ 1 A σ ( 3 ) = ρ G 1 G 2 G 3 G 4 ρ G 2 G 3 G 4 = 0.16
.
Thus, matrix of fuzzy measure is shown in Table 2:
Step 4. (Suppose n = 1 ) For patients x i ( i = 1 , 2 , 3 , 4 ) , utilize the series of P F P C A operators to aggregate the all the attributes values. Therefore, we can get overall values p i of patients.
Without loss of generality, utilizing the P F P C A F α , β , γ n operator in Theorem 6 to aggregate we get
p 1 = ( 0.547 , 0.201 , 0.281 ) ,   p 2 = ( 0.592 , 0.227 , 0.168 ) ,
p 3 = ( 0.462 , 0.283 , 0.314 ) ,   p 4 = ( 0.411 , 0.212 , 0.398 ) .
Step 5. Calculating scores of p i ( i = 1 , 2 , 3 , 4 ) according to Definition 3, we can get
s ( p 1 ) = 0.266 , s ( p 2 ) = 0.424 , s ( p 3 ) = 0.148 , s ( p 4 ) = 0.013 .
Therefore, the rank of the overall values is p 2 p 1 p 3 p 4 .
Step 6. The rank about the patients’ conditions can be obtained according to the rank of p i , which is x 2 x 1 x 3 x 4 . Further, the patients can be classified according to the ranking.
Therefore, the condition of patient 2 is the most serious, which means that patient 2 should be treated in a grade III, class A hospital. Meanwhile, patient 4 should be treated in a local hospital since his condition is not so serious. On the basis of the availability of the ward, patients 1 and 3 can be referred to other different types of hospitals.
(2) The decision-making steps based on the series of GPFPCA operators
Step 1. As all the attributes are benefit attributes, the decision matrix does not need to be normalized.
Step 2. (Suppose n = 1, λ = 5) For patient xi (i = 1,2,3,4), utilize the series of GPFPCA operators to aggregate the all the attributes values. Therefore, we can get overall values p i of patients.
Without loss of generality, we utilize G P F P C A F α , β , γ n in Theorem 16 to aggregate, and get
p 1 = ( 0.575 , 0.443 , 0.265 ) , p 2 = ( 0.615 , 0.537 , 0.102 ) ,  
p 3 = ( 0.506 , 0.415 , 0.289 ) , p 4 = ( 0.506 , 0.265 , 0.366 ) .
Step 3. Calculating scores of pi (i = 1,2,3,4) according to Definition 9, we can get
s ( p 1 ) = 0.311 , s ( p 2 ) = 0.513 , s ( p 3 ) = 0.217 , s ( p 4 ) = 0.140 .
Therefore, the rank of the overall values is p 2 p 1 p 3 p 4 .
Step 4. The rank about the patients’ conditions can be obtained according to the rank of p i , which is x 2 x 1 x 3 x 4 . Further, the patients can be classified according to the ranking.
Obviously, the ranking results of the two methods proposed in this paper are same. Thus, patient 2 should be treated in a grade III, class A hospital and patient 4 should be treated in a local hospital.

6.2. The Influence of the Parameter Vector λ on the Final Result

The prominent characteristic of our proposed operators is not only efficiently control the certainty degree of PFS given by doctor but also model the medical diagnosis problem in a more flexible manner using an additional parameter. To reflect the influences of different values of parameter λ on the results, we utilize different parameter values of λ to rank the patients by the proposed G P F P C A F α , β , γ n operator, and results are shown in Table 3.
From Table 3, we can know that performance of the patients counts on the values of the parameters λ and the score values may be different for different parameters of λ in the G P F P C A F α , β , γ n operator. However, the ranking result for patients’ conditions is always x 2 x 1 x 3 x 4 . By further analysis, we can easily find that the score values by the G P F P C A F α , β , γ n operator become bigger and bigger with increasing values of λ .
Figure 1 illustrates the scores of the patients’ conditions by the G P F P C A F α , β , γ n operator as assigned different values. From Figure 1, we can find that the scores of patient 3 are very close to those of patient 1 with increasing values of λ . It is easily find that that condition of patient 2 is always the most serious though the ranking results may be different for different parameters. Moreover, as the values of λ become greater and greater, the score values of the patients’ conditions are very close to the fixed values whatever the value of λ is.

6.3. Comparative Analysis

In order to further illustrate the effectiveness and advantages of the proposed method in this paper, we can evaluate the performance of our methods and compare it with the existing methods.

6.3.1. Validity Test

We can employ the following testing criteria [36] to evaluate the performance of our methods.
Test criterion 1: The most serious patient will not change when one of the serious patients is replaced with another slightly serious patient without changing the relative importance of each decision criteria;
Test criterion 2: An effective MCDM method should follow a transitive property;
Test criterion 3: A combined ranking of the patients should be the same as the ranking of original problem when we decompose an MCDM problem into smaller problems.
Thus, in order to test the stability of the patients’ ranking under test criterion 1, we use a new decision matrix where diagnostic values of patients p 1 , p 3 , and p 4 are replaced by less serious ones denoted as p ˜ 1 , p ˜ 3 , and p ˜ 4 respectively, which are determined by subtracting 0.1 from the original positive degree and adding 0.1 to negative degree of patients p 1 , p 3 , and p 4 , respectively. The new decision matrix now is shown in Table 4:
We get
p ˜ 1 = ( 0.448 , 0.201 , 0.386 ) , p ˜ 2 = ( 0.592 , 0.227 , 0.168 ) ,
p ˜ 3 = ( 0.359 , 0.283 , 0.419 ) , p ˜ 4 = ( 0.319 , 0.213 , 0.489 )
and s ( p ˜ 1 ) = 0.062 , s ( p ˜ 2 ) = 0.424 , s ( p ˜ 3 ) = 0.06 , s ( p ˜ 4 ) = 0.169 .
According to the scores of p ˜ i , we can get x 2 x ˜ 1 x ˜ 3 x ˜ 4 . Therefore, patient 2 is in the most serious condition again, which means that our method does not change the indication of the most serious patient when non-severe patients are replaced by other more serious patients.
On the other hand, we can also add a new alternative p 5 that is in less serious condition than p 2 to test criterion 1. The alternative p 5 is obtained by subtracting 0.1 from the original positive degree and adding 0.1 to negative degree of patient p 2 . Then we use the proposed method in Theorem 6 to rank p i (i = 1, 2, 3, 4, 5). We get p 5 = ( 0.271 , 0.328 , 0.232 ) and s ( p 5 ) = 0.039 . According to the scores of p 5 , we can get x 2 x 1 x 3 x 5 x 4 , which means that the indication of patient remains unchanged when adding a new alternative.
Hence, the proposed approach is effective under test criterion 1.
In order to test the validity of our method under criterion 2 and test criterion 3, we decompose the original problem into three subproblems { p 1 , p 2 , p 3 } , { p 2 , p 3 , p 4 } and { p 1 , p 3 , p 4 } . Then we use same approach to solve subproblems and get the ranking corresponding to each subproblem is x 2 x 1 x 3 , x 2 x 3 x 4 and x 1 x 3 x 4 respectively. Then a combined ranking of the patients’ condition is x 2 x 1 x 3 x 4 which is identical to the original overall ranking of the un-decomposed problem. Therefore, the proposed method is valid under test criterion 2 and test criterion 3.

6.3.2. The Advantages of the Proposed Method

To verify the effectiveness and advantages of our methods, we solve the same illustrative example by using different MADM methods including the PFWA operator, the picture fuzzy hybrid averaging (PFHA) operator, the picture fuzzy Einstein weighted average (PFEWA) operator in [13], and the picture fuzzy weighted geometric (PFWG) and picture fuzzy hybrid geometric (PFHG) operators in [12]. The ranking results are shown in Table 5.
From Table 5, we find that there are two final rankings obtained for the different eight operators; seven of them (PFHA, PFWG, PFHG, PFEWA, P F P C A F α , β , γ n , P F P C G F α , β , γ n , G P F P C A F α , β , γ n ) produce the same ranking result, i.e., x 2 x 1 x 3 x 4 and two of them (PFWA and G P F P C G F α , β , γ n ) produce a different ranking result, i.e., x 2 x 3 x 1 x 4 . Through the ranking results may be slightly different by the above MADM methods, the condition of patient 2 is always the most serious and thus should be treated in a grade III, class A hospital, and patient 4 should be treated in a local hospital. This fact verifies that the new method we proposed is effective.
In the following, we compare our proposed operators with the existing operators, such as PFWA, the picture fuzzy order weighted average (PFOWA) operator, PFHA, PFEWA, picture fuzzy Hammer weighted averaging (PFHWA) operator in [13], and PFWG, the picture fuzzy order weighted geometric (PFOWG) operator, PFHG, in [12], and the results are shown in Table 6.
From Table 6, we can find our operators have the following superiorities compared with the existing operators introduced in [12,13]:
(1)
Compared with PFWA and PFWG in [12,13], the new proposed methods are more general and flexible than those provided in the existing literature.
(2)
It should also be noted that the methods introduced in [12,13] are only based on the original information, and thus cannot control the certainty degree, while the new proposed methods can redistribute the membership or non-membership in PFNs according to different principles and thus can get more intensive information from the original PFS.
(3)
From Table 6, it can be concluded that aggregation operators introduced in [12,13] cannot consider correlations among arguments, but the proposed aggregation operators can efficiently take the various interactions among the decision criteria into account. Furthermore, when changing the parameter λ, different scores are acquired shown as in Table 3, which makes decision making more flexible and can meet the needs of different types of decision makers.
Based on the comparisons and analysis above, the methods proposed in this paper can not only control the certainty of the decision data, but also deal with situations where the decision data are correlative. Thus, the proposed approaches in this paper are superior compared with other methods.

7. Conclusions

The hierarchical medical treatment system is an efficient way to integrate all levels of medical service system resources and release the pressure of large hospitals in China. In order to divide patients under different conditions into different levels of hospitals in the hierarchical medical treatment system, doctors need to make scientific assessments of patient’s condition based on their personal experience. In this paper, we proposed a framework to the MADM problem under the picture fuzzy environment. The analytical results show that the new approach is applicable and operational in medical diagnosis, which provides convenience for supporting the hierarchical medical treatment system in China.
There are several main contributions of this work. Firstly, it defined some point operators under picture fuzzy environment to reduce the uncertainty of doctor’s diagnosis data. By the point operators, the PFNs are translated into other PFNs, which can express more intensive information. Secondly, this paper proposed a new class of picture fuzzy point–Choquet operators including the PFPCA, PFPCG, GPFPCA, and GPFPCG operators, which can not only reduce the uncertainty of doctor’s diagnosis data, but also deal with these situations where the decision criteria (symptoms) are correlative. Thirdly, it provided a novel approach to divide patients under different conditions into different levels of hospitals based on the developed operators. Fourthly, the proposed approach has a wide range of applications and can be further applied to other MADM problems, such as pilot hospital selection, supplier selection, emergency decision making, and so on.
Although the proposed model serves better from both theoretical and practical perspectives, this research still has some deficiencies. For instance, there are four parameters in the proposed model and different parameters may lead to different results, while this paper only discusses the influences of different values of parameter λ on results. Thus, an interesting topic worthy of further study in the future is to determine the optimal parameter combination of the proposed operators. Patient condition partially depends on doctor’s judgment, which may be a little subjective when making decisions. Thus, in future study, we will consider utilizing probability theory, such as the probabilistic Pythagorean fuzzy set [37], the proportional hesitant fuzzy set [38], and the probabilistic interval-valued intuitionistic hesitant fuzzy set [39] to express doctor judgments.

Author Contributions

R.Z. provided overall instructions of manuscript drafting and revised the manuscript. Y.X. proposed the idea and drafted the initial manuscript. X.S. and J.W. did the relevant literature review, and revised the manuscript. X.Z. finally checked the manuscript. All authors read and approved the final manuscript.

Funding

This work was partially supported by a Key Project of National Natural Science Foundation of China (NSFC) with grant number 71532002 and 61702023, Fundamental Research Funds for the Central Universities (No. 2017YJS067).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Proof of Theorem 2

Proof. 
We prove Equation (28) holds for all n, and the others can be proved analogously.
From Definition 5, we get F α , β , γ ( p ) , transform a PFN into another PFN, and
π F α , β , γ n = 1 ( μ p + α π p ) ( η p + β π p ) ( v p + γ π p ) = ( 1 μ p η p v p ) ( 1 α β γ ) = π p ( 1 α β γ )
Since F α , β , γ 2 ( p ) = F α , β , γ ( F α , β , γ ( p ) ) , we have
F α , β , γ 2 ( p ) = F α , β , γ ( F α , β , γ ( p ) ) = { ( μ p + α π p + α π F α , β , γ ) , ( η p + β π p + β π F α , β , γ ) , ( v p + γ π p + γ π F α , β , γ ) } .
Thus
π F ξ , ζ 2 ( γ ) = ( 1 [ μ p + α π p + α π F α , β , γ ] [ η p + β π γ q + β π F α , β , γ ] [ v p + γ π γ q + γ π F α , β , γ ] ) = ( ( 1 μ p η p v p ) ( α + β + γ ) π p ( α + β + γ ) ( 1 α β γ ) π p ) = π p ( 1 α β γ ) 2
.
Since F α , β , γ 3 ( p ) = F α , β , γ ( F α , β , γ 2 ( p ) ) , we have
F α , β , γ 3 ( p ) = F α , β , γ ( F α , β , γ 2 ( p ) ) = { ( μ p + α π p + α ( 1 α β γ ) π p + ) + α ( 1 α β γ ) 2 π p , ( η p + β π p + β ( 1 α β γ ) π p ) + β ( 1 α β γ ) 2 π p , ( v p + γ π p + γ ( 1 α β γ ) π p ) + γ ( 1 α β γ ) 2 π p }
and π F α , β , γ 3 ( p ) = π p ( 1 α β γ ) 3 .
Similarly, we have
F α , β , γ n ( p ) = F α , β , γ ( F α , β , γ n 1 ( p ) ) = { ( ( μ p + α π p + ( 1 α β γ ) α π p ) + ( 1 α β γ ) 2 α π p + + ( 1 α β γ ) n 1 α π p ) , ( η p + β π p + ( 1 α β γ ) β π p ) + ( 1 α β γ ) 2 ϕ π p + + ( 1 α β γ ) n 1 β π p , ( ( v p + γ π p + ( 1 α β γ ) γ π p ) + ( 1 α β γ ) 2 γ π p + + ( 1 α β γ ) n 1 γ π p ) } = { μ p + α π p 1 ( 1 α β γ ) n α + β + γ , η p + β π p 1 ( 1 α β γ ) n α + β + γ , v p + γ π p 1 ( 1 α β γ ) n α + β + γ } .
.
Therefore, Equation (28) holds, which completes the proof. ☐

Appendix B. Proof of Theorem 6

Proof. 
We prove the Equation (47) holds for all m, and the others can be proved analogously.
(i) We first prove P F P C A F α , β , γ n is also a PFN.
As 0 μ i , η i , v i 1 , 0 μ i + η i + v i 1 , and
μ F α i , β i , γ i n ( p σ ( i ) ) = μ p σ ( i ) + α i π p σ ( i ) 1 ( 1 α i β i γ i ) n α i + β i + γ i ,  
η F α i , β i , γ i n ( p i ) = η p σ ( i ) + β i π p σ ( i ) 1 ( 1 α i β i γ i ) n α i + β i + γ i ,
v F α i , β i , γ i n ( p i ) = v p σ ( i ) + γ i π p σ ( i ) 1 ( 1 α i β i γ i ) n α i + β i + γ i .
Thus 0 i = 1 m ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i 1 , we get
0 ( 1 i = 1 m ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i ) 1   and   0 i = 1 m v F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i 1
Again
1 i = 1 m ( 1 μ F α i , β i , γ i n ( p i ) ) ω ˜ i + i = 1 m η F α i , β i , γ i n ( p i ) ω ˜ i + i = 1 m v F α i , β i , γ i n ( p i ) ω ˜ i 1 i = 1 m ( 1 η F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i ν F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i ) ω ˜ i + i = 1 m η F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i + i = 1 m ν F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i = 1
.
Hence, P F P C A F α , β , γ n is a PFN according to Definition 5.
(ii) Next, we prove Equation (47) by using mathematical induction on m .
When m = 2 , by the operational law (18) and (19) in Section 2.2, we have
P F P C A F α , β , γ n ( p 1 , p 2 ) = ω ˜ 1 p 1 ω ˜ 2 p 2 = ( ( 1 ( 1 μ F α 1 , β 1 γ 1 n ( p 1 ) ) ω ˜ 1 ) , η F α 1 , β 1 , γ 1 n ( p 1 ) ω ˜ 1 , v F α 1 , β 1 γ 1 n ( p 1 ) ω ˜ 1 ) ( ( 1 ( 1 μ F α 2 , β 2 , γ 2 n ( p 2 ) ) ω ˜ 2 ) , η F α 2 , β 2 , γ 2 n ( p 2 ) ω ˜ 1 , v F α 2 , β 2 , γ 2 n ( p 2 ) ω ˜ 2 ) = { ( 1 ( 1 μ F α 1 , β 1 γ 1 n ( p 1 ) ) ω ˜ 1 + 1 ( 1 μ F α 2 , β 2 , γ 2 n ( p 2 ) ) ω ˜ 2 ( 1 ( 1 μ F α 1 , β 1 γ 1 n ( p 1 ) ) ω ˜ 1 ) ( 1 ( 1 μ F α 2 , β 2 , γ 2 n ( p 2 ) ) ω ˜ 2 ) ) , η F α 1 , β 1 , γ 1 n ( p 1 ) ω ˜ 1 η F α 2 , β 2 , γ 2 n ( p 2 ) ω ˜ 1 , v F α 1 , β 1 , γ 1 n ( p 1 ) ω ˜ 1 v F α 2 , β 2 , γ 2 n ( p 2 ) ω ˜ 2 } = ( 1 ( 1 μ F α 1 , β 1 γ 1 n ( p 1 ) ) ω ˜ 1 ( 1 μ F α 2 , β 2 , γ 2 n ( p 2 ) ) ω ˜ 2 , η F α 1 , β 1 , γ 1 n ( p 1 ) ω ˜ 1 η F α 2 , β 2 , γ 2 n ( p 2 ) ω ˜ 1 , v F α 1 , β 1 , γ 1 n ( p 1 ) ω ˜ 1 v F α 2 , β 2 , γ 2 n ( p 2 ) ω ˜ 2 } ) = ( 1 i = 1 2 ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i , i = 1 2 η F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i , i = 1 2 ν F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i ) .
Thus, result is true for m = 2 .
If Equation (47) holds for m = k , that is
P F P C A F α , β , γ n ( p 1 , p 2 , , p k ) = ( 1 i = 1 k ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i , i = 1 k η F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i , i = 1 k v F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i ) ,
then, when m = k + 1 , by the operational laws (7) and (9) in Section 2.2, we have
P F P C A F α , β , γ n ( p 1 , p 2 , , p k + 1 ) = P F P C A F α , β , γ n ( p 1 , p 2 , , p k ) ω ˜ k + 1 F α , β , γ n ( p k + 1 ) = ( 1 i = 1 k ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i , i = 1 k η F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i , i = 1 k v F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i ) ( 1 ( 1 μ F α k + 1 , β k + 1 , γ k + 1 n ( p σ ( k + 1 ) ) ) ω ˜ i , η F α k + 1 , β k + 1 , γ k + 1 n ( p σ ( k + 1 ) ) ω ˜ k + 1 , v F α k + 1 , β k + 1 , γ k + 1 n ( p σ ( k + 1 ) ) ω ˜ k + 1 ) = { ( 1 i = 1 k ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i + 1 ( 1 μ F α k + 1 , β k + 1 , γ k + 1 n ( p k + 1 ) ) ω ˜ k + 1 ( 1 i = 1 k ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i ) ( 1 ( 1 μ F α k + 1 , β k + 1 , γ k + 1 n ( p k + 1 ) ) ω ˜ k + 1 ) ) , i = 1 k η F α i , β i , γ i n ( p i ) ω ˜ i η F α k + 1 , β k + 1 , γ k + 1 n ( p σ ( i ) ) ω ˜ k + 1 , i = 1 k v F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i v F α k + 1 , β k + 1 , γ k + 1 n ( p σ ( i ) ) ω ˜ k + 1 } = ( ( 1 i = 1 k + 1 ( 1 μ F α i , β i , γ i n ( p σ ( i ) ) ) ω ˜ i ) , i = 1 k + 1 η F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i , i = 1 k + 1 v F α i , β i , γ i n ( p σ ( i ) ) ω ˜ i ) ,
Thus, Equation (47) holds for m = k + 1 .
Therefore, Equation (47) holds for all m , which completes the proof. ☐

References

  1. Tang, M.; Liao, H.C.; Li, Z.M.; Xu, Z.S. Nature disaster risk evaluation with a group decision making method based on incomplete hesitant fuzzy linguistic preference relations. Int. J. Environ. Res. Public Health 2018, 15, 751. [Google Scholar] [CrossRef] [PubMed]
  2. Qi, X.W.; Zhang, J.L.; Zhao, S.P.; Liang, C.Y. Tackling complex emergency response solutions evaluation problems in sustainable development by fuzzy group decision making approaches with considering decision hesitancy and prioritization among assessing criteria. Int. J. Environ. Res. Public Health 2017, 14, 1165. [Google Scholar] [CrossRef] [PubMed]
  3. Lin, H.; Zhou, J.W. Linguistic multi-attribute group decision making with risk preferences and its use in low-carbon tourism destination selection. Int. J. Environ. Res. Public Health. 2017, 14, 1078. [Google Scholar] [CrossRef] [PubMed]
  4. Chatterjee, P.; Mondal, S.; Boral, S.; Banerjee, A.; Chakraborty, S. A novel hybrid method for non-traditional machining process selection using factor relationship and Multi-Attributive Border Approximation Method. Facta Univ. Ser. Mech. Eng. J. 2017, 15, 439–456. [Google Scholar] [CrossRef]
  5. Roy, J.; Adhikary, K.; Kar, S.; Pamučar, D. A rough multicriteria approach for evaluation of the supplier criteria in automotive industry. Decis. Mak. Appl. Manag. Eng. 2018, 1, 82–96. [Google Scholar]
  6. Vasiljević, M.; Fazlollahtabar, H.; Stević, Ž.; Vesković, S. A rough strength relational DEMATEL model for analysing the key success factors of hospital service quality. Decis. Mak. Appl. Manag. Eng. 2018, 1, 121–142. [Google Scholar]
  7. Atanassov, K.T. Intuitionistic fuzzy sets. Fuzzy Sets Syst. 1986, 20, 87–96. [Google Scholar] [CrossRef]
  8. Cuong, B.C. Picture fuzzy sets. J. Comput. Sci. Cybern. 2014, 30, 409–420. [Google Scholar]
  9. Wei, G.W.; Alsaadi, F.E.; Hayat, T.; Alsaedi, A. Projection models for multiple attribute decision making with picture fuzzy information. Int. J. Mach. Learn. Cybern. 2018, 9, 713–719. [Google Scholar] [CrossRef]
  10. Wei, G.W. Picture 2-tuple linguistic Bonferroni mean operators and their application to multiple attribute decision making. Int. J. Fuzzy Syst. 2016, 19, 1–14. [Google Scholar] [CrossRef]
  11. Wei, G.W. Picture uncertain linguistic Bonferroni mean operators and their application to multiple attribute decision making. Kybernetes 2017, 46, 1777–1800. [Google Scholar] [CrossRef]
  12. Wang, C.; Zhou, X.; Tu, H. Some geometric aggregation operators based on picture fuzzy sets and their application in multiple attribute decision making. Indian J. Pure Appl. Math. 2017, 37, 477–492. [Google Scholar]
  13. Garg, H. Some picture fuzzy aggregation operators and their applications to multi-criteria decision-making. Arab. J. Sci. Eng. 2017, 42, 1–16. [Google Scholar] [CrossRef]
  14. Abbas, M.; Nilashi, M.; Kazimieras, E.; Rahmah, S.; Zare, H.; Mohd, N. Decision making methods based on fuzzy aggregation operators: Three decades review from 1986 to 2017. Int. J. Inf. Technol. Decis. 2018, 17, 391–466. [Google Scholar]
  15. He, Y.D.; Chen, H.Y.; He, Z.; Zhou, L.G. Multi-attribute decision making based on neutral averaging operators for intuitionistic fuzzy information. Appl. Soft Comput. 2015, 27, 64–76. [Google Scholar] [CrossRef]
  16. Choquet, G. On a class of set-functions. Bull. Am. Math. Soc. 1954, 60, 64. [Google Scholar]
  17. Zhou, L.G.; Chen, H.Y. Continuous generalized OWA operator and its application to decision making. Fuzzy Sets Syst. 2011, 168, 18–34. [Google Scholar] [CrossRef]
  18. Meng, F.Y.; Zhang, Q. Induced continuous Choquet integral operators and their application to group decision making. Comput. Ind. Eng. 2014, 68, 42–53. [Google Scholar] [CrossRef]
  19. Xu, Z.S. Choquet integrals of weighted intuitionistic fuzzy information. Inf. Sci. 2010, 180, 726–736. [Google Scholar] [CrossRef]
  20. Yager, R.R. On Using the Shapley Value to Approximate the Choquet Integral in Cases of Uncertain Arguments. IEEE Trans. Fuzzy Syst. 2018, 26, 1303–1310. [Google Scholar] [CrossRef]
  21. Liu, B.; Fu, M.; Zhang, S.; Xue, B.; Zhou, Q.; Zhang, S. An interval-valued 2-tuple linguistic group decision-making model based on the Choquet integral operator. Int. J. Syst. Sci. 2018, 49, 407–424. [Google Scholar] [CrossRef]
  22. Wen, X.Q.; Yan, M.C.; Xian, J.Y.; Yue, R.; Peng, A.H. Supplier selection in supplier chain management using Choquet integral-based linguistic operators under fuzzy heterogeneous environment. Fuzzy Optim. Decis. Mak. 2016, 15, 307–330. [Google Scholar] [CrossRef]
  23. Qu, G.H.; Li, Y.J.; Qu, W.H.; Li, C.H. Some new Shapley dual hesitant fuzzy Choquet aggregattion operators and their applications to multiple attribute group decision making-based TOPSIS. J. Int. Fuzzy Syst. 2017, 33, 2463–2483. [Google Scholar]
  24. Joshi, D.; Kumar, S. Interval-valued intuitionistic hesitant fuzzy Choquet integral based TOPSIS method for multi-criteria group decision making. Eur. J. Oper. Res. 2016, 248, 183–191. [Google Scholar] [CrossRef]
  25. Peng, X.D.; Yang, Y. Pythagorean fuzzy Choquet integral based MABAC method for multiple attribute group decision making. Int. J. Intell. Syst. 2016, 31, 989–1020. [Google Scholar] [CrossRef]
  26. Atanassov, K.T. Intuitionistic Fuzzy Sets; Physica-Verlag: Heidelberg, Germany; New York, NY, USA, 1999. [Google Scholar]
  27. Liu, H.W.; Wang, G.J. Multi-criteria decision-making methods based on intuitionistic fuzzy sets. Eur. J. Oper. Res. 2007, 179, 220–233. [Google Scholar] [CrossRef]
  28. Xia, M.M.; Xu, Z.S. Generalized point operators for aggregating intuitionistic fuzzy information. Int. J. Intell. Syst. 2010, 25, 1061–1080. [Google Scholar] [CrossRef]
  29. Xia, M.M. Point operators for intuitionistic multiplicative information. J. Intell. Fuzzy Syst. 2016, 28, 615–620. [Google Scholar]
  30. Peng, X.D.; Yang, Y. Fundamental properties of Pythagorean fuzzy aggregation operators. Int. J. Intell. Syst. 2016, 147, 415–446. [Google Scholar] [CrossRef]
  31. Peng, X.D.; Yuan, H.Y. Fundamental properties of interval-valued Pythagorean fuzzy aggregation operators. Int. J. Intell. Syst. 2016, 31, 444–487. [Google Scholar] [CrossRef]
  32. Xing, Y.P.; Zhang, R.T.; Xia, M.M.; Wang, J. Generalized point aggregation operators for dual hesitant fuzzy information. J. Intell. Fuzzy Syst. 2017, 33, 515–527. [Google Scholar] [CrossRef]
  33. Grabisch, M. Fuzzy integral in multi-criteria decision making. Fuzzy Sets Syst. 1995, 69, 279–298. [Google Scholar] [CrossRef]
  34. Guan, W.J.; Zheng, X.Y.; Chung, K.F.; Zhong, N.F. Impact of air pollution on the burden of chronic respiratory diseases in China: time for urgent action. Lancet 2016, 388, 1939–1951. [Google Scholar] [CrossRef]
  35. Yang, G.; Wang, Y.; Zeng, Y.; Gao, G.F.; Liang, X.; Zhou, M.; Wan, X.; Yu, S.; Jiang, Y.; Naghavi, M.; et al. Rapid health transition in China,1990–2010: findings from the Global Burden of Disease Study 2010. Lancet. 2013, 381, 1987–2015. [Google Scholar] [CrossRef]
  36. Wang, X.; Triantaphyllou, E. Ranking irregularities when evaluating alternatives by using some ELECTRE methods. Omega 2008, 36, 45–63. [Google Scholar] [CrossRef]
  37. Zeng, S. Pythagorean fuzzy multi-attribute group decision making with probabilistic information and OWA approach. Int. J. Intell. Syst. 2017, 32, 1136–1150. [Google Scholar] [CrossRef]
  38. Xiong, S.H.; Chen, Z.S.; Chin, K.S. A novel MAGDM approach with proportional hesitant fuzzy sets. Int. J. Comput. Intell. Syst. 2018, 11, 256–271. [Google Scholar] [CrossRef]
  39. Zhai, Y.L.; Xu, Z.S.; Liao, H.C. Measures of probabilistic interval-valued intuitionistic hesitant fuzzy sets and the application in reducing excessive medical examinations. IEEE Trans. Fuzzy Syst. 2018, 26, 1651–1671. [Google Scholar] [CrossRef]
Figure 1. Score values when λ [ 5 , 50 ] by the GPFPCAF α , β , γ n operator.
Figure 1. Score values when λ [ 5 , 50 ] by the GPFPCAF α , β , γ n operator.
Ijerph 15 01718 g001
Table 1. The picture fuzzy decision matrix.
Table 1. The picture fuzzy decision matrix.
G1G2G3G4
x1(0.6, 0.1, 0.2)(0.5, 0.3, 0.1)(0.5, 0.1, 0.3)(0.2, 0.3, 0.4)
x2(0.4, 0.4, 0.1)(0.6, 0.3, 0.1)(0.5, 0.2, 0.2)(0.7, 0.1, 0.2)
x3(0.2, 0.2, 0.3)(0.6, 0.2, 0.1)(0.4, 0.1, 0.3)(0.4, 0.3, 0.3)
x4(0.6, 0.1, 0.3)(0.1, 0.2, 0.6)(0.1, 0.3, 0.5)(0.2, 0.3, 0.2)
Table 2. Matrix of fuzzy measure.
Table 2. Matrix of fuzzy measure.
ρ A σ ( 1 ) ρ A σ ( 2 ) ρ A σ ( 1 ) ρ A σ ( 3 ) ρ A σ ( 2 ) ρ A σ ( 4 ) ρ A σ ( 3 )
x10.30.1860.3540.16
x20.30.1860.1770.337
x30.30.1860.1770.337
x40.40.1860.1740.246
Table 3. Rankings with different value of parameter.
Table 3. Rankings with different value of parameter.
Parameters Score   Value   of   p i   ( i   =   1 , 2 , 3 , 4 ) Ranking Results
λ = 5 s ( p 1 ) = 0.311 s ( p 2 ) = 0.513 s ( p 3 ) = 0.217 s ( p 4 ) = 0.14 x 2 x 1 x 3 x 4
λ = 10 s ( p 1 ) = 0.350 s ( p 2 ) = 0.542 s ( p 3 ) = 0.297 s ( p 4 ) = 0.204 x 2 x 1 x 3 x 4
λ = 15 s ( p 1 ) = 0.379 s ( p 2 ) = 0.562 s ( p 3 ) = 0.348 s ( p 4 ) = 0.231 x 2 x 1 x 3 x 4
λ = 20 s ( p 1 ) = 0.401 s ( p 2 ) = 0.576 s ( p 3 ) = 0.381 s ( p 4 ) = 0.246 x 2 x 1 x 3 x 4
λ = 30 s ( p 1 ) = 0.429 s ( p 2 ) = 0.594 s ( p 3 ) = 0.418 s ( p 4 ) = 0.262 x 2 x 1 x 3 x 4
λ = 40 s ( p 1 ) = 0.446 s ( p 2 ) = 0.606 s ( p 3 ) = 0.438 s ( p 4 ) = 0.271 x 2 x 1 x 3 x 4
λ = 50 s ( p 1 ) = 0.457 s ( p 2 ) = 0.614 s ( p 3 ) = 0.451 s ( p 4 ) = 0.276 x 2 x 1 x 3 x 4
Table 4. The new picture fuzzy decision matrix.
Table 4. The new picture fuzzy decision matrix.
G1G2G3G4
x1(0.5, 0.1,0.1)(0.4, 0.3,0.2)(0.4, 0.1, 0.4)(0.1, 0.3, 0.5)
x2(0.4, 0.4, 0.1)(0.6, 0.3, 0.1)(0.5, 0.2, 0.2)(0.7, 0.1, 0.2)
x3(0.1, 0.2, 0.4)(0.5, 0.2, 0.2)(0.3, 0.1, 0.2)(0.3, 0.3, 0.4)
x4(0.5, 0.1, 0.4)(0.1, 0.2, 0.7)(0.1, 0.3, 0.6)(0.1, 0.3, 0.3)
Table 5. Comparison of rankings with different aggregation operators.
Table 5. Comparison of rankings with different aggregation operators.
ApproachesScore Value of Xi (i = 1, 2, 3, 4)Ranking
Approach based on the PFWA operator [13] s ( p 1 ) = 0.199 s ( p 2 ) = 0.463
s ( p 3 ) = 0.222 s ( p 4 ) = 0.061
x 2 x 1 x 3 x 4
Approach based on the PFHA operator [13] s ( p 1 ) = 0.239 s ( p 2 ) = 0.407
s ( p 3 ) = 0.171 s ( p 4 ) = 0.108
x 2 x 1 x 3 x 4
Approach based on the PFWG operator [12] s ( p 1 ) = 0.243 s ( p 2 ) = 0.402
s ( p 3 ) = 0.103 s ( p 4 ) = 0.061
x 2 x 1 x 3 x 4
Approach based on the PFHG operator [12] s ( p 1 ) = 0.327 s ( p 2 ) = 0.413
s ( p 3 ) = 0.279 s ( p 4 ) = 0.133
x 2 x 1 x 3 x 4
Approach based on the PFEWA operator [13] s ( p 1 ) = 0.261 s ( p 2 ) = 0.425
s ( p 3 ) = 0.139 s ( p 4 ) = 0.005
x 2 x 1 x 3 x 4
Approach   based   on   P F P C A F α , β , γ n
operator (in this paper)
s ( p 1 ) = 0.265 s ( p 2 ) = 0.424
s ( p 3 ) = 0.148 s ( p 4 ) = 0.013
x 2 x 1 x 3 x 4
Approach   based   on   the   P F P C G F α , β , γ n
operator (in this paper)
s ( p 1 ) = 0.204 s ( p 2 ) = 0.399
s ( p 3 ) = 0.095 s ( p 4 ) = 0.051
x 2 x 1 x 3 x 4
Approach   based   on   the   G P F P C A F α , β , γ n
operator (in this paper)
s ( p 1 ) = 0.310 s ( p 2 ) = 0.513
s ( p 3 ) = 0.217 s ( p 4 ) = 0.140
x 2 x 1 x 3 x 4
Approach   based   on   the   G P F P C G F α , β , γ n
operator (in this paper)
s ( p 1 ) = 0.095 s ( p 2 ) = 0.358
s ( p 3 ) = 0.111 s ( p 4 ) = 0.091
x 2 x 3 x 1 x 4
Table 6. The comparison of different operators.
Table 6. The comparison of different operators.
Aggregation OperatorsWhether It Can Consider Correlations among ArgumentsWhether It Can Control the Certainty of PFNsFlexible (Whether There Is a Parameter to Reflect Preferences)
PFWA [13]NoNoNo
PFOWA [13]NoNoNo
PFHA [13]NoNoNo
PFWG [12]NoNoNo
PFOWG [12]NoNoNo
PFHG [12]NoNoNo
PFEWA [13]NoNoNo
PFHA [13]NoNoNo
PFPCAYesYesYes
PFPCGYesYesYes
GPFPCAYesYesYes
GPFPCGYesYesYes

Share and Cite

MDPI and ACS Style

Zhang, R.; Xing, Y.; Wang, J.; Shang, X.; Zhu, X. A Novel Multiattribute Decision-Making Method Based on Point–Choquet Aggregation Operators and Its Application in Supporting the Hierarchical Medical Treatment System in China. Int. J. Environ. Res. Public Health 2018, 15, 1718. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15081718

AMA Style

Zhang R, Xing Y, Wang J, Shang X, Zhu X. A Novel Multiattribute Decision-Making Method Based on Point–Choquet Aggregation Operators and Its Application in Supporting the Hierarchical Medical Treatment System in China. International Journal of Environmental Research and Public Health. 2018; 15(8):1718. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15081718

Chicago/Turabian Style

Zhang, Runtong, Yuping Xing, Jun Wang, Xiaopu Shang, and Xiaomin Zhu. 2018. "A Novel Multiattribute Decision-Making Method Based on Point–Choquet Aggregation Operators and Its Application in Supporting the Hierarchical Medical Treatment System in China" International Journal of Environmental Research and Public Health 15, no. 8: 1718. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph15081718

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

Article Metrics

Back to TopTop