# Ranking All DEA-Efficient DMUs Based on Cross Efficiency and Analytic Hierarchy Process Methods

Authors: Dariush Akbarian
Journal: Journal of Optimization
 (2015)

---

## Abstract

The aim of this paper is to present an original approach for ranking of DEA-efficient
DMUs based on the cross efficiency and analytic hierarchy process (AHP) methods. The
approach includes two basic stages. In the first stage using DEA models the cross efficiency
value of each DEA-efficiency DMU is specified. In the second stage, the pairwise comparison
matrix generated in the first stage is utilized to rank scale of the units via the one-step
process of AHP. The advantage of this proposed method is its capability of ranking extreme
and nonextreme DEA-efficient DMUs. The numerical examples are presented in this paper
and we compare our approach with some other approaches.


---

## Body

## 1. Introduction

Data envelopment analysis (DEA) is a tool for evaluation and measuring of the efficiency of a set of decision making units (DMUs) that consume multiple inputs and produce multiple outputs, first introduced by Charnes, Cooper, and Rhodes (CCR) [1] and extended by Banker, Charnes, and Cooper (BCC) [2]. DEA is used to establish a best practice group from among a set of observed units and to identify the units that are inefficient when compared to the best practice group [3–6]. One important issue in DEA which has been studied by many DEA researchers is to discriminate between efficient DMUs. Several authors have proposed methods for ranking the best performers ([7–14] among others). Perhaps the most widely known and applied ranking method is the super-efficiency DEA model. Developed by Andersen and Petersen [13], the main idea of this approach is to evaluate DMU after this performer itself is excluded from the reference set. The problem with super-efficiency DEA is that under certain conditions infeasibility occurs which limits the applicability of the technique (for details see [15]). Wu [4] utilized the concept of fuzzy for ranking of DMUs in the traditional DEA models. For this purpose, firstly the DMUs are evaluated by the CCR and the cross efficiency models. Secondly, a fuzzy preference relation is established. Finally, a priority vector of the preference is constructed and is used for ranking DMUs. Jahanshahloo et al. [8] proposed a ranking system based on changing reference set. In the proposed ranking system, the evaluation for efficient DMUs is dependent on the efficiency changes of all inefficient units due to its absence in the reference set while the estimate for inefficient DMUs depends on the influence of the exclusion of each efficient unit from the reference set. One popular method to rank multiple alternatives with respect to multiple criteria is analytical hierarchy process (AHP) method developed by Saaty [16]. The analytic hierarchy process is a multicriteria decision making technique and was used for complete ranking of DMUs by Sinuany-Stern et al. [17] and Jablonsky [18, 19] and references cited therein. Sinuany-Stern et al. [17] presented an AHP/DEA methodology for ranking DMUs, which integrates DEA and the AHP. In this AHP/DEA methodology, DEA was first run for each pair of DMUs separately to create a pairwise comparison matrix, which was then solved using the eigenvector method to produce a full ranking for all DMUs, efficient and inefficient. Their method has some problems, such as illogic comparing of two DMUs in a DEA model (see [20] for details). Alirezaee and Sani [20] extended Sinuany-Stern’s approach and proposed a new two-stage AHP/DEA methodology for ranking DMUs that removes these problems. They used together with the efficiency of units the influence rate of each unit to the other units in order to determine the elements of pairwise comparison matrix. Jablonsky [18] used an AHP model with interval pairwise comparisons for the evaluation and classification of efficient units and compares the results with super-efficiency DEA scores. In [19], Jablonsky presented an original procedure for ranking of DMUs in DEA models based on combination of AHP and DEA principles. He/she considered the DMUs as alternatives and the ratios of output/input as criteria. An extensive review of ranking models in DEA is given in Adler et al. [12]. The aim of this paper is to propose a method for ranking extreme and nonextreme efficient DMUs using cross efficiency and AHP methods. Our approach includes two stages. In the first step using DEA models the cross efficiency value of each DEA-efficiency DMU is specified. In the second stage, the pairwise comparison matrix generated in the first step is utilized to rank scale of the units via the AHP. This paper is organized as follows. Section 2 presents some basic DEA models and reviews some ranking method. In Section 3 we briefly review the AHP. Section 4 gives a brief introduction to the cross efficiency evaluation. The new ranking method for DEA-efficient DMUs is developed in Section 5. Numerical examples are examined in Section 6. The paper is concluded in Section 7.


## 2. Data Envelopment Analysis Models

Consider a set ofn DMUs which is associated with m inputs and s outputs. Particularly, DMUj (j∈J=1,…,n) consumes amount xij of input i and produces amount yrj of output r. Let Xj=(x1j,…,xmj) in which Xj≥0 and Xj≠0 and Yj=(y1j,…,yrj) in which Yj≥0 and Yj≠0. The production possibility set (PPS) of CCR model is defined as follows:
(1)TCCR={(X,Y)∣X≧∑j=1n‍λjXj,Y≦∑j=1n‍λjYj,λj⩾0,j∈J}
and similarly the production possibility set of BCC model is defined as follows:
(2)TBCC=(X,Y)∣X≧∑j=1n‍λjXj,Y≦∑j=1n‍λjYj,mmm∑j=1n‍λj=1,λj⩾0,j∈J.The input-oriented BCC and input-oriented CCR models, corresponding toDMUp, p∈J, are given by (3) and (4), respectively:
(3)min⁡ccθ-ϵ(∑i=1m‍si-+∑r=1s‍sr+)s.t.c∑j∈J‍λjyrj-sr+=yrp,r=1,…,smmcc∑j∈J‍λjxij+si-=θxip,i=1,…,mmmcc∑j∈J‍λj=1mmccλj≥0,j∈Jmmccsi-≥0,i=1,…,mmmccsr+≥0,r=1,…,smmccθfree(4)min⁡θ-ϵ(∑i=1mti-+∑r=1str+)s.t.c∑j∈Jλjxij+ti-=θxip,i=1,…,mmmcc∑j∈Jλjyrj-tr+=yrp,r=1,…,smmccλj≥0,j∈Jmmccti-≥0,i=1,…,mmmcctr+≥0,r=1,…,smmccθfree,
where ϵ is non-Archimedean small and positive number and si-, sr+, ti-, and tr+, i=1,…,m, r=1,…,s, are called slack variables belonging to R≥0. Note that si- and ti- represent input excesses; also sr+ and tr+ represent output shortfalls. Also θ and λj, j∈J, are real numbers and λj∈R≥0. Models (3) and (4) are called envelopment forms (with non-Archimedean number).DMU
p is said to bestrong efficient (CCR-efficient) if and only if θ*=1 and t+*=0 and t-*=0. Where the superscript (*) indicates optimality. In similar manner theBCC-efficient DMUs can be defined.The dual of model (4) (without ϵ, i.e., ϵ=0), which is called multiplier form (corresponding to DMUk), is as follows:
(5)max⁡θkk=∑r=1surkyrks.t.cc∑r=1surkyrj-∑i=1mvikxij≤0,j=1,…,nmmccc∑i=1mvikxik=1,mmcccurk≥0,r=1,…,smmcccvik≥0,i=1,…,m.The AP model is as follows [13]:
(6)min⁡θs.t.c∑j∈J-{p}λjyrj≥yrp,r=1,…,smmcc∑j∈J-{p}λjxij≤θxip,i=1,…,mmmccλj≥0,j∈J-{p}mmccθfree.
Jahanshahloo’s method corresponding inefficient DMUa is as follows [8]:
(7)min⁡θ-ϵ(∑i=1msi-+∑r=1ssr+)s.t.c-∑j∈J-{b}λjxij+θxia-si-=0,i=1,…,mmmcc∑j∈J-{b}λjyrj-sr+=yra,r=1,…,smmccλj≥0,j∈J-{b}mmccsi-≥0,i=1,…,mmmccsr+≥0,r=1,…,smmccθ,free.
The efficiency of strong efficiency DMUb will be denoted by Ω and will be given by
(8)Ωb=∑a∈Jn∂a,bn~,
in which Jn is the set of nonstrong efficiency DMUs and n~ is the number of nonstrong efficiency DMUs.Jahanshahloo et al. [21] used l1-norm in order to rank the extremely efficient DMUs in DEA models with constant and variable returns to scale, and the proposed method can remove the difficulties arising from AP and MAJ [11] models. Their proposed model is as follows:
(9)min⁡Γpc(X,Y)=∑i=1mxi-xip+∑r=1syr-yrps.t.∑j∈J-{p}λjxij≤xi,i=1,…,mmmc∑j∈J-{p}λjyrj≥yr,r=1,…,smmcxi≥0,i=1,…,mmmcyr≥0,r=1,…,smmcλjp≥0,j∈J-{p}.
As mentioned earlier, Sinuany-Stern’s approach [17] is a two-stage model for fully ranking model. In the first stage Sinuany-Stern run the following CCR model for each pair of DMUs A and B:
(10)max⁡EAA=∑r=1suryrAs.t.m∑r=1surkyrB-∑i=1mvixiB≤0,j=1,…,nmmmc∑i=1mvixiA=1,mmmc∑r=1suryrA≤1,mmmcur≥0,r=1,…,smmmcvi≥0,i=1,…,m.EAA is the efficiency score of unit A. The cross efficiency of DMU A is as follows:
(11)EBA=∑r=1sur*yrB∑i=1mvi*xiB
in which ur* and vi* are the optimal weights of model (10). In the similar manner EAB and EAA can be computed. The pairwise comparison matrix A=[aij] (i,j=1,…,n) is constructed as follows:
(12)aij=Eii+EijEjj+Eji,aii=1.
In the second stage using a one-step process (Step 3 in Section 3), the vector weights of the pairwise comparison matrix A are computed. Based on these weights, the efficient and inefficient DMUs are ranked.Alirezaee and Sani [20] proposed a two-stage AHP/DEA methodology for ranking DMUs. In the first stage they create the pairwise comparison matrix as follows.First, using CCR model (4) the relative efficiency of each DMU is determined. Then, each DMU is evaluated by the following model:
(13)min⁡θs.t.∑j=1,j≠knλjyrj≥yrp,r=1,…,smmc∑j=1,j≠knλjxij≤θxip,i=1,…,mmmcλj≥0,j=1,…,n,j≠kmmcθfree.
This model measures efficiency DMUp after excluding DMUk from PPS. The pairwise comparison matrix A is constructed as follows:
(14)A=[aij],aij=TEi+{TEij-TEj}TEj+{TEji-TEi},i,j=1,…,n,
in which TEi and TEji are the objective function of (4) and (13) at optimality, respectively.In the second stage, similar to AHP/DEA method, using a one-step process of AHP, the DMUs are ranked. In this paper we rank DMUs in the CCR model; in a similar way one can also rank DMUs in the BCC model.


## 3. Analytical Hierarchy Process (AHP)

AHP is a powerful tool for analysis of complex decision problems. AHP is employed for ranking a set of alternatives or for the selection of the best in a set of alternatives with respect to multiple criteria. In this section a brief discussion of AHP is provided. For more details see Saaty [16].Application of AHP to a decision problem involves four steps [23, 24].Step 1 (structuring of the decision problem into a hierarchical model).
The first step is to decompose a decision problem into elements. In its simplest form (Figure1), this structure comprises a goal at the topmost level, criteria (and subcriteria) at the intermediate levels, while the lowest level contains a set of the alternatives. A simple AHP model (Figure 1) has three levels (goal, criteria, and alternatives).Figure 1
Three-level hierarchy.Step 2 (making pairwise comparisons and obtaining the judgment matrix).
In this step, the elements of a particular level are compared with respect to a specific element in the immediate upper level. The more important (or more attractive) element receives higher rating (between 1 and 9) than another one that has the less important element [16]. The comparisons are used to form a matrix of pairwise comparisons called the judgment (pairwise comparison) matrix A. Each entry aij of the judgment matrix satisfied the three rules: aij>0; aij=1/aji; and aii=1 for all i, j. In fact, the number aij shows the relative importance of alternative Ai in comparison with alternative Ai with respect to criterion Cj or the relative importance of criteria Ci in comparison with criterion Cj with respect to goal. If aij=aik·akj, for all the entries of the matrix, then the matrix is said to be consistent. If the property does not hold for all the entries, the level of inconsistency can be captured by a measure called consistency ratio [16] (see next step).Step 3 (local weights and consistency of comparisons).
In this step, local weights of the elements are calculated from the judgment matrices using the eigenvector method (EVM). The normalized eigenvector corresponding to the principal eigenvalue of the judgment matrix provides the weights of the corresponding elements. Though EVM is followed widely in traditional AHP computations, other methods are also suggested for calculating weights, including the logarithmic least-square technique (LLST) [25, 26], goal programming [27], and recently DEAHP method [15, 16, 20, 22, 24–33]. When EVM is used, consistency ratio (CR) can be computed. The consistency index (CI) is calculated as follows:
(15)CI=λmax⁡-1n-1
in which n is the number of criteria and λmax⁡ is the maximum eigenvalue of the judgment matrix. Then, the consistency ratio (CR) is calculated as the ratio of consistency index and random index (RI); that is,
(16)CR=CIRI.
Random index is the consistency index of a randomly generated reciprocal matrix from the scales 1 to 9. Table 1 shows the value of RI sorted by the order of matrix. For a consistent matrix CR=0 and λmax⁡=n. A value of CR less than 0.1 is considered acceptable because human judgments need not be always consistent, and there may be inconsistencies introduced because of the nature of scale used. In the inconsistency case, λmax⁡≥n. If CR for a matrix is more than 0.1, judgments should be elicited once again from the DM till he∖she gives more consistent judgments.Table 1
Random index.
Matrix order
1
2
3
4
5
6
7
8
9
RI
0.00
0.00
0.58
0.9
1.12
1.24
1.32
1.41
1.45Remark 1.
If the pairwise comparisons matrix is as follows:(17)A=w1w1w1w2⋯w1wnw2w1w2w2⋯w2wn⋮⋮⋮⋮wnw1wnw2⋯wnwn,
then any weight derivation method, when applied to the above matrix A, should generate the true weights w1,w2,…,wn.Step 4 (aggregation of weights across various levels to obtain the final weights of alternatives).
In order to obtain final weights of the decision alternatives (elements at the lowest level), the local weights of elements of different levels are aggregated. For example, the final weight of alternativeAp is computed using the following hierarchical (arithmetic) aggregation rule in traditional AHP:
(18)FinalweightofAp=∑j‍(LocalweightofApwithrespecttoCj)mmmmmc×(LocalweightofcriterionCj).
By definition, the weights of alternatives and importance of criteria are normalized. The final weights represent the rating of the alternatives in achieving the goal of the problem.


## 4. Cross Efficiency Evaluation

The cross efficiency evaluation was first proposed by Sexton et al. [30] and was later investigated by Doyle and Green [34, 35]. It is often calculated as a two-stage process. In the first stage, the optimal weights of inputs and outputs are calculated for each DMU using the classical DEA formulation, for example, CCR model (5). Let urk*(r=1,2,…,s) and vik*(i=1,…,m) be the optimal solution to the LP (5). Then, θkk*=∑r=1surk*yrk is referred to as the optimistic efficiency or CCR efficiency of DMUk, which reflects the self-evaluation of DMUk, whereas θjk is referred to as a cross efficiency value of DMUj, which reflects the peer-evaluation of DMUk to DMUj(j=1,…,n; j≠k). Consider
(19)θjk=∑r=1surk*yrj∑i=1mvik*xij,k=1,…,n,j=1,…,n.
Therefore, each DMU will have multiple different efficiency scores, whose average reflects the overall performance of the DMU. Based on average cross efficiencies, DMUs can be compared and ranked.In general, the optimal weights obtained using classical DEA in the first stage are multiple solutions. Therefore, the valuesθkj* will change depending on these values in the second stage. In the second stage, these drawbacks are reduced and a suitable set of weights preserving the efficiency values obtained by DEA is selected for each DMU. To reduce this undesirable case, there are some model suggestions for preserving the self-efficiency scores, θkk*, obtained for each DMU. Sexton et al. [30] suggested introducing a secondary goal to optimize the input and output weights while keeping the CCR efficiency unchanged. According to Sexton et al. [30], the secondary goals could be, for example, either aggressive or benevolent. In the aggressive context, DMUk chooses the optimal solutions such that it maximizes self-efficiency and at a secondary level minimizes the other DMUs cross efficiency levels. The benevolent secondary objective would be to maximize all DMUs cross efficiency rankings. The most widely used secondary goals were suggested by Doyle and Green [34]. Wang and Chin [31] propose several new data envelopment analysis (DEA) models for cross efficiency evaluation by introducing a virtual ideal DMU (IDMU) and a virtual anti-ideal DMU (ADMU). These new DEA models determine the input and output weights by minimizing the distance of a DMU from IDMU or maximizing the distance between the DMU and ADMU or both of them simultaneously, without the use of the input and output information of the other DMUs as a secondary goal. Örkcü and Bal [33] presented goal programming models that could be used in the second stage of the cross evaluation. See also [31] for new alternative models for DEA cross efficiency evaluation. In this paper we present another secondary goal to the choice of weights among the alternative optimal solutions. For this aim, we first find the alternative optimal solutions corresponding to all extreme DEA-efficient DMUs using the algorithm proposed by Jahanshahloo et al. [36]. Then, the any nonnegative convex linear combination of these alternative optimal solutions can be considered as the chosen weights corresponding to DMUk (in this paper we consider the arithmetic mean of the alternative optimal solutions). The advantage of our approach is that in second stage it yields the positive weights (see Remark 2). We emphasize that one may choose another secondary goal to determine the input and output weights.Remark 2.
Any convex linear combination of optimal solutions of the following LP is optimal too:(20)max⁡cxs.t.mAx≤bmmcmcx≥0
in which A is the m×n matrix and x, c, and b are the n×1 and m×1 vectors, respectively.


## 5. Proposed Method

In this section we combine the cross efficiency method with AHP technique and provide a two-stage method for ranking all DEA-efficient DMUs using cross efficient approach and AHP technique. In the first stage, we obtain the CCR-efficient DMUs and the optimal solutions of the CCR model (5) corresponding to them. Note that model (5) may have the multiple solutions. In this case, we apply the proposed algorithm in [36] as secondary goal. Suppose (u1k*,u2k*,…,usk*,v1k*,v2k*,…,vmk*) are the optimal solutions (weights) of LP (5) corresponding to DMUk. The efficiency score of CCR-efficient DMUj using the weights of DMUk is as follows:
(21)θjk=∑r=1surk*yrj∑i=1mvik*xij,k=1,…,n,j=1,…,n.
Also, θjk is referred to as a cross efficiency value of DMUj, which reflects the peer-evaluation of DMUk to DMUj (j=1,…,n; j≠k). Furthermore, if the input-output weights of the LP (5) are not unique, the arithmetic mean of the alternative optimal solutions corresponding to DMUk is considered as the optimal weights of DMUk in order to obtain positive weights (see Remark 2). In the second stage, using a one-step process of AHP (Step 3 in Section 3), the vector weights of the pairwise comparison matrix, generated in the first stage, are calculated. The DEA-efficient DMUs are ranked based on these weights. The pairwise comparison matrix A is generated by θij as follows:
(22)A=aijinwhich,aij=θii+θijθjj+θji,i,j=1,…,n,θii=θii*.
Moreover, since the pairwise comparison matrices are constructed from objective data in our case, we are not concerned about the consistency of A.Remark 3.
It is evident thataii=1, aij=1/aji, i,j=1,…,n.It should be noted that since we intend to rank CCR-efficient DMUs, it is enough to compute the values ofaij for each pair of CCR-efficient DMUs DMUi and DMUj. After determining the pairwise comparison matrix A, AHP technique is applied for ranking CCR-efficient DMUs. In this paper we use the eigenvector method (EVM) for determining eigenvalue and eigenvector of the pairwise comparison matrix A.


## 6. Numerical Examples

In this section, two numerical examples are examined to illustrate the proposed method. Comparisons with other existing procedures will also be made. In Example2 we compare the proposed method with two other ranking methods which are based on AHP.Example 1.
We evaluated with our method the data of 20 branch banks of Iran. This data was previously analyzed by Amirteimoori and Kordrostami [22] and is listed in Table 2. Running the DEA model (4) will result in seven efficient units as 1, 4, 7, 12, 15, 17, and 20. Using the above-mentioned method the optimal weights corresponding to CCR-efficient DMUj, j=1,4,7,12,15,17,20, are chosen. The results are summarized in Table 3. The cross efficiency and pairwise comparison matrices of seven CCR-efficient DMUs are shown in Tables 4 and 5, respectively. For this pairwise comparison matrix, λmax⁡=7.0111825, CI=0.00186375, and CR=0.00141<1. So, pairwise comparison matrix is consistent enough. The corresponding eigenvector (weights) is shown in the last column of Table 5. Table 6 shows a comparison of our proposal and some other ranking approaches. All these approaches are implemented in input-oriented version under the condition of CRS. As reported in Table 6, DMU4 is the most efficient one in our method. According to the results, the rankings of DMUs by the four methods are almost similar; in particular, the results of our method are more similar to the method [21].Table 2
Example1: DMUs’ data (extracted from [[22], p. 689]).
Branch
Input
Output
Staff
Computer terminals
Space m2
Deposits
Loans
Charge
1
0.9503
0.70
0.1550
0.1900
0.5214
0.2926
2
0.7962
0.60
1.0000
0.2266
0.6274
0.4624
3
0.7982
0.75
0.5125
0.2283
0.9703
0.2606
4
0.8651
0.55
0.2100
0.1927
0.6324
1.0000
5
0.8151
0.85
0.2675
0.2333
0.7221
0.2463
6
0.8416
0.65
0.5000
0.2069
0.6025
0.5689
7
0.7189
0.60
0.3500
0.1824
0.9000
0.7158
8
0.7853
0.75
0.1200
0.1250
0.2340
0.2977
9
0.4756
0.60
0.1350
0.0801
0.3643
0.2439
10
0.6782
0.55
0.5100
0.0818
0.1835
0.0486
11
0.7112
1.00
0.3050
0.2117
0.3179
0.4031
12
0.8113
0.65
0.2550
0.1227
0.9225
0.6279
13
0.6586
0.85
0.3400
0.1755
0.6452
0.2605
14
0.9763
0.80
0.5400
0.1443
0.5143
0.2433
15
0.6845
0.95
0.4500
1.0000
0.2617
0.0982
16
0.6127
0.90
0.5250
0.1151
0.4021
0.4641
17
1.0000
0.60
0.2050
0.0900
1.0000
0.1614
18
0.6337
0.65
0.2350
0.0591
0.3492
0.0678
19
0.3715
0.70
0.2375
0.0385
0.1898
0.1112
20
0.5827
0.55
0.5000
0.1101
0.6145
0.7643Table 3
Example1: the optimal weights obtained by proposed method.
DMU
u
1
u
2
u
3
v
1
v
2
v
3
1
2743104820616200
7
1145504368295450
7
203344811696175
7
86942634280800
7
86942634280800
7
6674020361903369
7
4
10557780612150395395
2
7137450867780576185
6
8614624057487049745
3
15106565712832042015
6
15589621335236010925
6
99578464810982125259
20
7
192476232067007701700
53
3007266838159896400
31962169004909158050
53
129835340277577355800
53
159334283058868224250
53
35696050063119836479
53
12
4090713978428108425
17
4510281972189568425
17
228787181284072450
17
3200631113384367575
17
1036164623791167125
17
12048173865672775077
34
15
4221115845669960050
47
1647910063938525750
47
1461021310131305000
47
2289837544488904650
47
2192290568305410050
47
2546169695477270243
47
17
43696890644910227500
27
32066814438604442500
27
3549650609396840000
27
217430715355767500
2798202153611424125
3
76056230540954147200
27
20
65310366646698125
10656566179635625
70175639377223750
104635477353783125
5446944927631181.25
6814920193906181.25Table 4
Example1: cross efficiency matrix for CCR-efficient DMUs.
DMU
Target DMU
1
4
7
12
15
17
20
1
1
0.49504029
0.53660363
0.67987789
0.50878889
0.74467050
0.36909990
4
0.95525399
1
0.81873718
0.84047040
0.89121676
0.92400072
0.94263179
7
0.68426795
0.80156597
1
0.98769702
0.85624403
0.85682199
0.88697995
12
0.83162917
0.71092883
0.87483547
1
0.75165252
0.91722504
0.68657179
15
0.97375065
0.91292479
0.92595846
0.89933688
1
0.84347305
0.93897325
17
0.94554521
0.41786480
0.78137007
0.97194231
0.54846678
1
0.25490089
20
0.33804338
0.65282037
0.79256331
0.62361928
0.68029207
0.49238820
1Table 5
Example1: pairwise comparison matrices of seven CCR-efficient DMUs.
DMU1
DMU4
DMU7
DMU12
DMU15
DMU17
DMU20
Weight
DMU1
1
0.7646271
0.9123273
0.9171495
0.7644273
0.8967514
1.0232103
0.33493
DMU4
1.3078269
1
1.0095312
1.0757141
0.9886519
1.3569705
1.1753435
0.42144
DMU7
1.0960978
0.9905587
1
1.0601981
0.9638027
1.0423561
1.0526712
0.38601
DMU12
1.0903347
0.9296149
0.9432199
1
0.9222442
0.9722520
1.0387729
0.36912
DMU15
1.3081688
1.0114783
1.0375567
1.0843114
1
1.1905150
1.1539501
0.41546
DMU17
1.115136
0.7369356
0.9593650
1.0285398
0.8399725
1
0.8408676
0.34737
DMU20
0.9773161
0.8508150
0.9499641
0.9626742
0.8665885
1.1892478
1
0.36284Table 6
Example1: results for comparing.
DMU
k
Proposed method
Model (6)
Model (6) (AP)
Model (9)  l1-norm
Weight
Rank
Ω
b
Rank
θ
*
Rank
Γ
(
X
,
Y
)
Rank
1
0.33493
7
0.675
6
1.1009
7
0.0156
7
4
0.42144
1
0.708
2
1.9336
2
0.4802
2
7
0.38601
3
0.698
3
1.1737
5
0.1915
3
12
0.36912
4
0.677
5
1.1095
6
0.0589
6
15
0.41546
2
0.790
1
4.9139
1
0.7965
1
17
0.34737
6
0.677
5
1.3484
3
0.1760
4
20
0.36284
5
0.680
4
1.1849
4
0.1077
5Example 2.
Consider a performance assessment problem investigated by Sinuany-Stern et al. [17], which includes eight DMUs to be evaluated in light of two inputs and two outputs. The data set is provided in Table 7, together with the CCR-efficiencies of the eight DMUs, where the CCR efficiencies evaluate DMUs A, B, and D to be all DEA efficient and they need to be further differentiated.
Table8 shows the results obtained using the ranking methods of Alirezaee and Sani [20] and Sinuany-Stern et al. [17] and the proposed method. By our method DMU A has the best performance. In contrast, DMU A has the worst performance by Alirezaee’s method. In view of inputs and outputs of the efficient DMUs, it seems that our method gives more reasonable result than Alirezaee’s method.Table 7
Example2: DMUs’ data (extracted from [17]).
DMU
Input
Output
CCR efficiency
A
20
12
60
36
1
B
10
15
30
45
1
C
15
12
30
36
0.999
D
5
70
15
80
1
E
3
9
3
9
0.510
F
9
18
1
18
0.402
G
63
19
8
19
0.330
H
22
73
1
3
0.022Table 8
Example2: the results of ranking methods.
Sinuany’s method
Alirezaee’s method
Proposed method
B
0.20994
B
0.1993
A
0.591675
A
0.153001
D
0.1929
B
0.582900
D
0.14140
A
0.1883
D
0.556908


## 7. Conclusion

In this paper we propose a two-stage method for ranking all extreme and nonextreme efficient DMUs based on cross efficiency and AHP technique. In the first stage, cross efficiency score of each DEA-efficient DMU is calculated. In general, the optimal input-output weights obtained by classical DEA are not unique; therefore, the cross efficiency scores depending on these weights are also not unique. To remove this problem we first obtain all alternative optimal solutions of DEA model (5) and then the arithmetic mean of them is considered as a secondary goal. In the second stage, using a one-step process of AHP, the DEA-efficient DMUs are ranked. It seems that the proposed method gives more reasonable results, because we consider all alternative weights corresponding to DEA-efficient DMUs. Also, our method can be extended to interval cross efficiency and interval AHP. We suggest as future works a deeper analysis in this subject. Finally, in this paper we rank CCR-efficient DMUs; in a similar way one can also rank BCC-efficient DMUs.



---

*Source: 594727-2015-01-05.xml*