基于Shephard距离函数的DEA TOPSIS方法研究
第30卷 第5期运 筹 与 管 理
Vol.30,No.5
2021年5月OPERATIONSRESEARCHANDMANAGEMENTSCIENCE
May.2021
收稿日期:2019 05 18
基金项目:国家自然科学基金资助项目(72071192,71671172);国家自然科学基金重点项目(71631006);中科大创新体培育项目(WK2040000027);2020年安徽省级质量工程教学研究项目(2020jyxm2279);安徽省校企合作实践教育基地项目(2019sjjd02);中科大教学研究项目(2019xjyxm019,2020ycjg08)
作者简介:李勇军(1982 ),男,安徽无为人,副教授,硕士生导师,博士,研究方向:数据包络分析;蔡华权(1993 ),男,福建泉州人,硕士研究生,研究方向:数据包络分析。
基于Shephard距离函数的DEATOPSIS方法研究
李勇军, 蔡华权
(中国科学技术大学管理学院,安徽合肥230026)
摘 要:本文通过对Shephard距离函数的引入,正式构建了DEATOPSIS决策单元排序方法的框架。本文首先定义
了正(负)理想决策制定单元(DMU)以及相应的(反)生产可能集,然后在考虑正(负)理想DMU的条件下分别给
出D
MU的(反)效率评价模型以及对应的Shephard距离函数,然后基于评价对象到理想DMU相对接近度这一综合评价值给出了D
MU的一个完全排序。最后,本文通过算例分析说明了该方法的有效性和实用性。关键词:DEA;TOPSIS;排序;Shephard距离函数;生产可能集;效率前沿面中图分类号:C934 文章标识码:A 文章编号:1007 3221(2021)05 0116 06 doi:10.12005/orms.2021.0153
AStudyonIntegratingDEAandTOPSISBaseduponShephardDistanceFunction
LIYong jun,CAIHua quan
(SchoolofManagement,UniversityofScienceandTechnologyofChina,Hefei230026,China)
Abstract:TraditionalDEAmodelsmeasureDMU’sefficiencyusingasetofspecificweightswhichismostbene
ficialtothatDMU’sefficiencyscore.Somescholarspointoutthatitcan’tguaranteetheaccuracyofefficiencymeasure.TOPSISprovidesagoodideatosolvethisproblembecauseitconsidersbothidealandanti idealsitua tions.ByintroducingShepharddistancefunction,thispaperformallydevelopsageneralframeworkofintegrating
DEAandTOPSISforranking.Thepaperfirstdefinesideal
(anti ideal)DMUsandcorresponding(anti )pro ductionpossibilityset,thenbasedupontheideal(anti ideal)DMUs,elicitstheefficiencymeasurementmodelsandShepharddistancefunctionforeachDMU,andshowsacompleterankbyacomprehensiveindexcalledtherelativeclosenesstotheidealDMU.Attheend,thenumericalexperimenthasshownthatproposedmethodisrelativelyvalidandpracticable.ThispapercombinesadvantagesofDEAandTOPSIS,considersidealDMUsandanti idealDMUs,generatesfairerresults,andimproves2defectsofWL’smodel.Keywords:DEA;TOPSIS;rank;Shepharddistancefunction;productionpossibilityset;efficientfrontier
0 引言
数据包络分析(DEA)方法的一个主要作用就
是对决策制定单元(DMU)的排序[1~3]。但是,直接用传统D
EA模型(CCR[4]、BCC模型[5]
)中的多指标效率值进行评价存在诸多缺陷,主要有:(1)传统DEA模型只能分辨出决策单元是DEA有效还是非有效,不具备对决策单元进行分级、排序的
能力[1]
;(2)DMU之间排序取决于DMU的相对效
率值大小,传统DEA模型中计算效率值的权系数
只在对被评价单元最有利(使其效率值最大)的特定范围内取值,容易形成夸大长处、回避缺陷。因此,传统DEA模型中DMU效率值的这种乐观(
optimistic)确定方法很难保证效率评估结果的准确性[6]
由于存在以上的问题,学者们开始对DEA模型加以适当的改造和完善。Andersen和Petersen
提出的超效率模型[7]
,是解决传统模型缺陷(1)的
一个经典代表,该方法的基本思想是在评价某个决
策单元时,将其排除在决策单元的集合之外。这种方法在很多领域有着广泛的应用[8~10]。尽管如此,在某些情况下该模型可能无解以及有效DMU的超效率值可能不稳定[11~13]。类似的方法还有标杆(Benchmark)排序模型[14~16]。尽管如此,这些模型均没有涉及缺陷(2)的改进。
与之对应的是,Sexton等人提出的交叉效率评价方法[17],该方法主要是改变传统的自评体系(乐观方法)而采用自互评体系,从而达到纠正传统DEA模型缺陷(2)的目的,也有着广泛应用[18,19]。但是该方法的缺点也很明显,就是一般情况下交叉效率值不唯一[20]。从而由该模型得到的排序并不具有实际意义。Entani[21]等则首先从乐观(opti mistic)角度给出DMU的最大相对效率值,然后从悲观(pessimistic)角度计算了DMU的最小效率值,从而得到了DMU的一个效率区间。但是,显然该方法增加了在DMU排序中的难度。
此外,Wang和Luo(下文简记为WL)则[22]通过在生产可能集中引入虚拟的正负理想DMU,分别计算出DMU和正负理想DMU的最优或最劣的相对效率值,然后将DMU的最优(最劣)效率值与正(负)理想DMU的效率值差额分别取代TOP SIS[23]概念下的离正(负)理想解的距离,最后通过相对接近度这一综合指标得出DMU之间的排序。尽管如此,该模型仍明显存在以下缺点:(1)评价对象到正负理想解的距离数量级不同[24]。到正理想DMU的距离为小于1,而到负理想DMU的距离却大于1。(2)在评价的过程中,正负理想DMU的地位不同。主要表现在它们的效率值不等。(3)模型以效率差额作为评价对象到正负理想点的距离,混淆了效率差额和距离这两个概念。
Shi[25]等人在交叉效率模型的基础上引入了(负)理想前沿面,在确定权重时尽可能使DMU与(负)理想前沿面的偏差最小(大)化,使得权重的计算更中立和具有逻辑性。DEA与TOPSIS结合方法除上述外还有另一种方式。杜涛等[26]在评价医院效率时使用DEA方法计算出各个决策单元的松弛改进量,以松弛改进量构建TOPSIS方法的决策矩阵为决策单元排序;王钢和石奇[27]在评价商业银行经营效率时使用了TOPSIS DEA方法,同时借鉴了“垂足距离”的概念计算决策点到(负)理想决策点的距离;王春豪等[28]使用了非径向超效率DEA模型和TOPSIS效率测度模型评价中国区域公路物流的效率;Chitnis和Vaidya[29]先使用DEA方法出决策单元中的(负)理想决策单元,以此作为TOPSIS中的(负)理想点,计算印度银行
的效率。这类方法与本文的方法在思路上存在明显区别:在这类方法中,可以将DEA方法视为TOPSIS方法前的某种数据预处理步骤;在本文的方法中,DEA方法被融入到TOPSIS的方法中,成为TOPSIS方法中的一个步骤———计算当前DMU到(负)理想DMU的距离,成为其不可分割的一部分。
本文是在WL模型[22]的正负理想DMU基础上,通过引入Shephard距离函数、(反)生产可能集,给出了一种结合了DEA和TOPSIS新的有效排序方法。该方法能同时改进传统DEA两点缺陷,以及WL模型的三点缺陷。在后文中,本文将予以详细说明。
1 CCR效率指数与Shephard距离函数
设DMU
的投入向量x
∈Rm
和产出向量y
Rs
,j∈N={1,2,…,n}。那么,在评价DMU
的相对效率时,其最优CCR效率值可由以下模型求得:
   θ
=minθ
   s.t.∑n
j=1
λ
θx
∑n
j=1
λ
λ
0,j=1,2,…,n;θ,free(1)
其中,θ
是DMU
的最大CCR效率值。实际上,模型(1)可以简写为:
θ
=min{θ|(θx
,y
)∈T
CCR
}(2)
其中生产可能集T
CCR
是:
CCR
={(x,y)|∑n
j=1
λ
x,∑
j=1
λ
y,λ
0,j∈N}(3)Shephard距离函数这一重要的概念是最早由Shephard[30,31]提出,后由Banker等引入DEA研究
领域[5]。在生产可能集上T
CCR
,Shephard投入距离函数可定义为:
(x
,y
)=sup{γ∈R
|(x
/γ,y
)∈T
CCR
}(4)显然,有以下性质:
性质1 CCR模型的效率指数与Shephard的距离函数互为倒数,即:
(x
,y
)=1/θ
(5)
性质2 D
(x
,y
) 1成立的充分必要条件
是(x
,y
)∈T
CCR
在模型(1)中,若(x
)∈T,则显然有θ
(0,1],那么D
(x
,y
=1/θ
1,反之亦然。
第5期       李勇军,等:基于Shephard距离函数的DEATOPSIS方法研究
2 DEATOPSIS模型
本节将以Shephard距离函数为桥梁整合DEA和TOPSIS方法。为此,本文从定义理想DMU、负理想D
MU出发,在考虑(负)理想DMU的条件下分别给出DMU的(反)效率评价模型以及对应的Shephard距离函数,最后给出评价对象到理想DMU相对接近度的计算公式。其中,WL
[22]
首先
引入了虚拟的正(负)理想DMU的概念,本文在此基础上整理得到相关定理与性质,并引入Shephard距离函数。
2.1 离正理想DMU的Shephard距离函数
定义1
[22]
 正理想决策单元(IDMU)为DMUn+1,则其投入要素向量和产出要素向量分别为xn+1=(min1 j n
{xij}), i,yn+1(max1 j n
{yrj
}), r(6)那么,在考虑正理想DMU条件下,评价DMUd
的CCR相对效率模型为:
    θ+
=maxUT
yd
VT
xd
s.t.UTyj
VTxj
1,j=1,2,…,n+1
UT,VT 0
(7)
模型(7)中,U,V是权重向量,θ+
d是DMUd最优相对效率值。设模型(7)对应的最优解是U Td,V Td。
则有以下定理:
定理1 正理想决策单元DMUn+1是所有DMU的标杆(Benchmark),即DMUn+1在上述模型最优解
下的效率值恒有:θd
n+1
=U T
dyn+1
V Tdx
n+1=1, d∈N∪{
n+1}。证明 反证法:假设有θd
n+1<
1。根据上述模型的不等式,显然有:θdj
=U T
yj
V Tdx
U Tdyn+1V Tdxn+1
=θd
n+1<
1,j∈N\{d}。若令UTd=U Td/θdn+1,VTd=V T
d,
则(UTd,VTd)必然是模型(7)的可行解。是因为有:
    θdj
=UTdyjVTdxj<UTdyn+1VTdxn+1=U T
dy
n+1/θdn+1V T
dxn+1
=θdn+1/θ
n+1
=1,j∈N
而在可行解(UT
,VTd
)下,目标函数值为θd′=UT
dy
VT
dxd
=U T
dyd/θdn+1
V T
dx
=θ+
d/θdn+1。显然有:θd
′=θ+
d/θ+n+1
θ+d。则与θ+d是模型(7)的最优值(θ+
d θd
′)矛盾。因此假设不成立,又因d的任意性,
定理获证。依据C
C变换,模型(7)又可以转换为    maxuT
yd
    s.t.uTyj-vTxj
0,j=1,2,…,n+1(8)vT
xd=1uT,vT 0
模型(8)的对偶模型可简写为:
     θ+d=min{θ|(θxd,yd)∈T+
CCR}
(9)
其中T+CCR是包含正理想D
MU下的生产可能集,即:  T
CCR
={(x,y)|∑n+1
j=1
xjλj x,∑n+1
j=1
jλj y,λj 0,j=1,2,…,n+1}(10)那么,依据性质1,DMUd到正理想DMU(ID MU
)的Shephard距离函数是:D+d(xd,yd)=1/θ+
(11)
性质3 θ+n+1=1,θ+
d 1
,d∈N。2.2 离负理想DMU的Shephard距离函数
定义2[22] 负理想决策单元(ADMU)为DMU0
,则其投入要素向量和产出要素向量分别为xn+1=(max1 j n
{xij}), i,yn+1(min1 j n
{yrj
}), r(12)
定义3 DMUd的反效率值为θ-
=VT
xd
UT
yd
。从公式来看,反效率值的定义与DEA效率值
定义相反,投入要素的权重和为分子,而产出要素的权重和为分母。因此,投入要素量越大,反效率值越大,反之,产出要素量越大,反效率值越小。那么,在考虑ADMU条件下,评价DMUd的相对反效率值模型为:
     θ-
=maxVT
xd
UT
yd
     s.t.VTxj
UTyj
1,j=0.1,…,n
(13)
UT,VT 0
模型(13)中假设所有DMU的反效率值小于1,类似于传统的CCR模型中的效率值的评价。
定理2 在反效率评估中,负理想决策单元为DMU0是所有DMU的标杆(Benchmark)。(证明类似定理1)
上述模型可通过C-C转换后经取对偶模型可改为:
θ-d=min{θ|(xd,θyd)∈T-
CCR}
(14)
其中,考虑负理想点的反生产可能集T-
CCR为:
-CCR
={(x,y)|∑n
j=0
jλj x,∑n
j=0
yjλj
,λj
0,j=0,1,…,n}(15)811运 筹 与 管 理           2021年第30卷
同理,依据性质1,DMUd离ADMU的Shephard距离函数为:
D-d
(xd,yd
)=1/θ-d
(16)
性质4 θ-
=1,θ-d
1,d∈N
。根据性质3和4显然有:θ+0=θ-
0=
1,即正负理想D
MU的各自效率值或反效率值是相等的,这显然改进了WL模型的缺陷(2)。
此外,根据性质3和4有:θ+d 1,θ-d 1
, d∈N,所有DMU的效率值或反效率值均在0到1之间,因此,他们的数量级是相对的。显然,本文的模型改进了WL模型的缺陷(1)。
为了理解该反生产可能集与T+
CCR
,以及所构成的反效率前沿面与对应的效率前沿面的关系,我们以一个投入要素和一个产出要素为例,见图1
图1 正负理想DMU、(反)生产可能集以及
对应的(反)效率前沿面关系图
  在图1中,圆圈代表DMU,其中,理想DMU为DMUn+1,经过坐标原点O和DMUn+1的实心直线就是效率前沿面,而效率前沿面和X轴夹角的区域就是对应的生产可能集。相反,负理想DMU为DMU0,经过坐标原点O和DMU0虚线则是反效率前沿面,而该前沿面与Y轴交角的区域就是反生产可能集。
2.3 基于Shephard距离函数与正理想DMU的
相对接近度
    根据TOPSIS定义,DMUd到理想DMU的相对接近度为
RCd=D-
d(xd,yd
)D-d(xd,yd)+D+d(xd,yd)=θ+
θ+d+θ-
(17)
3 算例分析与实证
我们采用的样本数据来自于文献[7],见表1。表1中给出了5个DMU,每个DMU有两种投入要素,一种产出要素。为了便于比较,由本文方法所得的排序结果以及WL、Wu模型的结果分别放在
表2的二、三、四三列。
表1 样本数据蔡国权简介
DMUInput1Input2Output1
CCRSuperCCR
1212111228111.3158355111.24104111.25510610.75
0.75
IDMU241ADMU
10
12
  首先,从IDMU以及ADMU的效率值比较来看,只有本文给定的方法下两个理想DMU的效率值是相等的,而WL、Wu模型下的结果均不相等。因此,只有本文给出的正负理想DMU在效率评估中的地位是对等的。其次,在对5个决策单元评价中,Wu模型下就有3个DMU(DMU1、DMU4和DMU5)的RC值均是0,因此,Wu模型的排序能力
值得商榷。而在WL模型中,θ-d 1而θ+
d 1
,d∈N
,显然这两种效率值的数量级是不同的。Wu在文中也注意到了这个缺陷,相比之下,本文方法下的两
种效率值均是大于0而小于1的。这也正好验证了性质3和4的正确性。此外,注意文献[24]中指出DMU1和DMU4既在传统的效率前沿面上又在反效率前沿面上,同样在本文中,考虑IDMU以及ADMU下的效率前沿面、反效率前沿面上也均有DMU1和DMU4,见图2。因此,DMU1和DMU4是同等“有效”的,它们的排名也应该是等同的。正如本文的结果所示为RC1=RC4=
0.5。9
11第5期       李勇军,等:基于Shephard距离函数的DEATOPSIS方法研究
表2 三种不同模型下的排序结果
DMUOurmethod
WL
Wu
θ+
θ-
RCdθ+
θ-
RCdθ+
dθ-
dRCd1110.50.71429
10.244180.71429
0.20210.66670.611.42110.5222410.30.1304330.80.50.6153811.54290.560610.40.23077
4110.50.714291.17390.335850.714290.2050.66671
0.4
0.6251
0.22803
0.6250.2
IDMU1
1.667
1.667
ADMU
10.6920.2
  在表3给出了本文方法、WL模型[22]
、Wu模型[
24]以及超效率模型[7]
下基于表1数据的排序相关分析。结果显示与超效率模型的相关性最强的是本文方法(相关系数值为0.84147)
图2 基于正负理想DMU下的(反)效率前沿面
表3 四种模型下排序结果的相关分析
Ourmethod
WLWuSuperCCROurmethod
10.92885
0.851580.84147WL/10.93342
0.76317Wu//10.51812
SuperCCR
4 结论与展望
本文以Shephard距离函数为桥梁构建了一个基于DEA和TOPSIS新的排序方法。该方法改进了传统DEA理论的两个重要缺陷。此外,该方法既结合了DEA在效率评价中不需要计算要素权重系数的优势,而又结合了TOPSIS方法在排序过程中的多角度(乐观和悲观)考察对象的优点,从而在一定程度上保证了排序结果的公正性。最后,通过具体的算例说明了该方法的有效性。因此,该方法有望为决策者提供较为有效的决策支持。
本文给出的DEATOPSIS排序方法是以CCR模型下的生产可能集为例,该方法也可轻易推广到
BCC模型下的生产可能集,即在公式(3)中加入
∑n
j=1
λj
=1约束条件。由于篇幅所限,
作者就不在此详述,感兴趣的读者不妨亲自尝试。参考文献:
[1]AdlerN,FriedmanL,Sinuany SternZ.Reviewof
rankingmethodsinthedataenvelopmentanalysiscontext[J].Europeanjournalofoperationalresearch,2002,140(2):249 265.[2]JahanshahlooGR,JuniorHV,LotfiFH,etal.Anew
DEArankingsystembasedonchangingthereferenceset[J].EuropeanJournalofOperationalResearch,2007,181(1):331 337.[3]EmrouznejadA,YangGL.Asurveyandanalysisofthe
first40yearsofscholarlyliteratureinDEA:1978~2016[J].Socio EconomicPlanningSciences,2018,61:4 8.[4]CharnesA,CooperWW,etal.Measuringtheefficiency
ofdecisionmakingunits[J].Europeanjournalofopera
tionalresearch
,1978,2(6):429 444.[5]BankerRD,CharnesA,CooperWW.Somemodelsfor
estimatingtechnicalandscaleinefficienciesindataenvel
opmentanalysis
[J].Managementscience,1984,30(9):1078 1092.[6]ChambersRG,F reR,JaenickeE,etal.Usingdomi
nanceinformingboundsonDEAmodels:thecaseof
experimentalagriculturaldata
[J].JournalofEconomet rics
,1998,85(1):189 203.[7]AndersenP,PetersenNC.Aprocedureforrankingeffi
cientunitsindataenvelopmentanalysis[J].Management
science
,1993,39(10):1261 1264.[8]NahraTA,MendezD,AlexanderJA.Employingsuper
efficiencyanalysisasanalternativetoDEA:anapplica
tioninoutpatientsubstanceabusetreatment
[J].Europe anJournalofOperationalResearch,2009,196(3):1097 1106.[9]AvkiranNK.AssociationofDEAsuper efficiencyesti
mateswithfinancialratios
:investigatingthecaseforChinesebanks[J].Omega,2011,39(3):323 334.[10]DuJ,WangJ,ChenY,ChouSY,etal.Incorporating
healthoutcomesinpennsylvaniahospitalefficiency:an
021运 筹 与 管 理           2021年第30卷

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。