Ranking Efficient Decision Making Units (DMUs) are an important issue in Data Envelopment Analysis (DEA). This is one of the main areas for the researcher. Different methods for this purpose have been suggested. Appearing nonzero slack in optimal solution makes the method problematic. In this paper, we modify the nonradial supper efficiency model to remove this difficulty. Some numerical examples are solved by modified model.
1. Introduction
Data Envelopment Analysis is a mathematical programming technique which evaluates the relative efficiency of DMUs. DEA classifies the DMUs into two different classes, called set of efficient DMUs and the set of inefficient DMUs.
Conventional DEA model cannot differentiate the efficient DMUs whose efficiency value is one. Toward this end, different methods are suggested; see [1–3]. One of the important models is AP-Model which was proposed by Andersen and Petersen [4] and also see [5, 6].
This model is widely used and the results are almost satisfactory. The main deficiencies of AP model are being as follows:
infeasible for some kind of data,
unstable in the sense that a small variation in data causes big increase (degrease) in the result,
not taking into account the nonzero slacks which appear in optional solution (projection of omitted DMU is weak efficient in new PPS).
For removing these difficulties, many researches have suggested different models; for example, see [7–9].
Tone [7] suggested the nonradial supper efficiency model. This model fails to remove the 3rd deficiency. In this paper, we modify the nonradial supper efficiency model that takes into account nonzero slack that appears in optimal solution.
The rest of the paper is organized as follows. In Sections 2 and 3, AP-model and Tone’s model are discussed. Section 4 contains the modified model. In Sections 5 and 6, input and output oriented models are proposed. Discussion and conclusion cover Section 7.
2. Anderson Peterson (AP) Model
Consider n decision making units DMUj(j=1,…,n) which consumes 0≠xj≥0 and xj∈Rm(j=1,…,n) vector as input to produce output vector 0≠yj≥0 and yj∈Rs(j=1,…,n). The supper efficiency model may be written as follows:
(1)θp*=Minθs.t.∑j=1,j≠pnλjxij+si-=θxipi=1,…,m111∑j=1,j≠pnλjyrj-sr+=yrpr=1,…,s111λj≥0j=1,…,n,j≠p111si-≥0i=1,…,m111sr+≥0r=1,…,s.
The following example which has been taken from [7] shows the deficiency in case of having nonzero slacks in optimal solution.
Example 1.
Consider the data given by Table 1.
DMU
A
B
C
D
E
F
Input 1
2
2
5
10
10
3.5
Input 2
12
8
5
4
6
6.5
Output 1
1
1
1
1
1
1
By using AP-Model,
(2)SupA=θA*=1S2-*=4SupB=θB*=1.260AllslacksarezeroSupC=θC*=1.133AllslacksarezeroSupD=θD*=1.250S1-*=7.5SupE=θE*=0.750SupF=θF*=1.
In other words, the supper efficiencies are as follows. (3)θA*=1-4ɛ(*)θB*=1.260θC*=1.133θD*=1.25-7.5ɛ(*)θE*=0.750θF*=1(*).
Note. These values are not correct in [7], because the result shown in the paper has taken from original paper in which F is not included. From Figure 1 it can be seen that the radial supper Efficiency model is not able to rank F which is non-extreme efficient and also is not able to rank D and A, with nonzero slack in optimal solutions.
3. Nonradial Supper Efficiency Model (Tone’s Model)
In 2002, Tone proposed the following nonradial supper efficiency model:
(4)ρp*=Min(1/m)∑i=1m(x¯i/xip)(1/s)∑r=1s(y¯r/yrp)s.t.∑j=1,j≠pnλjxij≤x¯ii=1,…,m11∑j=1,j≠pnλjyrj≥y¯rr=1,…,s111x¯i≥xipi=1,…,m111y¯r≤yrpr=1,…,s111y¯r≥0r=1,…,s111λj≥0j=1,…,n,j≠p.
Based on the SBM model (4), the nonradial supper efficiency model should be as follows:
(5)γp*=Min1+(1/m)∑i=1m(si+/xip)1-(1/s)∑r=1s(sr-/yrp)s.t.∑j=1,j≠pnλjxij-si+=xipi=1,…,m111∑j=1,j≠pnλjyrj+sr-=yrpr=1,…,s111λj≥0j=1,…,n,j≠p111si+≥0i=1,…,m111sr-≥0r=1,…,s.
We will show that (4) and (5) are not equivalent, in the sense that their optimal solutions are different.
Example 2.
Consider the eight Decision Making Units (DMU), with two inputs and one output (see Table 2).
DMU
A
B
C
D
E
F
G
H
Input 1
4
7
8
4
2
10
12
10
Input 2
3
3
1
2
4
1
1
1.5
Output 1
1
1
1
1
1
1
1
1
First the following are defined:
(6)PPS={(x,y)∣x≥∑j=1nλjxj,y≤∑j=1nλjyj,λj≥0,j=1,…,n},PPS¯={(x,y)∣x≥∑j=1,j≠pnλjxj,y≤∑j=1,j≠pnλjyj,λj≥0,j=1,…,n},PPS¯¯=PPS¯∩{(x,y)x≥xp,y≤yp}.
In above-mentioned example and shown in Figure 2,
(7)PPS=REDCK=SPPS¯=R′E′DCK=S¯R′E′K′¯¯=S¯¯.
We can see that S¯¯ is a proper subset of S¯ and S¯ is a proper subset of S; that is,
(8)S¯¯⊂S¯⊂S.
In this example, the optimum value of objective function in (4) is
(9)ρE*=12(42+44)=1.5
and the optimal value of objective function in (5) is
(10)γS¯*=12(42+24)=12(2+12)=54=1.25,
so (4) and (5) are not equivalent. In other words, the constraints x¯≥xp and y¯≤yp and y¯≥0 should be omitted from (4). Now we show that the nonradial supper efficiency model has the same difficulty as radial supper efficiency model in treating nonzero slacks in optimal solution. Using (5), the supper efficiency of E, D, C may be evaluated as follows:
(11)γE*=Min(1+12∑i=12si+xiE)s.t.4λA+7λB+8λC+4λD+10λF+12λG+10λH-s1+=23λA+3λB+λC+2λD+λF+λG+1.5λH-s2+=4λA+λB+λC+λD+λF+λG+λH+s-=1λA,λB,λC,λD,λF,λG,λH,s1+,s2+,s-≥0,γE*=1+12(22+04)=1.5.
The projection of E is on weak frontier
(12)γD*=Min(1+12∑i=12si+xiD)s.t.4λA+7λB+8λC+2λE+10λF+12λG+10λH-s1+=23λA+3λB+λC+4λE+λF+λG+1.5λH-s2+=4λA+λB+λC+λE+λF+λG+λH+s-=1λA,λB,λC,λE,λF,λG,λH,s1+,s2+,s-≥0,γD*=1+12(24+03)=1.25.
The projection is on strong frontier
(13)γC*=Min(1+12∑i=12si+xiC)s.t.4λA+7λB+4λD+2λE+10λF+12λG+10λH-s1+=83λA+3λB+2λD+4λE+λF+λG+1.5λH-s2+=1λA+λB+λD+λE+λF+λG+λH+s-=1λA,λB,λD,λE,λF,λG,λH,s1+,s2+,s-≥0,γC*=1+12(28)=1.125.
The projection is on strong frontier.
In summery,
(14)γC*=1.125γD*=1.25γE*=1.5.
In the case the projection of omitted DMUp lied on weak frontier, nonzero slacks appear in optimal solution and the following approach is suggested.
4. Modified Nonradial Supper Efficiency Model
First solve the following model:
(15)γp*=Min1+(1/m)∑i=1m(si+/xip)1-(1/s)∑r=1s(sr-/yrp)s.t.∑j=1,j≠pnλjxij-si+=xipi=1,…,m∑j=1,j≠pnλjyrj+sr-=yrpr=1,…,sλj≥0j=1,…,n,j≠psi+≥0i=1,…,msr-≥0r=1,…,s.
The model, (15), can be linearized by the suggested method in [7] and solve by the simplex method.
Suppose that
(16)(x^,y^)=(xp+s+*,yp-s-*)
is the projection of DMUp on the frontier of PPS¯ and now solve the following model
(17)Wp*=Max(1/m)∑i=1m(s^i-/x^i)(1/s)∑r=1s(s^r+/y^r)s.t.∑j=1,j≠pnλjxij+s^i-=x^ii=1,…,m111∑j=1,j≠pnλjyrj-s^r+=y^rr=1,…,s111s^r+≥0r=1,…,s111s^i-≥0i=1,…,m,111λj≥0j=1,…,n,j≠p.
Set
(18)μp*=γp*-Wp*(supperefficiencyofDMUp).
It is evident that if (x^,y^) is strongly efficient, then Wp*=0 and μp* is supper efficiency score of DMUp.
5. Modified Input Oriented Nonradial Supper Efficiency Model
First the following model is solved:
(19)ξp*=Min(1+1m∑i=1msi+xip)s.t.∑j=1,j≠pnλjxij-si+=xipi=1,…,m1111111∑j=1,j≠pnλjyrj-sr+′=yrpr=1,…,s1111111si+≥0i=1,…,m1111111sr+′≥0r=1,…,s1111111λj≥0j=1,…,n,j≠p.
Let
(20)(x~,y~)=(xp+s+*,yp+s+*′).
Now solve the following problem. (21)Wp*=Max(1/m)∑i=1m(ti-/x~i)(1/s)∑r=1s(tr+/y~r)s.t.∑j=1,j≠pnλjxij+ti-=x~ii=1,…,m11111∑j=1,j≠pnλjyrj-tr+=y~rr=1,…,s11111ti-≥0i=1,…,m11111tr+≥0r=1,…,s11111λj≥0j=1,…,n,j≠p.
Let
(22)μp*=ξp*-Wp*.
6. Modified Output Oriented Nonradial Supper Efficiency Model
Consider the following model:
(23)ηp*=Min(11-(1/s)∑i=1s(sr-/yrp))s.t.∑j=1,j≠pnλjxij+si-′=xipi=1,…,m∑j=1,j≠pnλjyrj+sr-=yrpr=1,…,ssi-′≥0i=1,…,msr-≥0r=1,…,sλj≥0j=1,…,n,j≠p.
Let
(24)DMUpprojection=(x¯,y¯)=(xp-s-′*,yp-s-*).
Now solve the following problem:
(25)Wp*=Max(1/m)∑i=1m(s^i-/x¯i)(1/s)∑r=1s(s^r+/y¯r)s.t.∑j=1,j≠pnλjxij+si-=x¯ii=1,…,m∑j=1,j≠pnλjyrj-sr+=y¯rr=1,…,ssi-≥0i=1,…,msr+≥0r=1,…,sλj≥0j=1,…,n,j≠p.
Let
(26)ρp*=ηp*-Wp*.
7. Conclusion
In this paper, it has been shown that both AP-Model and nonradial supper efficiency model are not able to rank DMUs if the projection of omitted DMU is weak efficient in PPS¯. The new method removes this difficulties and the example which is solved by using the new method (modified method) confirms the validity.
The results for comparing the methods are shown in Table 3.
DMU
C
D
E
F
G
AP
1.14
1.25
2.00
1.00
1.00
SBM
1.125
1.25
1.50
1.00
1.00
Modified SBM
1.125
1.25
1.25
0.90
0.73
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
KhodabakhshiM.AryavashK.Ranking all units in data envelopment analysis201225122066207010.1016/j.aml.2012.04.019MR2967790ZBL1260.90113KhodabakhshiM.AsgharianM.GregoriouG. N.An input-oriented super-efficiency measure in stochastic data envelopment analysis: evaluating chief executive officers of US public banks and thrifts2010373209220972-s2.0-7044952542810.1016/j.eswa.2009.06.091JahanshahlooG. R.KhodabakhshiM.Using input-output orientation model for determining most productive scale size in DEA20031462-384985510.1016/S0096-3003(02)00636-7MR2008592ZBL1101.90366AndersenP.PetersenN. C.A produce for ranking efficient units in data envelopment analysis19933912611264KhodabakhshiM.A super-efficiency model based on improved outputs in data envelopment analysis2007184269570310.1016/j.amc.2006.06.110MR2294934ZBL1149.90077KhodabakhshiM.An output oriented super-efficiency measure in stochastic data envelopment analysis: considering Iranian electricity distribution companies20105846636712-s2.0-7795059367810.1016/j.cie.2010.01.009ToneK.A slacks-based measure of super-efficiency in data envelopment analysis20021431324110.1016/S0377-2217(01)00324-1MR1922620ZBL1073.90520KhodabakhshiM.Super-efficiency in stochastic data envelopment analysis: an input relaxation approach2011235164576458810.1016/j.cam.2010.03.023MR2805700ZBL05930981Ramezani-TarkhoraniS.KhodabakhshiM.MehrabianS.Nuri-BahmaniF.On ranking decision making units with common weights in DEAApplied Mathematical Modelling. In press