Probabilistic Hesitant Intuitionistic Fuzzy Linguistic Term Sets and Their Application in Multiple Attribute Group Decision Making

.


Introduction
Multi-attribute group decision-making (MAGDM) is usually described as the process of selecting the optimal option that has the highest degree of satisfaction from a series of finite alternatives, which are evaluated by multiple experts over multiple attributes.A complete MAGDM process can be divided into three parts: (1) expression of the evaluation values provided by a group of experts; (2) determining the weights of attributes; (3) aggregating the individual evaluations and obtaining the final ranking [1,2].Over the past few decades, many researchers have investigated the theories and methods of MAGDM problems, and a great deal of research has been done in the fields of operational laws [3,4], aggregation operators [5][6][7], distance measures and similarity measure [8,9], and outranking relations [10][11][12].
In many practical situations, to describe complex and imprecise cognition, it is more convenient to provide the decision-making information in linguistic terms, such as "good" quality or "very bad" credit level, rather than precise numerical values, owing to the unquantifiable and ambiguity of the real-life condition.Zadeh [13] proposed the fuzzy linguistic approach, which adopts linguistic terms to depict the intensity of membership between an element and a set.However, experts show some hesitancy among several evaluation values, as the fuzzy linguistic approach only allows singleton linguistic term to represent the value of a linguistic variable, and thus it may not express the preference of the decision maker accurately.To improve the elicitation of linguistic information, Rodriguez et al. [14] extended the fuzzy linguistic approach and introduced hesitant fuzzy linguistic term sets (HFLTSs), which allow experts to provide their evaluations in several linguistic term values.For example, the evaluation of the new product's market prospect can be "between medium and high", and the evaluation of a supplier's service capability can be "at least good".In recent years, HFLTSs have received a great deal of attention from researchers and theories and its methods have been widely discussed [15][16][17][18][19][20][21].
Nevertheless, experts assume equal weight when employing each linguistic term to describe the preference intensity over HFLTSs, which may not reflect the actual conditions in some cases.For example, an expert can use 1 2 3 ( , , ) s s s to express his opinions and assign a weight to each linguistic term as (0.5,0.3,0.2) , respectively.The HFLTS fails to describe such a situation accurately due to reasons of limitations in expression.Thus, Pang et al. [22] presented probabilistic linguistic term sets (PLTSs), which not only describe the hesitant information with more than one linguistic term, but also reflect the evaluation values by assigning a weight to each linguistic term.Due to the characteristics of more flexibility and comprehensiveness, many researches have investigated the related calculation methods and extended theories.Pang et al. [22] defined the basic operations and aggregation operators for PLTS and applied them to deal with the MAGDM problems.Wu et al. [23] presented an aggregation method to obtain the collective judgement by integrating the individual opinions for PLTS.Zhang et al. [24] introduced the concept of the probabilistic linguistic preference relation (PLPR) and applied it to an investment risk assessment problem.Zhang et al. [25] investigated the consensus reaching process for MAGDM with PLPRs.Liao et al. [26] established the linear programming model of probabilistic linguistics to handle the MADM and gave an application example about evaluating the hospital levels in China.However, PLTSs also have a big limitation for the reason of regarding the membership degree of a linguistic value to be 1 by default, and do not consider the degree of non-membership and hesitancy of a linguistic variable.However, in fact, for the reasons of cognitive limitations and knowledge gaps, this information is indispensable.A simple example of evaluating the sales of a new product can be "60% sure high and 40% sure very high, and it cannot be bad".Therefore, Malik et al. [27] proposed probabilistic hesitant intuitionistic fuzzy linguistic term sets (PHIFLTSs), which facilitated experts to express their judgment about membership and non-membership information in linguistic terms simultaneously to cope with such complicated and ambiguous situations in the realworld.
Nevertheless, since PHIFLTSs have been put forward only a few years ago, many issues need to be solved to enrich its methodology framework.
(1) The existing PHIFLTSs proposed in [27] are based on the direct calculation between subscripts of linguistic terms and their associated probabilities.However, a great deal of important defects for these operations appear in some situations.For example, let s s , but also loses the associated probability information.(2) The existing forms of PHIFLTSs are based on a qualitative scale mapped to a sequence of adjacent integers that are equally distributed.In fact, unbalanced conditions are very common if we take the psychology of experts into account [28].Torra [29] discussed the unbalanced semantics for an ordered set of linguistic terms for the first time.Afterward, many achievements were obtained about the unbalanced linguistic terms [30][31][32][33][34][35][36][37].However, Malik et al. [27] did not consider the situation of unbalanced linguistic terms on the PHIFLTS environment.(3) The distance measure defined in [27] is also based on the subscripts of each linguistic term and the associated probability.However, on the one hand, it cannot handle the conditions of unbalanced linguistic terms over PHIFLTSs; on the other hand, unreasonable results may be produced in some specific situations.(4) Compared with the current well-known ranking techniques including TOPSIS, VIKOR, ELECTRE, PROMETHEE, AHP, and MOORA, the multiplicative multi-objective optimization by ratio analysis (MULTIMOORA), which is specific to peculiarities of three aggregation models from the aspects of fully compensatory, non-compensatory, and incompletely compensatory has separately verified the superiority concerning time consumption, robustness, simplicity, and effectiveness [38][39][40][41].However, advances in the state of the art have shown that no studies applied the MULTIMOORA method to the probabilistic hesitant intuitionistic fuzzy linguistic term sets environment.
To overcome the above-stated limitations and enrich the theories and methods of the PHIFLTSs, this paper focuses on the following four main aspects: (1) This paper proposes some novel operational laws for the probabilistic hesitant intuitionistic fuzzy linguistic term sets to enrich the computation between PHIFLTSs and to improve the applicability and methodology in multi-attribute group decision-making models.(2) This paper establishes the distance and correlation measures for the probabilistic hesitant intuitionistic fuzzy linguistic term sets, which make up for the shortage of the current distance measures.
(3) This paper takes the unbalanced linguistic terms of the probabilistic hesitant intuitionistic fuzzy linguistic term sets environment into account to better describe the difference in cognitive information of experts under different situations.(4) This paper presents a MAGDM model based on the MULTIMOORA approach with the use of the developed novel operational laws and correlation measures.
Based on these four innovative strands of research, the following main goals of this paper are two-fold: (1) with the use of the proposed novel operational laws, distance, and correlation measure, and the unbalanced linguistic terms PHIFLTSs, to provide new theories and methods for the PHIFLTSs to overcome the existing defects; (2) with the use of the MULTIMOORA ranking approach, to propose a new multi-attribute group decision-making method for the PHIFLTSs to obtain superiority concerning time consumption, robustness, simplicity, and effectiveness at the decision making level.
The structure of this paper is as follows: Some concepts about PHIFLTSs, the normalized PHIFLTSs, intuitionistic fuzzy sets, and linguistic scale function are briefly reviewed in Section 2. Section 3 discusses the comparison between PHIFLEs.The novel operational laws for PHIFLTSs are presented in Section 4. Section 5 discusses the distance and correlation measure of the PHIFLTSs.The MAGDM method based on the MULTIMOORA approach under the PHIFLTSs environment are proposed in Section 6.An application example is presented in Section 7 to illustrate the proposed approach.Section 8 presents a comparison with some existing related methods.Our conclusions are presented in Section 9.

Preliminaries
In this section, the basic definitions and operations about probabilistic hesitant intuitionistic fuzzy linguistic term sets, intuitionistic fuzzy sets, correlation measures, and correlation coefficient are introduced.

The Normalization of PHIFLEs
Normally, the sum of the membership probabilities and the non-membership probabilities are both 1 in PHIFLE.However, if either of them has a sum less than 1, the normalization process should be employed to resolve it.Definition 2. [27] Let can be defined as follows:

Intuitionistic Fuzzy Sets and the Basic Operational Laws
Definition 3. [42] Let X be a non-empty universe of discourse and an intuitionistic fuzzy set (IFS) A in X is given by

{ , ( ), ( ) : }
 .The numbers of ( ) A x  and ( ) A v x denote the degree of membership and non-membership of x to the set A, respectively. ( ) IFS X is used to denote the set of all the IFSs in X.In addition, for each The concept of linguistic scale function (LSF) can be utilized to different types of semantics of linguistic terms.Normally, the LSF f is a monotonically increasing function, where : . According to the different semantic distributions represented by linguistic terms, the authors in [43] defined three different types of LSFs, as follows: (1) If the semantics of linguistic terms are divide equally, the LSF (the first type) is defined as follows: ( ) , where 0,1, ..., (2) If the semantics of any adjacent linguistic term are unequally divided, the deviation values in every side of /2 g s both present the tendency of increase; the LSF (the second type) is defined as follows: If the scale level is .
(3) If the semantics of any adjacent linguistic term are unequally divided, the deviation value in every side of /2 g s present the tendency of decrease; the LSF (the third type) is defined as follows: where  , (   are determined by the decision makers based on their psychological attitudes.

The Comparison between PHIFLEs
With the aim of comparing the PHIFLEs, authors in [27] proposed the score function and the deviation degree of PHIFLTSs.However, the value of the score function defined in [27] is a linguistic term, while the score should be a crisp number, and it also cannot handle the unbalanced linguistic terms PHIFLTSs.To overcome these defects, we propose a new score function and variance value function of the PHIFLTSs in this section.
Based on Equations ( 2) and (3), the linguistic scale functions f can be utilized to handle the conditions of unbalanced linguistic terms on the PHIFLTSs environment.For two PHIFLEs For two PHIFLEs 1 ( ) s E P and 2 ( )

( ( )) ( ( ))
S S S E P S E P  , we thus cannot distinguish these two PHIFLEs.Therefore, utilizing Definition 6, we can obtain

Some Novel Operational Laws of the PHIFLEs
As aforementioned, the existing operations over PHIFLEs encounter many important defects due to the reasons of combining two different dimensions (linguistic terms and their associated probabilities) together.To consider eliminating the above-stated limitations, it is necessary to propose some novel operational laws for PHIFLEs.Inspired by [42], the linguistic scale functions f and f −1 can be viewed as a transformations tool, which make the equivalent transformations between the PHIFLEs and the IFVs possible.Utilizing the linguistic scale function f, the PHIFLEs are transformed to the IFVs, and then we can use the operations (Definition 4) to calculate these transformed IFVs.Furthermore, the inverse functions f −1 can be used to transform these calculation results to the PHIFLEs equivalently.Therefore, the novel operational laws for PHIFLEs are defined: ( ) ( ), ( )

S S S E P M p N p 
 be two normalized PHIFLEs, f and f −1 be the linguistic scale functions and the inverse functions, and  be a positive real number.Then (1) be an LTS, ( ) ( ), ( ) and  , 1  , 2  be three positive real numbers.Then (1) (1) and ( 2) are obvious, and the proofs are omitted. (3): ) ( : ) ( ) { (0.5), (0.5)},{ (0. is selected as the equivalent transformation functions f. Then: Utilizing the linguistic scale functions f, we can obtain . Then, based on Definition 7, we can obtain the final results as follows: (1)

Distance and Correlation Measure of the PHIFLTSs Based on Adjusted PHIFLEs
Distance measure can reflect the specific differences between two evaluations or sets.The existing distance measure is based on the subscripts of each linguistic term and the associated probability.However, it cannot handle the unbalanced linguistic terms PHIFLTSs and also may produce unreasonable results in some specific situations.Therefore, inspired by the method proposed in [23], we propose a new distance and correlation measure for PHIFLTSs based on adjusted PHIFLE.

The Adjusted PHIFLEs with Different Probability Distributions
Usually, this is a common situation for two PHIFLEs . In such situations, the existing method needs to increase the linguistic terms (the probabilities associated with the added linguistic terms are "0") in the relatively shorter PHIFLE to equalize their lengths, because the PHIFLEs with different lengths bring important problems in distance measure or aggregation operators.However, on the one hand, multiplying the probabilities by the subscripts is unreasonable; on the other hand, unreasonable results may be produced.
Wu et al. [23] proposed a novel method to adjust the linguistic terms with different probability distributions.With this method, the PLEs can be adjusted to the same probability set and the operations reduced to the dimension of the linguistic terms only.Therefore, we utilized it to adjust the probability distribution on PHIFLEs.Definition 8. [23] Analogously, The processes of adjusting the PHIFLEs are described in Figure 1.

The Distance Measure of the PHIFLTSs Based on the Adjusted PHIFLEs
The authors in [27] presented the distance measure between PHIFLEs, which is based on the aggregation of the subscript of each linguistic term and the associated probability.Nevertheless, on the one hand, it cannot handle the situations of unbalanced PHIFLTSs; on the other hand, it also encounters defects in some specific situations.Let us give an example to reveal it.Proof.
(1) It is obvious and the proof is omitted.
(2) Since , and thus , that is (3) If E P E P  .□

The Correlation Measure and Correlation Coefficients of PHIFLEs
Contrary to the distance measure, the correlation measures of PHIFLEs can be regarded as the closeness degree between two variables.Normally, there are two types of methods to resolve the correlation measures, the statistics-based point of view and the information energy-based point of view.A large number of samples are the basis for the former method; however, a relatively small set of PHIFLTSs are more common for us to deal with.Therefore, we carry out the correlation measure from the views of information energy.Motivated by the studies in FSs [44], IFSs [45], HFS [46], and HFLTS [19], we propose the concept of information energy in the PHIFLTSs, which can be defined as: be an LTS, ( ) {( , ( ), ( ) ) , the information energy of ( ) S E P can be defined as: With the use of the correlation and information energy, we propose the correlation coefficient between two PHIFLEs where , , 1,2,...

PHIFLTS s s P H I F L T S s P H I F L T S s C E P E P max E E P E E P
 .Since For convenience, the LSF ( ) , 0,1,...
is selected as the equivalent transformation functions f.The correlation coefficient between Finally, the correlation and correlation coefficient between ( ( ), ( )) 0.98

The MAGDM Method Based on the MULTIMOORA Approach under the PHIFLTSs Environment
Based on the novel operational laws and correlation measures, we propose a new MAGDM method for the probabilistic hesitant intuitionistic fuzzy linguistic term sets utilizing the MULTIMOORA approach in this section.

The Problem Statement
Here, let  .The multiple attribute group decision making problems considered in this paper can be represented as follows: where ( ) t ij E P represent the values associate with the criteria j c in alternative i a given by the decision maker t d , which is assessed as PHIFLEs.

Determine the Optimal Weights of Attributes Based on the TOPSIS and LP Optimization Method
The optimal weights of attributes  denote the similarity measures between the j-th elements in the i-th row of the collective decision matrix and the positive ideal decision H + and the negative ideal decision H − , respectively.It can be obtained as follows: where where

The MULTIMOORA Approach under the PHIFLTSs Environment
The MULTIMOORA approach is selected as the ranking techniques to solve our problem.First of all, we can calculate the scores of R based on Definition 5, which is depicted as follows: The vector normalization of ( ) S E can be derived by 1. Probabilistic Hesitant Intuitionistic Fuzzy Linguistic Ratio System (PHIFLRS) model The PHIFLRS model is defined as the arithmetic weighted aggregation operator, and we can obtain the first subordinate utility value 1 ( ) i U a of the alternative i a ( in descending order, and we obtain the first subordinate ranks of alternatives as Rank a a a  .

Probabilistic Hesitant Intuitionistic Fuzzy Linguistic Reference Point (PHIFLRP) model
For the reasons of avoiding the selected alternative that has poor performance under some specific attributes, we define the PHIFLRP model, which is determined by the worst performance of alternatives i a with respect to different attributes.
where ( )  ( ) ( ) / ( ) in descending order, and we obtain the third subordinate ranks of alternatives as Rank a a a  .
At last, the above three subordinate ranks of alternatives need to be aggregated comprehensively to derive the final ranking.Since different conditions of problems are resolved by the three models, they can be regarded as three attributes, which are PHIFLRS (C1), PHIFLRP (C2), and PHIFLFMF (C3).
Then, two matrices can be built based on the utility value , , In the traditional MULTIMOORA approach, the dominance theory was adopted to carry out the aggregation.However, as per the analysis in [23], this technique not only fails to consider the utility value of each alternative associated with each model but also has an extremely complex and timeconsuming operation process.Therefore, inspired by the method proposed in [23], we utilized the improving Borda rule to derive the final ranking.
Firstly, the vector of these three kinds of subordinate utility value can be normalized and derived the vector Subsequently, we can translate the ordinal values to the scores for reasons of aggregating the cardinal of an alternative, which can be depicted as the weight of each cardinal value.The final ranking can be derived by the values of i FS in descending order, which is defined as To clarify the proposed MAGDM method based on the MULTIMOORA approach under the PHIFLTSs environment, we summarize the specific processes and describe them in Algorithm 1.

Algorithm 1 The MAGDM method based on the MULTIMOORA approach
Step 1. Collect the individual linguistic evaluations values described as PHIFLEs from experts and construct the individual decision matrices ( ( )) . Go to the next step.
Step 2. The collective decision matrix can be built based on the frequency or the aggregation of probabilities associated with each linguistic term that appear in each individual decision matrix t D .Go to the next step.
Step 3. The score function and the variance value function of PHIFLEs defined in Definitions 5 and 6 are utilized to derive the PID H + and the NID H -associated with each attribute in the collective decision matrix R. Go to the next step.
Step 4. The PHIFLEs in R should be adjusted to the same probability set based on Definition 8 to calculate the correlation coefficients between ( ( )) Step 5.The LP optimization method should be utilized to determine the optimal weights of each attributes based on Equations ( 14)-( 16).Subsequently, we can derive the vector normalization of the scores of the collective decision matrix R based on Equations ( 17) and (18).
Step 6.The MULTIMOORA approach is utilized to calculate the final ranking for all alternatives.Firstly, the subordinated methods as PHIFLRS are used to obtain the utility values and the ranks of alternatives as 1 R ank based on Equation (19).
Step 7. The second subordinated methods as PHIFLPR in the MULTIMOORA approach are used to obtain the utility values 2 ( ) Step 10.The improving Borda rules are utilized to obtain the final ranking of alternatives.Firstly, the normalized vectors of utility value N U are obtained based on Equation (22).Secondly, the ordinal values need to be translated to the scores based on the Borda rules for the reasons of aggregating the cardinal of each alternative.Finally, the final ranking can be derived by the values of i FS ( 1,2,..., ) i m  in descending order based on Equation (23).End.

Application Example
For reasons of comparison, this part applies the same application example as [27] and [47] to validate the proposed theory and decision making models.Seven experts, denoted by  , , , , , }  d d d d d d d , need to select the most profitable method for their savings.Through The LTSs used on all attributes are united to .All experts assess the performance of the alternatives with respect to each attribute based on the LTS, and their opinions in terms of HIFLTSs are listed in Tables 1-3.In addition, we need to consider the semantics of unbalanced LTSs, and the second type of LSF shown as Equation ( 4) is selected for all attributes.5 and 6.  7 and 8.    .Firstly, the normalized vectors of utility value N U are obtained based on Equation (22).
Secondly, the ordinal values need to be translated to the scores based on the Borda rules for reasons of aggregating the cardinal of each alternative.Finally, the final ranking can be derived by the values of i FS ( 1,2,..., ) i m  in descending order based on Equation (23), which are listed in Table 9.

Comparison with Other Related Research Works
In this section, we present a comparison with the existing related works to show the advantages and innovation of the proposed method.

Comparative Analysis from the Numerical Points with the existing method
The authors in [27] firstly defined the PHIFLTSs and the corresponding application in MAGDM.The proposed method in this paper differs from the approach in [27] with respect to the following points: (1) The operations over PHIFLTSs proposed in [27] are directly based on multiplying the probabilities by the subscript of the corresponding linguistic terms, which may lose the associated probability information.For example, the PIS and NIS obtained in [27]  {s (0.29),s (0.5),s (0.21)},... .Obviously, the former one loses the probability information, which may lead to the failure of the final decision-making; (2) The method in [27] does not consider the situations of unbalanced linguistic terms on the PHIFLTSs environment; however, this paper selects the second type of linguistic scale function for all attributes to deal with the semantics of unbalanced LTSs; (3) Table 10 shows the difference of the ranking orders between the work in [27] and this paper, which is due to the effect of the innovation in the basic operational laws, the unbalance LTSs, the distance measure and the MULTIMOORA approach proposed in this paper.It can be clearly seen that the method proposed in this paper is an advantages and innovation to solve MAGDM problems over the PHIFLTSs environment.

Comparative Analysis from the Numerical Points with HIFLTSs
For purposes of comparison, the application example presented in Section 7 will be solved based on the TOPSIS method with HIFLTSs, that is to say, the probabilities associated with each linguistic term will be omitted.Then, we can obtain the final ranking as 5 a a a a a     .Obviously, the difference is due to the effect of probabilities associated with each linguistic term, which verified clearly the importance of the probabilities in the whole process of decision making.

Conclusions
In this paper, we proposed a new MAGDM model based on the MULTIMOORA approach under the PHIFLTSs environment, which has the following superiorities: (1) Utilizing the new score function and variance value function of PHIFLTSs makes the comparison between PHIFLEs more accurate and rational; (2) The unbalanced linguistic term sets were considered during the whole decision processes; (3) The correlation coefficient of the PHIFLTSs was used to measure the similarity instead of the current distance function in order to alleviate the drawbacks; (4) The decision results are robustness, simplicity, and effectiveness based on the three aggregation subordinated methods with the use of the extension of the MULTIMOORA approach.
For future research, we shall further analyze the PHIFLTSs in MAGDM problems and extend the proposed methodology to make it apply widely.Modeling the dynamic environment that allows the individual sets of attributes and alternatives is an interesting topic for further discussion.We also want to extend the PHIFLTSs to the probabilistic hesitant intuitionistic fuzzy linguistic preference relation (PHIFLPR), which can be utilized in more extensive decision making situations.

.
2,..., i m, =1,2,..., j n , LP optimization method should be formulated to determine the * ( 1,2,..., ) The objective function of the optimization model aims at minimizing the relative similarity from the H − and maximizing the relative similarity from the H + simultaneously.
are benefit attributes and cost attributes, respectively.The alternatives i a ( 1,2,..., ) i m  are ranked by the value of 1 ( )( 1,2,..., ) i U a i m  the second subordinate ranks of alternatives as

3 .
Probabilistic Hesitant Intuitionistic Fuzzy Linguistic Full Multiplicative Form (PHIFLFMF) model We can obtain the third subordinate utility value 3 ( ) i U a of the alternative i a ( 1,2,..., ) i m  based on the geometric weighted aggregation operator.


and the PID H + and the NID H -.

Step 8 .Step 9 .
The third subordinated methods as PHIFLFMF in the MULTIMOORA approach are used to obtain the utility values 3 ( ) The above three subordinate ranks of alternatives aggregated to obtain the final ranking.Two matrices can be built based on the utility values interviews and surveys, four attributes, which are c1: Risk factor, c2: Growth, c3: Quick refund, c4: Complicated documents, are used to evaluate five alternatives: x1: Real estate, x2: Stock market, x3: T-bills, x4: National saving scheme, and x5: Insurance company.
decision matrix R, which are listed in Tables

Definition 5. Let
) ( ) { (0.2), (0.2), (0.5), (0.1),{ (0.3), (0.7)}} PHIFLEs.For convenience, the first type of linguistic scale function ( ) , 0,1,...Based onDefinition 5, we can obtain 1 ( ( )) 0.37 S S E P  , 2 ( ( )) 0.447 S S E P  , 3 ( ( )) 0.37 S S E P  .Since s E P is superior to 2 ( ) S E P .However, for the reason of 1 2 To overcome such limitations, we propose a new distance measure based on the adjusted PHIFLEs, which is defined as follows: s E P than 2 ( ) s E P 、 , which goes against the fact.S E P can be defined as Utilizing the proposed method in Example 4, we can obtain the Euclidean distance as s are the subscripts of the membership and the non-membership linguistic terms set, respectively.The linguistic scale function f is determined based on the practical situations.0   , if 1   , the distance measure proposed in Definition 9, is a Hamming distance measurement; if 2   , Definition 9 is an Euclidean distance measurement.sEP is closer to 1 ( ) s E P than 2 ( ) s E P 、, and this is more rational.(( ), ( )) ( ( ), ( ))

Table 1 .
The valuation values provided by decision makers

Table 2 .
The evaluation values provided by decision makers

Table 3 .
The evaluation values provided by decision makers

Table 4 .
The PHIFLEs collective decision matrix R .for reasons of comparison between PHIFLEs, the score function and the variance value function of PHIFLE defined in Definitions 5 and 6 are utilized to obtain the positive ideal decision H + and the negative ideal decision H -associated with each attribute {s (0.14),s (0.36),s (0.36),s (0.14)}, {s (0.16),s (0.42),s (0.42)}

Table 5 .
The score of the collective decision matrix R. c1:

Table 6 .
The PID H + and the NID H -associated with each attribute j c j  In Step 3, to calculate the correlation coefficients between the collective decision matrix R and the PID H + and the NID H -for all the alternatives ( 1,2,...,5) i a i  i Sim  and i Sim  between the adjusted PHIFLEs of the collective decision matrix R and the PID H + and the NID H -for all alternatives ( 1,2,...,5) i a i  , which are listed in Tables

Table 7 .
The similarity i Sim  between the adjusted PHIFLEs of R and the PID H + . c1:

Table 8 .
The similarity i Sim  between the adjusted PHIFLEs of R and the NID H -. c1:

Risk Factor c2: Growth c3: Quick Refund c4: Complicated Documents Requirement
In Step 8, the second subordinated methods as PHIFLPR in the MULTIMOORA approach are used to obtain the utility values 2 ( ) (0.1934,0.2431,0.1697,0.2263,0.3039) In Step 11, the improving Borda rules are utilized to derive the final ranking of alternatives i

Table 10 .
Comparison of results.