A Parametric Method for Ranking Intuitionistic Fuzzy Numbers and Its Application to Solve Intuitionistic Fuzzy Network Data Envelopment Analysis Models

In this study, a new parametric method is proposed to rank intuitionistic fuzzy numbers in a general form. One of the advantages of the proposed method is that the decision maker’s idea is taken into account by selecting appropriate amounts of decision level and hesitation degree parameters. In some illustrative examples, the superiority of the proposed method over some other approaches is demonstrated. Furthermore, to show the ability of the method to solve intuitionistic fuzzy optimization problems, the proposed method is applied to solve intuitionistic fuzzy network data envelopment analysis (IFNDEA) problems. Also, in three appropriate examples, the validity of the suggested method and its capacity to solve real-world problems are illustrated.


Introduction
Intuitionistic fuzzy sets (IFSs), first proposed by Atanassov [1,2], are a generalization of Zadeh's fuzzy sets [3,4] to model noncrisp and uncertain sets. For these sets, we define both membership and nonmembership functions. In this case, we will be able to define a hesitation function that is the difference between the "membership function" and "one minus nonmembership function." us, by IFSs, we can efficiently model imperfect information. Many researchers provided enormous works on IFSs in both theory and applications (see e.g., [5][6][7][8]). As an application, Akram et al. [9] proposed a novel decision-making method based on hypergraphs in the intuitionistic fuzzy environment and applied it to real-life problems.
In 1994, by proposing intuitionistic fuzzy numbers (IFNs), Burillo et al. [10] generalized a fuzzy numbers concept and paved the way for applying intuitionistic fuzzy logic in real-life problems such as decision-making and risk analysis. For example, Liu and Wang [11] and Wang and Liu [12] proposed new Schweizer-Sklar Maclaurin and Einstein operation rules for IFNs and applied them to solve real-life decision-making problems.
Ranking IFNs is one of the important subjects in IFS theory, and a considerable number of articles have been published in this field. In these papers, different approaches are employed to compare and rank IFNs. For example, Grzegorzewski [13] generalized the ranking approach of fuzzy numbers in Grzegorzewski [14] and suggested a new method for ranking IFNs by defining concepts of expected interval and expected value. Later, Ye [15] calculated the expected value of a TraIFN that was introduced by Grzegorzewski [13] and applied it to a trapezoidal intuitionistic fuzzy multicriteria decision-making problem. Using four pairwise ranking functions, Mitchel [16] defined the average intuitionistic fuzzy rank index and proposed an approach to compare each pair of IFNs. Using score and accuracy functions of trapezoidal intuitionistic fuzzy numbers (TraIFNs), Jianqiang and Zhong [17] proposed a stepwise ranking method and applied it to a multicriteria decision-making problem. Wang and Zhang [18] ranked TraIFNs by converting them into interval numbers. Nehi [19] proposed a generalization of the characteristic value first suggested by Chiao [20] and ranked IFNs. In the work of Nan et al. [21], by introducing average indexes of membership and nonmembership functions, TraIFNs are ranked, and then, the approach is applied to matrix games with payoffs. Later, in a note, Verma and Kumar [22] modified some mathematical incorrect assumptions in the work of Nan et al. [21]. Value and ambiguity concepts were first introduced by Delgado et al. [23], and employing these concepts is a famous approach for ranking IFNs (see e.g., [24][25][26][27][28]). For example, in the work of Li [26], value and ambiguity are used to rank nonnormal triangular intuitionistic fuzzy numbers (TriIFNs), and then, the proposed method is applied to a personnel selection problem. In addition, Chutia and Saikia [25] generalized value and ambiguity concepts for membership and nonmembership functions and proposed a ranking method for IFNs using α and β cuts. ey also applied the proposed method in a risk analysis problem with intuitionistic fuzzy data. Wan [29], using possibility mean, possibility variance, standard deviation, and variance coefficient of TriIFNs, introduced possibility variance coefficients of membership and nonmembership functions to create a lexicographic ranking method. Das and Guha [30] ranked TraIFNs by calculating the centroid points and compared them. Li and Chen [31] proposed a distance index between two arbitrary TraIFNs based on α-cuts, which was, then, applied to the intuitionistic fuzzy multicriteria group decision-making problem. In the work of Lakshmana Gomathi Nayagam et al. [32], eight score functions are defined and applied to rank TraIFNs. By appropriate examples, the method is compared with some other approaches in this field. Also, Nayagam et al. [33] generalized a lexicographic ranking method to propose a total ordering of IFNs using α-cut and β-cut sets. In another work, Nayagam et al. [28] proposed a parametric score function using improved value and ambiguity indexes to rank TraIFNs and applied it to solve a multicriteria decision-making problem. Aggarwal and Gupta [34] defined an index to rank generalized symmetrical TraIFNs and solved an intuitionistic fuzzy solid transportation problem. Prakash et al. [35] introduced a ranking method by calculating the centroid of TraIFNs. Singh and Yadav [36] defined score and accuracy indices for the LR type of normal IFNs and proposed an order relation. Furthermore, by applying these indices, they proposed a method for solving a fully intuitionistic fuzzy linear-programming problem. Later, Canedo and Morales [37] analyzed this approach and modified it to obtain unique optimal solutions. Darehmiraki [38] extended the method by Shureshjani and Daremiraki [39] to rank IFNs and applied it to the partner selection problem. Table 1 summarizes the abovementioned studies based on the relevant themes. Although many methods have been proposed, there is still not a method to rank different kinds of IFNs, which can also be easily employed in optimization problems such as Data Envelopment Analysis (DEA) under intuitionistic fuzzy environment via the decision-maker idea.
DEA is a nonparametric and linear-programming-based approach for evaluating the efficiency of decision-making units (DMUs) with multiple inputs and outputs data. From the first work by Charnes et al. [48], there are many studies on both theory and applications in this field. In DEA models, a DMU is considered as a black box with its interior structures not taken into account. However, in many applications, DMUs have important interior structures, and ignoring them has led to inappropriate results. Network DEA is an attempt to consider the interior structures of DMUs. Among all the proposed structures, the simplest and most basic network DEA model is two-stage network DEA models. Also, in conventional DEA models, all the data should be crisp and certain which is not always possible in real applications. To solve these problems, fuzzy numbers and IFNs, fuzzy DEA, and intuitionistic fuzzy DEA (IFDEA) models are proposed, respectively (see Emrouznejad et al. [49] for a literature review on fuzzy DEA models and its applications until 2013). Fuzzy DEA can efficiently model DEA problems with noncrisp data, but when we face both noncrisp and uncertain data, IFDEA models are the best alternative. Although many papers have been published on fuzzy DEA, the proposed works in the IFDEA field are limited. For example, Xu et al. [40] introduced a fuzzy superefficient cross-DEA model to measure the efficiency of logistics enterprises using an interval-valued intuitionistic fuzzy Bayesian network. Using a weighted aggregation operator, Razavi Hajiagha et al. [41] proposed an approach to evaluate the efficiency of DMUs in an intuitionistic fuzzy BCC model. Puri and Yadav [42], by applying the expected interval and expected value (Grzegorzewski [13]), proposed an index to estimate triangular intuitionistic fuzzy input and output data in an IFDEA model. Similar to the work of Puri and Yadav [42], by applying the expected value (Grzegorzewski [13]), Singh [43] transformed an intuitionistic fuzzy DEA/AR model into a conventional DEA/AR model. Using α-cut and β-cut concepts in IFS theory, Arya and Yadav [44] introduced lower and upper bounds for the efficiency of DMUs in intuitionistic fuzzy SBM and superefficiency intuitionistic fuzzy SBM models. In another paper, Arya and Yadav [45] proposed appropriate models based on α-cuts and β-cuts to provide lower and upper bounds of efficiency measures in IFDEA, and then, by applying the index proposed by Chen and Klein [50], they ranked DMUs. Also, their method was applied to evaluate and rank 16 hospitals in India. Moreover, Arya and Yadav [46] developed intuitionistic fuzzy BCC and intuitionistic fuzzy superefficiency BCC models with triangular IF inputs and outputs data. ey developed the proposed approach to determine the efficiencies and rankings of DMUs in the presence of infeasibility. Finally, in Ameri et al. [47], by substituting the IF inputs and IF outputs of a parallel IFNDEA model with their assigned expected values (Grzegorzewski [13]), the parallel IFNDEA model transforms into a linear program. In a case study, their method was applied for the self-assessment of an Iranian hospital (See Table 1).
Two important limitations of the abovementioned approaches in IFDEA are that, almost all of them work for a specific group of IFNs (triangular IFNs) only and the decision maker does not play any role in the decision-making process. In addition, they are not easily generalizable to intuitionistic fuzzy network DEA (IFNDEA) models. Despite a large number of papers published in the IFDEA field, there has been one published paper (by Ameri et al. [47]) in the IFNDEA field. e main motive of this study is to overcome the shortcomings listed above by assigning an appropriate parametric index to IFNs considering two main concepts in IFSs, i.e., alpha cut and hesitation degree. In this case, by proper selection of the decision level (alpha cut) and the hesitation degree, the decision maker's idea is accounted for in the decision-making process, and less information is lost in IFNs. is parametric index can be applied to compare and rank IFNs.
Also, by substituting the IF data with their assigned parametric indexes, an intuitionistic fuzzy decision-making problem transforms into a parametric decision-making problem that can be easily solved by proper selection of the parameters. In this study, a new parametric method is developed to rank IFNs. Also, some reasonable properties to rank IFNs are examined based on the work of Wang and Kerre [51]. To show the ability of the method in real-life problems, it is applied to solve the IFDEA problem in singlestage, series, and parallel structures. As will be seen, the proposed method is also easily generalizable to other kinds of IFNDEA problems. e remainder of this research is organized as follows. In Section 2, some preliminaries of IFSs, IFNs, and the related arithmetic operators are introduced. In Section 3, using decision level and hesitation degree concepts, a new parametric method is proposed to rank IFNs in a general form, by appropriate examples, the advantages of the proposed method are illustrated, and the obtained results are finally compared with some other approaches in this field. In Section 4, we applied our proposed method to solve the IFNDEA problem. In this method, IFDEA and IFN-DEA models are transformed into parametric DEA and parametric NDEA models, respectively. ree examples illustrate the method. Finally, Section 5 concludes the paper.

Preliminaries
is section includes some basic concepts and notions. To easily understand the equations in this paper, some mathematical symbols are summarized in Table 2.
Definition 1 (see [1,2]). Let X be a fixed universe. An IFS A in X is given as Ranking IFNs Studies to develop a ranking method using the expected value concept Grzegorzewski [13], Ye [15], and Jianqiang and Zong [17] Studies to develop a ranking method using the statistical viewpoint Mitchell [16] and Wan [29] Studies to develop a ranking method using average indexes of membership and nonmembership functions Nan et al. [21] and Verma and Kumar [22] Studies to develop a ranking method using the distance index Wang and Zhang [18], Aggarwal and Gupta [34], and Li and Chen [31] Studies to develop a ranking method using score and accuracy indices Lakshmana Gomathi Nayagam et al. [32], Nayagam et al. [33], Singh and Yadav [36], and Canedo and Morales [37] Studies to develop a ranking method using value and ambiguity concepts Li [26], Li et al. [27], Chutia and Chutia [24], Nayagam et al. [28], and Chutia and Saikia [25] Studies to develop a ranking method using a centroid concept Das and Guha [30] and Prakash et al. [35] Studies to develop a ranking method using an integral value of membership and nonmembership functions Nehi [19] and Darehmiraki [38] IFDEA models A superefficient cross-DEA model based on the Bayesian network in the interval-intuitionistic fuzzy environment Xu et al. [40] An intuitionistic fuzzy BCC model Razavi Hajiagha et al. [41] Optimistic and pessimistic IFDEA models with triangular intuitionistic fuzzy data Puri and Yadav [42] An intuitionistic fuzzy DEA/AR model with triangular intuitionistic fuzzy data Singh [43] Intuitionistic fuzzy SBM and superefficiency intuitionistic fuzzy SBM models with triangular intuitionistic fuzzy data Arya and Yadav [44] An intuitionistic fuzzy CCR model with triangular intuitionistic fuzzy data Arya and Yadav [45] Intuitionistic fuzzy BCC and intuitionistic fuzzy superefficient BCC models with triangular intuitionistic fuzzy data Arya and Yadav [46] IFNDEA models Study to develop a parallel intuitionistic fuzzy network DEA model Ameri et al. [47] Complexity 3 where μ A (x), v A (x) are the degrees of membership and nonmembership functions from X to [0, 1], respectively.
Definition 2 (see [1,2]). e hesitation function of an IFS A is defined as It is clear that 0 ≤ π A (x) ≤ 1. is function indicates the hesitation degree of an element x (indeterminacy) belonging or not belonging to A.
In fuzzy sets, we do not have a hesitation degree (hesitation degree is zero). us, we can consider each fuzzy set as a special case of an IFS as Definition 3 (see [5]). An IFN A � (μ A , v A ) is an IFS in the set of real numbers R with the membership function as and the nonmembership function as where f A and r A are nondecreasing continuous functions from R to [0, 1] and g A and h A are nonincreasing con- In this definition, w is the maximum degree of the membership function and u is the minimum degree of the nonmembership function. So, according to the Definition 1, 0 ≤ w, u ≤ 1 and 0 ≤ w + u ≤ 1 ( Figure 1). In Definition 3, if w � 1 and u � 0, then, we have a normal IFN.
Definition 4 (see [2]). α-cut, β-cut, and (α, β)-cut sets of an IFN A � (μ A , v A ) are crisp subsets of R that are defined as follows: α-cut set: Li et al. [27] proposed a new approach for defining a TriIFN, and here, we generalized their approach for TraIFNs as follows:

Definition 5. A TraIFN
A is an IFN with membership and nonmembership functions as follows:

Mathematical symbols Meanings A
An intuitionistic fuzzy number Hesitation function of A w Maximum degree of the membership function u Minimum degree of the nonmembership function e assigned function to A for the selected decision level (α) and hesitation degree (k) m(r), n(r), and n′(r) Bounded monotonic decreasing left continuous functions m(r), n(r), and n′(r) Bounded monotonic increasing left continuous functions θ o Intuitionistic fuzzy efficiency score of the decision-making unit under evaluation x ij e i th intuitionistic fuzzy input data of the j th decision-making unit y rj e r th intuitionistic fuzzy output data of the j th decision-making unit z dj e d th intuitionistic fuzzy intermediate data of the j th decision-making unit v i e given weight to intuitionistic fuzzy input data i u r e given weight to intuitionistic fuzzy output data r w d e given weight to intuitionistic fuzzy intermediate data d 4 Complexity In Definition 5, if w � 1 and u � 0, then, we have a normal TraIFN. Also, TriIFN is a special case of TraIFN in which b � c and b ′ � c ′ . Figure 3 shows the hesitation degree of belonging or not belonging an element x to TraIFN A.
By Definition 4, α-cut, β-cut, and (α, β)-cut sets of a TraIFN A � 〈(a, b, c, d; w), (a ′ , b ′ , c ′ , d ′ ; u)〉 are the following closed intervals, respectively, We can see that the (α)-cut set of a TraIFN A can be calculated from the intersection of the closed intervals obtained from the α-cut and β-cut sets of A. ese three cut sets (α-cut, β-cut, and (α, β)-cut sets) remove some of the elements of a given TraIFN, whose degrees of its membership and nonmembership functions do not satisfy the levels given as arguments of the corresponding operators.
Definition 6 (see [25]). Let Complexity two arbitrary TraIFNs and λ ∈ R; then, the arithmetic operations are stipulated as follows: Ma et al. [52] proposed a new approach for defining a fuzzy number. In this approach, the variable on the vertical axis in Figure 1 is considered as an independent variable. By generalizing this approach, we obtain the following definition for IFNs.
, m(r)) of functions m(r) and m(r) ; 0 ≤ r ≤ w which satisfy the following requirements: and v A ′ is a pair (n(r), n(r)) of functions n(r) and n(r); u ≤ r ≤ 1 which satisfy the following requirements: (1) n(r) is a bounded monotonic decreasing left continuous function (2) n(r) is a bounded monotonic increasing left continuous function (3) n(r) ≤ n(r), u ≤ r ≤ 1 So, by Definition 7, we can represent a TraIFN as Figure 4).
As mentioned above, in Definition 7, the variable on the vertical axis is considered as an independent variable. It provides the ability to obtain the location of the left and right parts of the membership and nonmembership functions of an IFN on the horizontal axis for different amounts of the independent variable. In the following, we will use this property to assign an appropriate parametric index to IFNs.

Ranking Intuitionistic Fuzzy Numbers (IFNs)
In real-world applications, dealing with imprecise information is a common problem, and the uncertainty is unavoidably involved in every real-world problem. To model these problems, we should correctly incorporate the uncertainty concept into the problem description. Because of the ability of the intuitionistic fuzzy numbers to express imprecise information (by applying the hesitation concept), they are a useful tool to model these problems. One of the first important issues after the definition of intuitionistic fuzzy numbers is how to compare and rank them. e comparison and ranking of intuitionistic fuzzy numbers are complicated, and different methods have been proposed to address this problem.
Engaging the decision maker in the decision process provides flexibility and is a fascinating phenomenon [24]. Using α-cuts, Shureshjani and Darehmiraki [39] proposed a parametric method to rank fuzzy numbers based on the decision maker's opinion. In addition to α-cuts, we have another important concept in IFSs, namely, hesitation. Hesitation in an IFN, as shown in Figure 3, is an area between "one minus nonmembership function" and "membership function." It is clear that, optimally, this area will be added to the membership function, and in the most pessimistic view, we should consider this area as a part of the nonmembership function. us, in ranking IFNs, it would be appropriate to consider a degree of optimism (or pessimism) on the part of the decision maker in the decisionmaking process.
In this section, to consider the abovementioned points, we generalize Shureshjani and Darehmiraki's method [39] and propose a new parametric method for ranking IFNs based on the decision level and hesitation degree parameters.

e Proposed Index. For an arbitrary IFN
we assign the following index: where n ′ (r) and n ′ (r) belong to 1 − v A ′ (r). e selection of α and k parameters depend on the decision maker's idea. From the S α,k index definition, we can see that, after the selection of α parameter, only the elements of an IFN with the membership and hesitation values of larger than or equal to α will be important. erefore, if we choose α parameter close to one, we have a "high-level decision," and a "low-level decision" is made when the selected α is close to zero. 6 Complexity If we rewrite the abovementioned formula, we will have It is clear that if the decision maker sets k � 0, we have a pessimistic decision because we have not considered the hesitation area in our decision (hesitation area is added to the nonmembership function). Also, if the decision maker sets k � 1, an optimistic decision has been made because all the hesitation areas are added to the membership function. It should be noted that any choice of k (0 < k < 1) reflects the desirable hesitation degree that is considered by the decision maker in ranking IFNs process.
So, the assigned S α,k index to a TraIFN will be as follows: (13) Figure 6 shows the geometric meaning of the S α,k index for an arbitrary TraIFN A with k � 0 and k � 1, respectively.
We can see that the value of the S α,k index is the summation of the marked areas with solid and empty dots. Complexity e following theorems demonstrate the reasonable behavior of the proposed index for specific types of IFNs. Proof 〉 be two arbitrary TraIFNs where w 1 � w 2 and u 1 � u 2 ; then, we have Proof. From Definition 6, we have Figure 6: e geometric representation of S α,k index with k � 0 (a) and k � 1 (b). 8 Complexity Proof. From Definition 6, we have Let λ ≥ 0; then, we have Now, if λ < 0, then, we have

Our Proposed Ranking
Approach. Based on the proposed S α,k index, we can give the following definitions to rank IFNs.
be two arbitrary IFNs; then, for a decision level higher than α and a selected hesitation degree k (0 ≤ α < min w A , w B , k ∈ [0, 1]), we have Wang and Kerre [51] suggested seven reasonable properties to evaluate the rationality of the proposed indices for ranking fuzzy numbers. Here, we examine the validity of our proposed index for ranking IFNs based on Wang and Kerre's [51] approach.
Let I be the set of intuitionistic fuzzy quantities and A be a finite subset of I. Proof. From A≿B and B≿C, we have ∀α, k: S α,k (A) ≥ S α,k (B) and S α,k (B) ≥ S α,k (C).
ese imply that ∀α, k: S α,k (A) ≥ S α,k (C) that means A≿C.  I and (A, B, C  In Example 1, we will see that the proposed method can efficiently compare intersected TraIFNs, and the decision maker can apply his/her preferences by choosing appropriate amounts of decision level (α) and hesitation degree (k). Similar to our proposed method, Darehmiraki [38] tried to extend Shureshjani and Daremiraki's [39] method to rank IFNs. However, Example 2 shows that there are some situations where this method led to inappropriate results. As mentioned in the introduction section, various methods have been proposed for ranking IFNs. All these methods can be classified into two categories: parametric and nonparametric. In the parametric approach, using value and ambiguity concepts is the most popular (see e.g. [25][26][27][28]). Among these papers, two of the newest are that of Nayagam et al. [28] and Chutia and Saikia [25]. In Example 3, our parametric method is compared with that proposed by Nayagam et al. [28] and Chutia and Saikia [25] and two other nonparametric methods (by Ye [15] and Nayagam et al. [32]).
In an attempt to generalize Shureshjani and Darehmiraki's [39] method to rank IFNs, Darehmiraki [38] proposed an index using α-and β-cuts. Practically, however, it is evident that this index leads to inappropriate results. e obtained results from Darehmiraki's [38] method for different amounts of the parameters are presented in Table 4. As observed, the results are not correct.
But from our proposed method, we have It is clear that, for all decision levels of α and hesitation degrees of k, we have (See Figure 9). So from Definition 9, we obtain A ≺ B that is a reasonable result.
Complexity ese parametric functions are compared and plotted at different decision levels of α (0 ≤ α < 1) and hesitation degrees of k (0 ≤ k ≤ 1) in Figure 10.
As can be seen from Figure 10, for different amounts of α and k, we obtain 5 different ranking results as follows: (1) If 0 ≤ α < 0.2, k ∈ [0, 1], then, S α,k (A) < S α,k (C) < S α,k (B), and from Definition 8, we have and so A ∼ α,k B ≺ α,k C (5) If 0.5 < α < 1, k ∈ [0, 1]: S α,k (B) < S α,k (A) < S α,k (C), and so B ≺ α,k A ≺ α,k C    By applying Chutia and Saikia's [25] method, after plotting value function (V α,β ) of A, B, and C, it can be seen that, for most values of α and β, we obtain V α,β (A) < V α,β (B) (see Figure 11), so for these values, we have A ≺ B. For specific amounts of α and β such as α � 0.8, β � 0.05, we have V α,β (B) < V α,β (A), so B ≺ A (see Figure 11). Also, for all amounts of α and β, we obtain V α,β (A) < V α,β (C), so we have A ≺ C (see Figure 11). Moreover, for different amounts of α and β, we obtain different ranking results between B and C (see Figure 11). So, by appropriate amounts of α and β, we can obtain the following ranking results: Ambiguity function (A α,β ) of A, B, and C is plotted in Figure 12 too. We can see that the equality condition between these IFNs will not occur by Chutia and Saikia's [25] method. erefore, unlike our method, the ranking results of A ≺ C ∼ B or A ∼ B ≺ C are not obtained by Chutia and Saikia's [25] method. Besides, unlike Chutia and Saikia's [25] method that apply a two-level approach to compare and rank IFNs, we assign a parametric index to IFNs which is dependent on the decision maker's idea. erefore, we can easily apply our method for solving optimization problems under intuitionistic fuzzy environments.
In another paper, Lakshmana Gomathi Nayagam et al. us, based on their index, we have C ≺ B ≺ A which is not achieved by Ye [15], Chutia and Saikia [25], and our proposed method.
From the abovementioned examples, we can see that the proposed parametric method has appropriate discrimination power to rank IFNs based on the decision level (α) and hesitation degree (k).

Application in Intuitionistic Fuzzy Network DEA (IFNDEA) Problems
In this section, to propose the advantages of the suggested index in application problems, we apply it to solve IFDEA and IFNDEA problems. One drawback is that they do not work with noncrispy and uncertainty data in real-life problems. As we know, IFNs can efficiently describe fuzzy appearing in the data in the presence of hesitation. By considering IFNs in DEA problems, we should solve an IFDEA problem. Figure 14 shows an IFDEA model with m IF inputs (x ij , i � 1, . . . , m) and s IF outputs (y rj , r � 1, . . . , s).
ere have been suggested various DEA models for evaluating the efficiency of DMUs. Among them, the Charnes, Cooper, and Rhodes (CCR) model proposed by Charnes et al. [48] is the dominant DEA model. Considering IF inputs and IF outputs, the basic CCR model is converted to the following IFCCR model: � 1, . . . , m, r � 1, . . . , s, (29) where θ o is the IF efficiency score of the DMU under evaluation (DMU o ) and v i and u r are the given weights to IF input i and IF output r, respectively.
By substituting the IF inputs and IF outputs with their assigned S α,k functions, the abovementioned IFCCR model is transformed into the following parametric programming which is dependent on decision level (α) and hesitation degree (k). 14 Complexity where 0 ≤ α < w(w � min w x ij , w y rj /i � 1, . . . , m, r � 1, . . . , s, j � 1, . . . , n}) and 0 ≤ k ≤ 1. Clearly, after determining the amounts of decision level (α) and hesitation degree (k), the abovementioned parametric programming model will be converted to a simple linear program. us, we can easily compute the efficiency of DMUs for different amounts of α and k. Generally, the amounts of α and k will be selected based on the decision maker's idea. If the decision maker considers α close to one, only the elements of an IFN with high membership values are of importance, and so, the obtained efficiency results are made with a high-level decision. Also, by assuming α close to zero, we obtain efficiency results under a low-level decision [53]. Moreover, if the decision maker considers k close to zero, DMUs will be evaluated with a low hesitation degree (pessimistic decision), whereas for a high hesitation degree (optimistic decision), k should be chosen near to one.
Although we display our method for the IFCCR model, it can be easily generalizable to other kinds of IFDEA models.
For more explanation, consider the following example that is proposed by Arya and Yadav [45].

Example 4.
e following data are collected from 16 public hospitals of Meerut in the Indian state of Uttar Pradesh (Arya and Yadav [45]). ere are two IF inputs and two IF outputs. e IF inputs are the total number of beds (x 1 ) and the sum of pathologists and laboratory technicians (x 2 ), and the IF outputs are the number of pathology operations (y 1 ) and the sum of plasters and tubal ligations (y 2 ) (see Table 5).
Using model (30), the efficiency of 16 public hospitals for different amounts of decision level (α � 0, α � 0.5, and α � 0.95) and hesitation degree (k � 0, k � 0.5, and k � 1) are calculated and shown in Table 6. As can be seen, different efficiency and ranking results are obtained for different amounts of α and k. For example, for a fixed decision level α � 0, if we set k � 0, then HCM will be efficient (θ � 1) and ranked 1, but for k � 1, HCM will be inefficient (θ � 0.852528) and ranked 2. Also, for a fixed hesitation degree k � 0.5, if we set α � 0, then the obtained efficiency result of GFD will be 0.517736 and ranked 8, but for α � 0.95, the efficiency result of GFD is 0.743246 and ranked 3. Also, from Table 6, for all the selected amounts of decision level (α) and hesitation degree (k), the BD, BLD, and GBD hospitals are recognized as efficient hospitals and ranked one. Inspired by the α-level-based approach in fuzzy DEA (see Emrouznejad et al. [49]), Arya and Yadav [45] developed an approach to solve IFDEA models based on α-cuts and β-cuts. e method provided lower and upper bounds of DMUs' efficiencies for different amounts of α and β. en, by applying Chen and Klein's [50] index, they combined the obtained lower and upper bounds to rank DMUs. is approach assigns a wide range of efficiencies to a DMU for different amounts of α and β. Specifically, when they combine these obtained efficiency ranges, it can lead to incorrect ranking results. To illustrate, consider the obtained ranking results for 16 public hospitals based on Arya and Yadav's [45] method that is shown in the last column of Table 6.
Evidently, there are fundamental differences between our proposed method and Arya and Yadav's [45] method in ranking DMUs.
For example, if we compare inputs and outputs of GBD and MC hospitals, the outputs of GBD will be larger than the outputs of MC, and the inputs of GBD are smaller than the inputs of MC. us, the obtained efficiency and ranking results of GBD should be better than that of MC. Our method confirms this result for different amounts of α and Complexity k. Despite that, MC is considered better than GBD by Arya and Yadav's [45] method (MC ranked 2 and GBD ranked 10). Similarly, the outputs of GBD are larger than the outputs of GC, and the inputs of GBD are smaller than the inputs of GC. Hence, GBD should be more efficient than GC. Our method confirms this result for different amounts of α and k, but by Arya and Yadav's [45] method, GC is considered better than GBD (GC ranked 1 and GBD ranked 10).

Intuitionistic Fuzzy Network DEA (IFNDEA) Problems.
Another drawback of traditional DEA models is that the structure of the decision-making units (DMUs) is assumed to be a black box, which takes no account of internal relations, operations, and structures. However, in real-world applications, we face DEA problems with important internal structures that cannot be ignored. Network DEA (NDEA) is an attempt to consider important interior structures of DMUs to DEA models. Assuming intuitionistic fuzzy data in an NDEA model, we will have an IFNDEA model. Among the various proposed structures, series and parallel are the two famous structures of NDEA models. But, although there are many papers on FDEA, IFDEA, and FNDEA fields, we found only one paper on IFNDEA (see [47]). In the work of Ameri et al. [47], IFNs are considered as inputs and outputs of a parallel NDEA model. In the following sections, we apply our proposed method in both series and parallel structures. Also, in a real-world application of the parallel IFNDEA models, the obtained results from our proposed method are compared with those of Ameri et al. [47].

Series Structure.
Among all the proposed NDEA models, two-stage network processes have the simplest structure.
Consider a simple two-stage IFNDEA model in series form with m IF inputs (x ij , i � 1, . . . , m), d IF intermediates (z dj , d � 1, . . . , D), and s IF outputs (y rj , r � 1, . . . , s) data ( Figure 15). Different approaches have been proposed to evaluate the efficiency of DMUs in a two-stage NDEA model (see Cook and Zhu [54]), the multiplicative method initially proposed by Kao and Hwang [55] being well-known. In this method, overall efficiency is defined as a production of the stage's efficiencies.  [45]).

DMUs
Inputs Outputs x 1 x 2 y 1 y 2 Table  6: Efficiency results of 16 public hospitals for different amounts of decision level (α) and hesitation degree (k).

DMUs Efficiency and ranking results
Ranking results from Arya and Yadav [45] α � If we consider IFNs for inputs, intermediates, and outputs data, the multiplicative two-stage IFNDEA model will be as follows: where θ o is the IF efficiency score of the DMU under evaluation (DMU o ) and v i , w d , and u r are the given weights to IF input i, IF intermediate data d, and IF output r, respectively. Now, by substituting the IF input, IF intermediate, and IF output data with their assigned S α,k functions, the twostage IFNDEA model will be converted to the following parametric program: where 0 ≤ α < w (w � min{w x ij , w z dr j , w y rj /i � 1, . . . , m, d � 1, . . . , D, r � 1, . . . , s, j � 1, . . . , n}) and 0 ≤ k ≤ 1.
Based on the decision maker's opinion, by appropriate selection of the decision level (α) and hesitation degree (k), model (32) will be converted to a linear program. e following numerical example is proposed to illustrate the abovementioned method.
Example 5. Consider five hypothetical DMUs with two IF inputs, two IF intermediate products, and two IF outputs in a two-stage NDEA structure ( Table 7).
As can be seen from Table 7, all the input, intermediate, and output data are assumed to be asymmetric TraIFNs. Expectedly, different selection of the decision level (alpha) and hesitation degree (k) led to different ranking and efficiency results. Using model (32), the efficiency results of DMUs are calculated and shown in Table 8 for different amounts of decision level (α � 0, α � 0.5, and α � 0.75) and hesitation degree (k � 0, k � 0.5, and k � 1). We can observe that different amounts of α and k have led to different efficiency and ranking results. For example, the efficiency score of DMU1 for α � 0, k � 1 (a low-level decision with a high hesitation degree) is 0.377979 and ranked 2, but for α � 0.75, k � 0 (a high-level decision with a low hesitation degree), the score is 0.114912 and it ranked 5. Also, the efficiency score of DMU4 for α � 0.75, k � 0 (a high-level decision with a low hesitation degree) is 0.498449 and ranked 1, but for α � 0.75, k � 1(a high-level decision with a high hesitation degree), the score is 0.223572 and it ranked 4. In most cases, DMU2 is the worst DMU and DMU3 is the best DMU. Also, DMU4 and DMU5 experience most changes in their rankings (from the first to fourth).
Since crisp or fuzzy numbers are special cases of IFNs, our proposed method can solve NDEA models with a collection of crisp, fuzzy, and intuitionistic fuzzy numbers as data too.

Parallel
Structure. Unlike series structures, there is no flow of intermediate products in the parallel structures. In parallel structures, the input and output data are divided among processes (Figure 16).
Consider a parallel IFNDEA model with q independent processes ( Figure 16). In this model, each process p uses IF inputs x (p) i , i ∈ I (p) to produce IF outputs y (p) r , r ∈ O (p) . In this system, the sum of the IF inputs used by all the q processes ( q p�1 x (p) ij � x ij ) is equal to the system IF inputs, and the sum of the IF outputs used by all the q processes ( q p�1 y (p) rj � y rj ) is equal to the system IF outputs. Based on Kao's [56] method, to evaluate the system and process efficiencies, first, the following IFNDEA model should be solved: ..., m, r � 1, ..., s.
By our proposed method, we should first substitute the IF inputs and IF outputs with their assigned S α,k functions. So, the abovementioned IFNDEA model will transform into   9, 11, 16; 0.9), (3, 9, 11, 17; Table  8: Efficiency results of five hypothetical DMUs for different amounts of decision level (α) and hesitation degree (k).
Consequently, after the selection of the decision level (α) and hesitation degree (k), we can solve the obtained linear program.
Let (v * i , u * r ) be an optimal solution of the model (34). We can calculate the process and system efficiencies by the following formulas: To apply our proposed method in a real-world application, we consider the following example proposed by Ameri et al. [47]. In Ameri et al.'s [47] method, by substituting the IF inputs and IF outputs with their assigned expected values (Grzegorzewski [13]), the model (34) is transformed into a linear program. Unlike our proposed method, in the work of Ameri [47], we are allowed only to use a specific group of IFNs (normal TraIFNs), and the decision maker's opinion is not considered in the decision-making process.
Example 6. Imam Khomeini is one of the most famous hospitals in Iran. For the purpose of self-assessment of the Imam Khomeini hospital, Ameri et al. [47] collected the hospital data in eight different periods. So, we have eight DMUs that are made by one hospital data in these different eight periods of time. e internal structure of the hospital department is considered as a parallel structure with four processes in the following sections: internal, general surgery, pediatric, and CCU. Table 9 shows the hospital data, and we can see that each process has three inputs and four outputs. e input data are the physician x 1 , the number of nurses x 2 , and the number of beds x 3 . Also, the output data are the bed occupancy percentage (BOP) y 1 , satisfaction with physician performance y 2 , satisfaction with nurse performance y 3 , and satisfaction with environmental conditions y 4 . As evident in Table 9, x 1 , x 2 , x 3 , and y 1 data are real numbers, while to control the uncertainty and for better modeling of the real data represented by organizations' decision makers, y 2 , y 3 , and y 4 data are represented by normal TriIFNs.
Using model 34, the efficiency results of the hospital and its sections are calculated for eight different periods of time (see Table 10). Unlike the previous example, in this example, there are no significant differences among the obtained efficiency results for the sections based on different values of the decision level (α) or hesitation degree (k). To explain, all the input data and an input datum are real numbers. Also, the other three output data are not represented by asymmetric or intersected TriIFNs. From Table 10, the greatest differences among the obtained efficiency results for different amounts of α and k is for the CCU section in part 6 (from 0.84959443 (α � 0, k � 1) to 0.91640810 (α � 0.95, k � 0)). We can see that there are not any efficient sections in parts 1, 2, 3, and 5. Also, the pediatrics section in parts 4 and 6, CCU section in parts 4 and 7, and general surgery section in part 8 are reported as efficient sections. e obtained efficiency results from Ameri et al.'s [47] method is shown in the last column of Table 10. As can be seen, although the problem is solved with two different approaches, the obtained efficiency results from our proposed method are nearly identical to those of Ameri et al. [47]. However, the general surgery section in part 6 from our proposed method is an inefficient section for all the selected amounts of α and k, but from Ameri et al.'s [47] method, it is an efficient section.

Process q
Process p Figure 16: A parallel system.

Complexity 21
As mentioned earlier, in our proposed method, the decision maker's idea is accounted for in the decisionmaking process by proper selection of the decision level (alpha cut) and the hesitation degree, and so, less information is lost in IFNs. But, in Ameri et al. [47], only a real number (expected value) is assigned to each IF data, and the decision maker does not play any role in the decisionmaking process.
Although we adopted our method in the multiplicative two-stage IFNDEA and parallel IFNDEA models, it will be easily generalizable to other kinds of IFNDEA models.
We can summarize the algorithm of the proposed method to solve IFNDEA models in real-world applications as follows: Step 1. An appropriate IFNDEA model is chosen based on the structure of the real-world problem Step 2. e input, intermediate, and output data of the problem are determined and represented by crisp, fuzzy, or intuitionistic fuzzy numbers Step 3. e appropriate decision level (α) and hesitation degree (k) are selected based on the decision maker's opinion Step 4. Input, intermediate, and output data are replaced with their assigned S α,k index Step 5. e transformed IFNDEA model is solved to obtain the efficiency score of the DMU under evaluation based on the selected decision level (α) and hesitation degree (k) As evident in the algorithm, after replacing input, intermediate, and output data with their assigned S α,k indexes, we should solve a transformed IFNDEA model that is a linear-programming model. Different methods are proposed for solving linear-programming problems such as the simplex method or Karmarkar method. us, the time complexity of our proposed method depends on the selected method for solving the obtained linear program.

Conclusions
IFNs can efficiently model noncrispy and uncertainty in real-world data. Ranking is the second most important issue after IFNs definition. In this study, a new parametric method was proposed for ranking IFNs, and by appropriate examples, its superiority over some other methods was explained. Engaging the decision maker in the decisionmaking process is an important advantage of the proposed method. It was observed that the decision maker's opinion about the decision level and hesitation degree parameters can affect the obtained ranking results. Many proposed methods to date rank only a special group of IFNs (e.g., TriIFNs, TraIFNs, or normal IFNs), but the proposed method is able to simultaneously rank all kinds of IFNs. Also, to show the ability of the proposed method in realworld applications, we applied our method to solve IFDEA, two-stage IFNDEA, and parallel IFNDEA models. e proposed approach is easily generalizable to solve other kinds of IFNDEA models too.
As future studies, we aim to apply our proposed method to multicriteria decision-making (MCDM) techniques, such as the analytic hierarchy process (AHP), best-worst method (BWM), and VIKOR, to solve intuitionistic fuzzy MCDM problems. Data Availability e data used to support the findings of this study are available from the corresponding author upon request.