(FDH)and - J-Stage

6 downloads 0 Views 897KB Size Report
Sep 20, 1994 - Part II:Iih]eeDisposal Hull (FDH)and. Russell Measure ... dominance as treat,edby Free DisposalHull and ..... P inte R. Except inthe special case.
The OpeiationsReseaich The Operations Research Society Society

of of Japan Japan

Journal of the Operations Research Societ.v of Japan Vel. 39, No. 3, September 1996

MODELS AND MEASURES Part II:Iih]eeDisposal Hull

FOR

EFFICIENCY DOMINANCE IN DEA (FDH)and Russell Measure (RM) Approaches'

W. F. Bowlin University of Northern Iowa

1. Bardhan fhe University of

Tkxas at Austin

W. W. Cooper T7ieUnit:ersity of 7'exasat Austin

T. Sueyoshi ScienccUniversity of fbkyo

September 20,1994; RevisedApril12,1995) {Received Abstfuct Models and measures of eMc:iency Measure approaches te eMciency evaluatioR of EraciencyDominance) in DEA. (Measllres

dominance are

examined

as as

treat,edby Free DisposalHull and Russell they relate te additive medels and MED

1. Introduction Mixed integerversions of the additive associated measures "rere intromodels of DEA with duced in [4] as a way to dealwith dominance". Other models and measures are exainined in this paper. [I]hefocusison (a)the DisposalHull" (FDH)and (b)the "Russell Measure" (RM)approaches, euch of which has formed the basisformajor research efforts the first by H. M. Tulkens and his associates at the University of Louvain and the second by C. A. K, Lovelland hisassociates at the Universityof Georgia.See,e,g., [11] and See also [22j-[25]. l13]--[16]. "eMciency

C`Free

2. FDH

Ii}rbee Disposal Hull VLrestart with FDH after definingXJ-,and Ylas input and output vectors, respectively, with components xiJ・, i 1,・・・,m and y,j,r 1,・・・,s,for each of a collection of DMUj -・ ・ Making where Then, we use DMU, to designate the DMUj (Decision Units) j' 1,・ to be evaluated and sa}r that the eMciency ef itsobserved is dominated by performance DMUk ifX. 2 Xk and }'bS }'li with strict inequa!ity holding!n at leastone inputor output ==

=

=

,n.

component, "free As noted in Bardhan et al. [4], the assumptiQn of disposal"givesrise to a danger of identifying a dominatedDMU as efflcient when non-zero slack ispresent. This dangeris avoided in the FDH approach, however,by a first-stage algorithm which uses paired vector comparisons to identify nondominated optima DMUs.i The possible presence of alternate can nonetheless lead to problems likethose we now describe; Fbr simplicity we followLovell[16, pp. 38-40] and restrict attentioll to output・s by assuming that all DMUs use only a single input in the same amount. Then, using yrj to "Acknowledgment

versions

of

this paper

isowed to R. M. Thrall fornumereus comments and suggestions he offered on earlier itspredecessors. Support forthisresearch from t,heIC2 Institute ofThe University

and

Texas at Austin isalso gratefu11yacknowledged. ISee [[hilkens and Vanden Eeckaut [24, between input and pp. 3-5]who furtherdistinguish "weak dominanee--and introduceadditional categories such as dominance;;,etc,, which we do sue herebecausethese refinements do not Ieadto importantdifferences in our results, of

output not

pur-

333

NII-Electronic Library Service

The OpeiationsReseaich The Operations Research Society Society

of of Japan Japan

334

I,Barthan,LV.E Bo-clin, LM W. (lbolper & T.Sucyoshi

represent

output

r

forDMVj, the

ip* m,ax.=m,ln,. ==

Here J' E D

eMciency

( £

measure

in FDH isobtained

used

)L, 1, A, E

;E!i,l"')Lj

=

{o,1},vv E

D)

from

==

,Z,D

;:i (1)

that the indexesare selected from the set D of non-domina,t;ed DMUs, as determinedin the first stage algorithm, so tha,ty.kly..2 1, Vr,with equality achieved when the DMU. to be evaluated iseMcient-as indicated by the choice of Aj 1 from the means

==

A,

vector

The followingexarnple

can

help to show

what

isinvolved,

di

max

subject

to

(2)

12ip f{l4Ai+24)L2, 12gbSl 15Ai+14)L2, 12ipE{ l4A,+14A,, l == A,+A2,

Ai and A2 must b,eeither zero or ene,2 We then choose between DMUi and DMU2 "rhich both dominate DMU., the DMU to be evaluated, with 12 units recorded foreach of its three outputs. To put this in the form of (1), we divideeach constralnt in (2) by the corresponding output of DMU. and obtain, where

A:

ips{

ll/IAi+?l/1 Il lg/1)Li+ IS){2,

The

measure

==

qs s{

L2,

1 and

)Li+tl/i)L2, qsf{

defined by (1) isassociated 1 both give

)L: =

with

(3)

i=)Li+)L2J

alternate

optima

fbr(3) since

the

cheices

14

(4)

ip: ipl=Tt, =:

Thisvalue of ip* also givesthe proportioninwhich all of DMU.'s to be augmented. With the choice of either DMUi or DMU2 as the evaluator we thus findthat DMU. is inerncient-i,e.,ip' > 1-----and this same value of ip' indicatesthat DMU. shouid have produced 14 units in each of itsthree outputs. Althougheach of these alternate optima yield the same value of di*, the two solutions are associated with amounts of slack. This isdealtwith by Tulkens et al. by quitedifferent simply listing the slacks but this isinadequatewhen alternate optima are present (as in the ifip* is to be used forranking, In addition, simple listing of non-zero slacks above example) isalso questionablewhen alternate optima may be present. Seethe example in Bardhan et aL [4]. Nor isthis the end of possible troublesfor,as can be seen, we need to only replace the output values of 14 units by 12 units in (2) to obtain a value of ip' 1 which accords DMU. The same measure of efficiency as isaccorded to the two DMUs which dominate it. Figure 1 below will help to illuminatewhat isinvolved,The solid Iinerepresents the to as the Free DisposalHu]1, This FDH is derivedfrom graph of a step functionreferred the observed values of yi and fbr each DMIJ via the followingcensiderations, All points y2 on or below the solid lineto the left of R are dominated by the observations recorded at R as yi 8,y2 2g. This solid lineisnot extrapolated to the vertical axis, however,because Q with yi 4 and y2 3 isnot dominated by R, or vice versa and so on, PointsQ and R, which form the vertices firom which thi$ FDH isgenerated, are obtained from pairedcomparisons in the firststage. These vertices correspond to the set D of nondominatedvectors used in (1). to this set, the point P will have a value di Relative as

t・hemeasure

outputs

of efficiency.

are

=

=

=

=

==

=

2The measure

non-dominated of

DMUs

assigned

to D in the first pass

ef

the

algorithrn

are

all

given an

i

erncienc}'

unity.

NII-Electronic Library Service

The OpeiationsReseaichSociety The Operations Research

Society of ofJapanJapan

ET}i'ciency Dominancc-Ibec DisposallllJil

335

Y2

)

Yl

Figure 1:FDHPertraval u therefore be ineficient. This same

;

;(2,2)

value of ip' isused to obtain (3,3), P'. This projectionto P' therefore results in a point which is on the Ftee DisposalHull.However, the slack value of sl 4- 3 1 isunaccounted fbrin this measure of eficiency, Under the assumption of disposal"this slack might be ignored butT'ulkensand hiscolleagues to at leastlist these values. However, they generally proceed have been unable te includethem in a more comprehensive measure--as Tulkens notes [23, 15]: doubt tha,t seeking fora single numerical index of eficiency which p, includesthe slacks would be a usefuI extension, Additionalcolumns exhibiting the values of a}1 slacks would be more informative".Such a simple listingcreates diMcultieswith Tespect to the rankings that Tulkens and hisassociates use, however, and issuescan also be ra-ised with respect- to alternate optima with differing amounts of slack while, finally, problems alse arise becausea value ef e' 1 may be assigned to a dominated DMU.

a,nd

the

coordinates

=

of

=

=

==

"free

"I

=

3."leMID

and

approach

MED

Measures these problems via the following model

S} i=! Xio

÷

aL

[4],

2 r=1SrYro

n

to:

subject

Bardhan et

s

s,+

max

taken from

xi.

Aj+si-

=Exij

,

i=1,・・・,m,

j:,1

yro ==

2 yrjAj- s.-,

(5) r=l?-・・,s,

j;tl

1

-:

£ Aj. j'-m1

"iethegaln some evalua,tor

.

perspectlve ofDMU.

as

by on the case where DMUk isused inthis the fbllowingfocusing expression forthe ith input and rth output

on

model

as of

NII-Electronic Library Service

The Opeiations The Operations Research Reseaich Society Society

of of Japan Japan

336

LBaidhan, LV.Pl Bowiin,;U UL Cbqpcr&T Sucyoshi

DMUo,

-*

+*

S,

+S. Yro

Xio

The first term

(6)

=Xio-Xik+Yrh-Yro,

Yro

Xio

the input excess

input used and the second represent the outpvt shoftfa1} relative to the out・put produced by DMU. compared to DMUk. Both terms are stated in limeasure, DMU. to be as istrue for all the other evaluated in (5). Itistherefore meaningfu1 to sum these measures and interpret the result as a. measure of relative distance firomwhat isneeded to arrive at the fullDEA eficiency. This same measure of relative distanceisused for each DMUj in (5).IUIIeMciency for any one of these DMUj can be achieved only by traversing this distance,A ranking based on thismeasure of relative distancecan thereforebe used. There are, of course, other criteria that can be used for ranking such as tot・alcost, total profitab!lity, etc., but the above measure, however,isinvariantto the measures used in the =i. and y... SeeBardhan represents

to the

relative

amount

of

[4].

et al.

The rankings obtained from (5) in thisfashionmay differ from thoseobtained from (1). The latter with a value of 1 S di'. A normed measure can also be yields a radial measure derivedforthe inefficiency of DMU. relative to the eutput of DMUk by focusingonly on the output terms in (6) and writing '

-*

l+ To

possiblevalues

avoid

exceeding

Yrk

Sr

=-.

Yro

unity,

we

Yro the

use

Yro

unity

and

average

measure

which

over

all

may

of

relative

to

for output

eficlency

to

obtain

r

reciprocal

forDMU. which cannot exceed Qf ip' in (1), Moreover, we can

obtain

£ o'< -

r-iYrolYrk g1

'

(7)

s

ifDMU. isfu11y eficient, than We illustrate with P in Figure l where we observe that (5) designates R rather for its evaluation, as indicated by line from P to S and then from S the dotted Q going to R---torepresent as the relative the tidistancemeasure that isused,3 This gives measure of eMciency in the licomponent foryi and as the eficiency measure fory2, also in limeasure, and the average of these two ratios is with

unity

attained

ifand

this expression '

Yrk

be interpretedas the

also

outputs

of

.- Yre

Yro+Sr-'

We then have a

reciprocal

only

g t =

tilt; g =

E+g 2

This provides the nance) in Bardhan et

measure aL

which

with [4]

was

3Also

the discussion of this and called

`Ccity-block7'

other

metric

40

labeledas MED

a complementary

21 ]・--=40 '

21

in the

operations

measure

19 40

of EMciencyDomi(=Measure of

inefllciency given by

'

research

fbra literature. See Appendix A in [6]

metrics.

NII-Electronic Library Service

The OpeiationsReseaich The Operations Research Society Society

of of Japan Japan

337

Efl)'cicncy DominEince-FteeDisposnJ HbN which

was

to as MID

referred

of lneficiency Dominance). (Measure

These meftsures are obtained from the DMU which ismost dominant in this timetric, by the choice of R rath ¢ r than Q in Figure L However, these measures do not result in a scalar which will projectP inte R. Except in the special case of fullyeMcient to the components in the MED or performance, such identificationswill require recourse MID measures. as exhibited

4. RussellMeasures C, A. K. Lovell[l6, p. 29] has expressed measures.

the foIIowingconcerns

with

some

of

the

available

'

de

"I

not

combining

likethe ideaof slacks

aggregating

witharadial

slacks,

eMciency

I likeeven

and

lessthe idea

of

'

score,"

Possiblyforthese reasons,4 he turned (with Friedand others) to a use of "Russell measures in a series of papers evaluating the performance ef U.S. Credit Unions. We therefdre investigate this approach and relate itto our precedingdiscussions via the fo11owing mode] "5

m

s

2

max

iprtO

SUbjeCt

£ et i=1

r=1

n

! Z)YrjA,i, iprYro

1i''',s,

rb =

(8)

j:i

2xi,iAj,

eixi.}l "-1i 1

i

==

1,・・・,m,

=EAj

j'=1

A.iG {O, 1}, Vi,and add the m + s additional conditions O f{lei E{ 1, 1 S ip, to ensure that the requirements for dominance are satisfied. Unlikethe situat・ion for (5), however,it isnow necessary to assume that all data are positivesince the possibi}ityof infinite solutions must otherwise be admitted because recourse to the Ali-Seiford translation discussedin Bardhan et al, [4] is not available forthese ¢ and ei values, The objective in the above formulationdiffers from the one in Fare,Grosskopfand Lovell as well as the one in FEre and Lovell However, itsuits our immediatepurposes [12] [13].6 to treat these other formulationsof the above objective laterin this paper. We can t・hen begin our comparisons herewith the above effciency measures and models by presenting an optimal solution to (8) in the fo11owingform where

we

also

require

.

di:yro=yr" r==

when

AX

=

1 isan

optimal

Yro=Yrk-Sr-', when 4In

An

=

1 isan

optimal

1,・・・,s

choice, r=1i'''is

choice

and

Similarlywe and

el-xio xih ==

i= 1,,・.,m,

(9)

write

xio=xik+s"・*,

i=l7・・・,7n,

(10)

for(5).

a subsequent clevelop a measure paper, Lovelland Pastor[17] in the concluding sections of Bardhan et al. [4], 5See Fare and Lovell [13]. See also Russell[21], 6See also Russell[21].

which

isanalogous

to the ones

'

described

NII-Electronic Library Service

The OpeiationsReseaich The Operations Research Society Society

of ofJapanJapan

338

l.Bardhan,iMEBowlin, ;U iM Cbopet'&TlSuqyoshi

We

now

also

(9)to definenew

use

--

SF sl

the terms

where

the

on

variables

non-negative

ip:Yro-Yro=Yrl-Yro,r=

right

variables

are

Yrh7

=

=

eiXio

We (9).

from

obtained

iprYro Yro +STL'

r

(11)

to introducethe new (10)

use

also

1)'''7si

=

(12) (5).

i-- 1,...,m,

xi.-s,+・'=xik,

:=

1,''',s,

i=1,.・・im,

=xi.-e,'・xi.=xi.-xit,

O g e,S 1, 1 g ip. so satisfy (8)issatisfied and, similarly, the expressions (11) From (11) we see that these s.- and s,+・ are associated with Ai 1 as a solution to (5), forwhich Al 1 isoptimai, while from (12) we see that these ip. and eiare associated with Ak fbrwhich Ar isoptimal. Hence suppose we could have fbr(8), where

=

=

=1

=1

s

m

£ di;

s

-2e,'・

r=1

the choices

with

Ar

=

1

Ak

and

=

>

i--1

m

-Ee, di.

£

(i3)

i=1

ri=1

1. Using (ll) and

however,this would (12),

give,

(i4) t/l.,ii/i.+l..,iift>t9.,; the

contradicts

which

direction, assume

optimality

that the

for A#

assumed

choices

A2

1 and At

==

=1

in give

:1

(5).Turning in the

opposite

(is, tf.ll,;'ll+:.,/fi>t9.,:/l for(5). However, again,

using

and we ebtain (11) (12) m

s

s

Edir£ ei>Zip.' -

iul

r=1

J

the optimality

contradict・s

which

Theorem By

A

1.

virtue

of

this theorem

we

can

the

AI,as

used

of

and

the

other

and

alternate

(5)

are

optima

DMUo・

ofit is optimal fbr(5).

to the

(17)

i=: 1,・・-,m,

same

optimai

aTe

present. As

(17)

may

One minimizes eficiency be present in evaluating

complementary,

ineficiencies that

maximizes

oniy

xio-e,{xie=s,+-*i

are applicable ip:,e,*-

the

aboverwhen

in (8) and

objectives a,nd

slacks

ifand

therefore have the foIlowing

write

di:yro-yro=s.LX r=17・・・,s where

fbrAr in (8). We

assumed

(16)

e,i i--1

r=1

to (8) isoptimal

solution

m

T2

now

choice-either makes

clear,

Al or the two

accomplishments

the performance

With (17) in hand, we

ca.n do more than has heretofore been done with Russell MeaConsider,forinstance,how one might treat variables, Various approaches are describedin Banker et al, 152-153]. Here we focus only on the [2,pp. Banker-Morey [3] approach rather, itsextension intothe presentcase) foradded insight intothe assumption of disposal".To adapt this approach to a use of Russell measures, ene simply assigns values of unity to the pertinent ip. or ei and carries these values "nondiscretionary"

sures.

"free

(or,

NII-Electronic Library Service

The OpeiationsReseaich The Operations Research Society Society

of of Japan Japan

ErficiencJt Dominance-Jin]eDisposafHljli

intowhatever

339

isused. As describedin Cooper et al. [10], the ei are then interpretable as di. goods" which do not enter intothe overall eMciency scores except insofaras they can improve the values of other in the efliciency scores. Turning to (5), on the other hand,we reca.11 that ipJor e,"・ the slack values in the objective (but not in the constraints) enter intoour eMciency scores. The zeros associated with (17) under a choice of e,'・ or ip: 1 are then interpret・ed to mean that・this is the value te be used in MED or the other ratings and evaluations discussedin This isachieved by assigning a zere coeflicient to these slacks in the objective. Hence, [41. in addition to being a freegood in the constraints, non-zero slacks are also to be regarded as not involving any cost in their disposalin a manner analogous to th¢ case when free disposalisassumed. Puttingthis all together, this study can see that we have genera}izedthe Banker-Morey treatment of non-discretionary slacks t-ocases in which some not all) of the inputsand (but outputs are non-discretionary. This applies to the BCC' and CCR models used by Banker and Morey and extends to other models as well. Indeed,as describedin Banker et al. [2], even more situations involving and generaltreatments are available to accomniodate "ceilings" iNrith ranges fbr which some of the variables are discretionary and other ranges non-zero

aggregate

slacks

eMciency

associated

measure

these

with

"free

or

=

"fioors"

they

when

We the

are nondiscretionary. inake

now

our

form

promised

I.i k

ir

e +

minimize and

a

awkward

'n

objective

used

0:' E]/'=i

S

in Fare

et at.

which [12],

has

i/di;. Z)"T=,i

¢

os) + , that is (a) bounded by unity and (b)attains this beund provides a measure if DMU. is eMcient, As seen on the right of (18), the FGL objective is to weight・ed average of arithmetic and harmonic means. It isdifficult to interpret to use so we repiaced itby the simpler objective in (8), and we now justify ,.i.

This objective if and only

to the

return

this choice

by the foIlowing,

Theorem

2.

(18),

An

optimat

=

+]Ii

solutien

,.+

to (8)isalso

eptimal

when

the

objective

isreplaced

[[bprove t-histheorem, lete,'・,ip: be optimal for (8) and letbi,di, be optimtt1 in (8) isreplaced with Now assume that we could have (18).

when

6y the

objective

£

(ig)

.,bi+te.,S. O, Vr, so,

a

fortiori, m

Ebi i=:1

s

-Z&. ・r=1

s

m

- 1, Vr. All + 11ip.* ip. di.' the assumtion that these O,'・ and ipIwere where

::::

must

(19)

replace

with

di,

Hewever, bi, were

holdand

being satisfied, we havea contradiction for (8), avoid this contradiction,

constraints

so, as

to

rlb

optimal

we

(23) #-i":'t9-ii*r`i=iO te be

assumed

the theorem

optimal

asserts,

the

under

Hence, equality (18).

objective

have

we

must

S.,ei']l.li,2・= £

(2`)

.,bi't9.,21':,'

By

of

virtue

two theorems,

above

(fA:

1.

Corollary

the

have the foilowillg,

alse

then fbr(5),

1 is optimat

=

we

it is atso optimal

for(8)with (18)as

its

objective.

from the preceding developments. As back and forth, or we may directly by (17) proceed setting e,'・ xiklxi. and ip.'y.kfy,.,etc. In any case, (8), and iseasier to compute simpler to int・erpret.7It is also more general.For instance,iteliminates the need forthe rn + s additional constraints required foreach DMU. to be evaluated in (6)to insurethat dominance isnot violated, and there isno need te assume that all input-s and outputs are for every DMU, as is the case for Finally, the use of us in contact with positive (8). (5)puts the even more general model" of DEA with to still are related properties which other DEA models, as describedin Ahn, Charnes and Cooper [1], and this makes itpessib}e to exploit featuresand extension of these models-which iswhat we did when we showed how to handlenon-discretionary inputsand outputs in eur discussionof (17). [[bt-akefulladvantage of what has just been said, we complete the present development by removing the restrict・ion to integerselutions so that we can extend our results to the additive model as first With thisrestriction removed general given in Charneset al, [9]. from both (5)and (8), we then have We

are

already

in position to reap

now

noted,

can

we

=

several

to

use

advantages

move

=

LLadditive

Corollary 2. Let Ak reforesent Xu・tion to (li?. (fAZisoptimalfor

for(8)using

optimal

A

the

a

(5) with (18)

objective

by

so-

then it' isalso

removed,

the integrutityrestriction'removed.

n

ZxiJA;

A:.+si" }il

=2x,j

formedfroma

components

writing

n

x,.

non-negative

the integratity restriction

with

proofcan be obtained

simple

with

vector

,,:i

yro =EYriAj'

i= 1,''',m?

=m:k,

(25)

J:i

-s.-*

g2yrp AJ

j=1

=

yr*ki

'=

1,'''iS'

j=1

We then have '

Xio-Xlk=S,+' 7Fare,

i= 1i'''iM

and Grosskopf Lovell[13] also

(only)

objective

Reference to slacks,

';

which

which

(]7),however, means

simplify

matters

the ei and an eutput shows that this iseqnivalent that these approaches are suited minimizes

y.'k-yre -- s,r'i

and

fbrsome

(only)

applications

oriented

objective

1,'''is・

r=

by

using

which

an

input・ oriented

maximizes

to maximizing only the inpllt slacks only to very special situations.

(26)

or

the

Or.

the output

NII-Electronic Library Service

The OpeiationsReseaich The Operations Research Society Society

of of Japan Japan

EfficiencyDon]i"ancc--P}ecDisposafHuli

We

have

also

eixio==//k, i-- 1,・・・,m The

341

(27)

r=1,-・・,s. ip:y.'.=y.'k,

and

the proof followsthe same route as before.That is,we show that these e,*・ and ip: are optimal for (8)and that they remain optimal when the objective isreplaced by We can then conclude this extension of the precedingdevelopmentsby noting that the same route can be fo11owedto either of the two classes of additive models by omitting remainder

of

(18).

that the solutions must satisfy. not intendedto mean that further research on uses of Russellmeasures should be abandoned. Indeed,such research can buildon what has already beenaccomplished. The choice of suitable measures since the measure providesan example in is not the only Others inight be prescribed explored as follows, (18) possibility, The component sums in the objective fbr(8) sat・isfy or

imposmg

a convexity

constraint

The preceding developments are

'

E]Y-i ipr s

su, combining

the two,

we

have

O, EI);・II , Sl,

)1'

(28)

.

sli 3,: M

(,,)