数据包络分析(DEA)方法

二、数据包络分析(DEA)方法数据包络分析(dataenvelopmentanalysis,DEA)是由著名运筹学家Charnes,Cooper和Rhodes于1978年提出的，它以相对效率概念为基础，以凸分析和线性规划为工具，计算比较具有相同类型的决策单元(Decisionmakingunit ，DMU)之间的相对效率，依此对评价对象做出评价 。DEA方法一出现，就以其独特的优势而受到众多学者的青睐， 现已被应用于各个领域的绩效评价中 [2],[3] 。在介绍 DEA方法的原理之前，先介绍几个基本概念 :决策单元一个经济系统或一个生产过程都可以看成是一个单位 (或一个部门)在一定可能范围内，通过投入一定数量的生产要素并产出一定数量的 “产品 ”的活动。虽然这种活动的具体内容各不相同，但其目的都是尽可能地使这一活动取得最大的 “效益”。由于从“投入”到“产出”需要经过一系列决策才能实现，或者说，由于出”是决策的结果，所以这样的单位 (或部门)被称为决策单元 (DMU)。因此，可以认为，每个 DMU(第DMU 常记作DMUi)都表现出一定的经济意义，它的基本特点是具有一定的投入和产出，并且将投入转化

“产i个成产出的过程中，努力实现自身的决策目标。在许多情况下，我们对多个同类型的 DMU 更感兴趣。所谓同类型的 DMU，是指具有以下三个特征的DMU 集合：具有相同的目标和任务；具有相同的外部环境；具有相同的投入和产出指标。生产可能集设某个DMU在一项经济(生产)活动中有m项投入，写成向量形式为x(x1,,xm)T；产出有s项，写成向量形式为y(y1,,ys)T。于是我们可以用(x,y)来表示这个DMU的整个生产活动。定义1.称集合T {(x,y)|产出y能用投入 x生产出来}为所有可能的生产活动构成的生产可能集。在使用 DEA方法时，一般假设生产可能集 T满足下面四条公理 :公理1(平凡公理):(xj, yj)T,j 1,2, ,n。n公理2(凸性公理):集合T为凸集。如果(xj,yj)T,j1,2,,n,且存在j0满足j1j1nn则(j1jxj,j1jyj)T。公理3(无效性公理)：若x,y????,则(x,y)T。，公理4(锥性公理):集合T为锥。如果 x,y T那么 (kx,ky) T对任意的k 0。若生产可能集Ｔ是所有满足公理 1,2,3和4的最小者，则 T有如下的唯一表示形式nn。Tx,y|xjjx,yjjy,j0,j1,2,,nj1j1技术有效与规模收益'

'(1)技术有效： 对于任意的

(x,y)

T，若不存在

y

y，且(x,y)T

，则称

(x,y)

T为技术有效的生产活动。(2)规模收益：将产出和投入的同期相对变化比值

k

y/

x称为规模效益。k若1，说明规模收益y

x递增，这时可以考虑增大投入；若

k 1

，说明规模收益递减，这时可以考虑减小投入；若

k 1

，说明规模收益不变，且称为规模有效。（一） DEA方法原理与 CCR模型DEA 方法的基本原理是：设有 n个决策单元 DMUj(j 1,2, ,n)，它们的投入，产出向量分别为：Xj (x1j,x2j, ,xmj)T 0,，Yj (y1j,y2j, ,ysj)T 0, j 1, ,n。由于在生产过程中各种投入和产出的地位与作用各不相同，因此，要对 DMU 进行评价，必须对它的投入和产出进行 “综合 ”，即把它们看作只有一个投入总体和一个产出总体的生产过程，这样就需要赋予每个投入和产出恰当的权重。假设投入、产出的权向量分别为v(v1,v2,,vm)T和u(u1,u2,,us)T，从而就可以获得如下的定义。suyTrrj定义2.称uYjr1n)为第j个决策单元DMUj的效率评价指数。jTm,(j1,2,vXjvixiji1根据定义可知，我们总可以选取适当的权向量使得 j1。如果想了解某个决策单元，假设为DMUo(o{1,2,,n})在这n个决策单元中相对是不是“最优”的，可以考察当u和v尽可能地变化时，o的最大值究竟为多少?为了测得o的值，Charnes等人于1978年提出了如下的CCR(三位作者名字首字母缩写)模型：sr1uryroMaximizemovixioi1suryrjsubjecttor11,j1,2,,n,(1)mvixiji1ur0,vi0,r,i.利用Charnes和Cooper(1962)[4]mvixio,提出的分式规划的Charnes-Cooper变换:t1/i1rtur,(r1,,s),itvi,(i1,,m)变换后我们可以得到如下的线性规划模型：sMaximize ryro o,r1msubjectto ixio 1,i1(2)smryrjixij0,j1,,n,r1i1r,i0,r1,,s;i1,,m.根据线性规划的相关基本理论，可知模型 (2)的对偶问题表达形式：Minimizeonxsubjectto ijj1

j o ,iox1,2i, , ,mn (3)yrj j y,ror1,2, s,,j 10,j1,2,n,.上述的模型是基于所有决策单元中 “最优 ”的决策单元作为参照对象，从而求得的相对效率都是小于等于1的。模型(2)或者(3)将被求解 n次，每次即得一个决策单元的相对效率。模型 (3)的经济含义是：为了评价 ( {1,2, ,})DMUoo n 的绩效， 可以用一组假想的组合决策单元与其进行比较。 模型(3)的第一和第二个约束条件的右端项分别是这个组合决策单元的投入和产出。从而，模型 (3)意味着，如果所求出的效率最优值小于 1，则表明可以找到这样一个假想的决策单元，它可以用少于被评价决策单元的投入来获取不少于该单元的产出，即表明被评价的决策单元为非 DEA有效。而当效率值为 1时，决策单元为 DEA有效。有关DEA有效根据松弛变量是否都为零还可以进一步分为弱 DEA有效与 DEA有效两类。即通过考察如下模型中的 si(i 1, m)与sr(r 1, ,s)的值来判别。m sMinimizeo(sisr)i1r1nsubjecttoxijjsioxio,i1,,mj1(4)nyrjjsryro,r1,,sj1j,si,sr0,i,j,r.其中 为非阿基米德无穷小量。根据上述模型给出被评价决策单元定义3.若模型(4)的最优解满足定义4.若模型(4)的最优解满足定义5.若模型(4)的最优解满足

DMUo(o{1,2, ,n}) 有效性的定义：*1，则称DMUo为弱DEA有效。*1，且有si0，sr0成立，则称DMUo为DEA有效。o*1，则称DMUo为非DEA有效。o对于非 DEA有效的决策单元，有三种方式可以将决策单元改进为有效决策单元：保持产出不变，减少投入；保持投入不变增大产出；减小投入的同时也增大产出。 CCR模型容许 DMU 在减小投入的同时也增加产出。对于 CCR模型，可以通过如下投影的方式将其投向效率前沿面，从而投影所得的点投入产出组合即为DEA有效。?*x*x**xsiio(1o)xiosixio,i1,,miooio?y*ysryro,r1,,s.roro上述投影所得值与原始投入产出值之间的差异即为被评价决策单元欲达到有效应改善的数值， 设投入的变化量为 xio，产出的变化量为 yro：xx?xx*xsi*ioioio(oio),i1,,m?ioyyys*yrs(ro),1,,.rorororro（二） BCC模型CCR 模型是假设生产过程属于属于固定规模收益，即当投入量以等比例增加时，产出量应以等比增加。然而实际的生产过程亦可能属于规模报酬递增或者规模报酬递减的状态。为了分析决策单元的规模报酬变化情况，Banker,Charnes与Cooper以生产可能集的四个公理以及Shepard距离函数为基础在1984年提出了一个可变规模收益的模型，后来被称为BCC的模型[5]。线性形式的BCC模型可表示为：sMaximizerryrouo,1msubjecttoiixio1,(5)1smr1ryrji1ixijuo0,j1,,n,r,i0,r1,,s;i1,,m.含松弛变量形式的BCC对偶模型msMaximizeo(sisr)i1r1nsubjecttoxijjsioxio,i1,,mj1nyrj j sr yro,r 1, ,s (6)1n1j1j,si,sr0,i,j,r其中 为非阿基米德无穷小量。 根据 BCC模型中的 uo的取值大小， Banker 和Thrall(1992) [6]提出如下判别方法来判断模型(5)的规模收益。定理1[6].假设含有投入产出组合(xo,yo)的DMUo是有效的，那么下面的条件可以判别模型(1)之下DMUo的规模收益：(i)对于投入产出组合(xo,yo)规模收益不变当且仅当在某个最优解情况下有uo*0；(ii)对于投入产出组合(xo,yo)规模收益递增当且仅当在所有最优解情况下都有uo*0；(iii)对于投入产出组合(xo,yo)规模收益递减当且仅当在所有最优解情况下都有uo*0。其中uo*代表模型(5)中的最优解。该定理的证明参见文献[6]。CCR模型或者 BCC模型计算出来的效率可能存在多个效率值为 1的情形，为了进一步区分这些有效决策单元，常用的方法有超效率模型，交叉效率模型以及双前沿数据包络分析模型。下面依次做个简单介绍。（三） 超效率模型CCR 模型在计算效率值时，经常会出现多个有效的决策单元 (效率值为 1)的情形，从而使得有效决策单元之间无法进行比较分析。 Andersen 和Petersen(1993) [7]为了实现决策单元的完全排序， 将被评价的决策单元从效率边界中剔除，以剩余的决策单元为基础，形成新的效率边界，计算剔除的决策单元到新的效率边界的距离。由于剔除的决策单元不被效率边界所包围，对于有效的决策单元而言，其计算出来的新效率值就会大于 1，而对于无效的决策单元而言，其所得的效率值不变，仍小于 1，从而使得全体决策单元可以实现完全排序。由于有效的决策单元效率大于 1，从而就有了超效率 (Super-efficiency) 的概念。基于CCR模型的超效率DEA模型为：Minimizensubjecttoxijjxio,i1,2,,m,1jon(7)yrjjyro,r1,2,,s,1oj0,jo.Banker和Chang(2006)[8]证实了超效率极易受离群值的影响，因此该方法可以用来检测数据集中是否存在离群值。（四）交叉效率模型为了解决 DEA 有效决策单元的排序和比较问题， Sexton等人(1986)[9]提出了交叉效率评价的概念。所谓交叉效率评价就是每个 DMU 分别确定一组输入输出权重，供所有的 DMUs 评价使用，其中：用 DMU自身确定的权重评价自己的效率，称为自我评价效率；用其它 DMU 确定的权重评价自己的效率，称为交叉效率或同行评价效率。以 表5—1为例，交叉效率评价的实质是对每个 DMU 同时进行自评和同行评价，这样不仅考虑 DMU 自评的最好相对效率，而且还考虑了 DMU 同行评价给出的交叉效率，利用自我评价和交叉效率的平均值作为衡量 DMU 绩效的综合指标， 该指标不仅较好地解决了 DMUs 间排序和比较问题，而且解决了 CCR模型由于输入输出权重不一致性导致的不可比较问题。Sexton等人(1986)通过引入二级目标来确定输入输出权重、消除权重的不唯一性。随后 Doyle 和Green(1994,1995) [10],[11]从同行评价的角度解释了交叉效率的含义，并给出了后来的到广泛引用的二级目标函数-攻击型计算方式和仁慈型计算方式， 下面两个模型依次为攻击型交叉效率模型和仁慈型交叉效率模型：表5—1交叉效率示意表交叉效率决策单元算术平均值12?n1?1n11121nj11jn2?1n21222nj12jnn?1nnjn1n2nnj1n攻击型交叉效率模型：snMinimizeurkyrjr1j1,jkmnSubjecttovikxij1,i1uyj1,jkvx0,(8)srkrkmiik*r1i1smurkyrjvikxij0,j1,,n;jk,r1i1u0,r1,,s,rkvik0,i1,,m.仁慈型交叉效率模型:snMaximizeurkyrjr1j1,jkmnSubjecttovikxij1,i1j1,jkvxsuym0,(9)rkrkiik*r1i1smurkyrjvikxij0,j1,,n;jk,r1i1urk0,r1,,s,vik0,i1,,m.然而，至今仍无一个准则来判别什么情况下使用攻击型或者是仁慈型。为了避免目标函数选择上的两难，Wang和Chin(2010a)[12]提出了一种中性交叉效率模型。其模型形式如下所示：uyMaximizeMinimizeroro,mr{1,,s}i1vioxios1uroyrosubjectto*r,oom1vioxiosyru1rorj(10)jom1,j1,,n;jo,1vioxijuro0,r1,,s,vio0,i1,,m.利用Charnes-Cooper 的变换公式，可得中性交叉效率模型的线性模型Maximizemsubjecttovioxio1,i1suy*,rorooor1srom(11)rjioijuyvx0,j1,,n;jo,r1i1uy0,r1,,s,rorovio0,i1,,m.0.交叉效率模型还有其他一些改进方式，例如：Liang等人(2008a)13年提出了3个可供选择的二级目标计算方式；Liang等人(2008b)[14]将非合作博弈理论与交叉效率评价方法结合起来，提出了博弈交叉效率的概念，并设计了算法求解博弈交叉效率值，同时证明了该博弈交叉效率值即为纳什均衡点；Wang和Chin(2010b)[15]提出了一些可选择性交叉效率评价模型；Wang和Chin(2011)[16]在交叉效率的研究中率先引入有序加权平均算子 (Orderedweightedaveragingoperator,OWA) ，很好的体现了决策者的各种偏好，尤其是对不合理的交叉效率评价值赋予较小的权重，从而使得最终的评价结果更为科学合理。有兴趣的读者可以进一步参阅其他有关交叉效率的相关论文。（五）几何平均效率模型为了区分有效决策单元的排序难问题， Wang等人(2007)[17]于2007提出了悲观效率模型，并将其与乐观效率模型相结合，提出了基于几何平均值的双前面数据包络分析方法。基于悲观前沿面的数据包络分析模型为：sMinimizer1ryro,msubjecttoixio1,(12)i1smryrjixij0,j1,2,,n,r1i1r,i0,r1,2,,s;i1,2,,m,其中r和vi是非负权重。模型(12)与模型(2)的区别在于：模型(12)计算所得效率均大于等于1，而模型(2)所得的效率值均小于等于1。基于几何平均值的双前沿数据包络分析方法就是将模型(12)所得的效率与模型(2)所得的效率通过几何平均的方式加以综合，即：***ooo其中*为综合后的DMU(o{1,2,,n})的效率值，而*和o*分别对应该决策单元在模型(2)与模型(12)下ooo的最优效率值。下图为有效前沿面和无效前沿面的一个演示图。产出2/投入8AE6 CG 有效前沿面4B D无效前沿面F H I1 J1 2 4 6 8 产出1/投入图5-1 决策单元的有效和无效前沿面（六）最优决策单元的选择在实际应用中，决策者有时候关心的是哪个方案或者哪个决策单元是最优的，而对于其他单元的排序并不在意。因此，如何利用 DEA模型直接寻求最优决策单元成为学者们所感兴趣的问题。 Amin和Toloo(2007)[18]提出了一个混合整数线性规划模型，采用两步法以期实现寻求最优决策单元。然而随后 Amin(2009)[19]发现这种两步法有时会产生两个或者两个以上的最优决策单元，因此他提出一个非线性混合整数模型。Foroughi(2011) [20]发现Amin的非线性规划模型在有些情况下是不可行的。 不过 Foroughi(2011) 的模型存在着一些冗余的约束且对输入输出权重给定了保证域，并且该模型易受离群值 (outliers)的影响，从而导致所选择的最优决策单元不正确。因此， Wang和Jiang(2012) [21]提出了三种混合整数线性规划模型来改进Foroughi(2011) 的模型中所存在的问题。这三种最优决策单元选择的模型分别为：基于不变规模收益的混合整数线性规划模型的最优决策单元选择方法m n x

s n yMinimizeviijurrji1j1r1j1smSubjecttouryrjvixijIj,j1,,n,r1i1nIj1,(13)j 1Ij {0,1}, j 1, ,n,ur1r1,,s,,(ms)max{yrj}j1,i1,,m,vi(ms)max{xij}j其中Ij(j1,,n)是二元变量,且只有一个变量可以取非零值1。如果Io1,那么约束条件smsmuryrjvixijIj对应的DMUo的约束为r1uryroi1vxiio1,即允许DMUo的效率值大于1，而r1i1smsm其余的DMU的约束uryrjvixijIj与原始的CCR模型的约束相同,也就是uryrjvixij0对于r1i1r1i1任意的j{1,,n}除了jo。因此,只有最有效的决策单元的效率值会大于1，而其余决策单元的效率均小于等于1。权重约束沿用(Sueyoshi,1999[22])提出的松弛变量模型中的形式，该约束形式在实际应用中被y广泛采用，即ur(1/((ms)max{rj}))(r1,,s);vi(1/((ms)max{xij}))对于任意的(i1,,m).jj基于投入导向的BCC模型的混合整数线性规划最优决策单元选择方法模型(13)是基于不变规模收益下的最优决策单元的选择方法。该方法可以拓展到可变规模收益的情形如下所示，该模型的形式是基于投入导向的 BCC模型下的形式 ：mnxnv0snyMinimizeviijurrji1j1r1j1s mSubjectto uryrj vixij v0 Ij, j 1, ,n,r1 i 1nIj1,(14)j 1Ij {0,1}, j 1, ,n,1ur,r1,,s,(ms)max{yrj}j1vi ,i 1, ,m,(m s)max{xij}jv0 无符号限制 .基于产出导向的BCC模型的混合整数线性规划最优决策单元选择方法同理可得，基于产出导向的可变规模收益的 BCC形式下的混合整数线性规划模型如下：mnxsnynu0Minimizeviijurrji1j1r1j1suyu0mxIj,j1,,n,Subjecttorrjviijr1i1suyrjrnIj1,(15)j 1Ij {0,1},j 1, ,n,ur1,r1,,s,s)max{yrj}(mj1,i1,,m,vi(ms)max{xij}ju0 isfreeinsign,suy其中约束条件rrju00(j1,,n)是为了保证全体产出是非负的，因为负的产出没有意义。r1（七）举例说明下面用 3个例子来说明 DEA方法的应用。例1：假设现有七个被评价的决策单元，投入、产出项各有一项，投入项为 X，产出项为 Y，输入如下表所示。此时七个决策单元的相对位置如图 5—2所示。在CCR模型下，连接原点与点 B的射线构成前沿面，如图中所示，其余的点均位于该前沿面的下方。表5—2七个决策单元的投入、产出数据DMUXYEfficiencyA210.5000B331.0000C860.7500D620.3333E540.8000F1060.6000G74.50.6429产出876 C F5G4 E3 B2D1 A1 2 3 4 5 6 7 8 9 10 11 投入图5—2七个决策单元的分布及其在生产前沿面上的投影从图2中可以看出， 只有决策单元 B位于生产前沿面上， 而其他所有决策单元均位于该生产前沿面的下方，即A,C,D,E,F,G 均为非DEA有效，从表 5—2最后一列的效率值大小也很容易得到确认。 为了使非 DEA有效决策单元为 DEA有效，可依图中箭头所示的方向将非 DEA有效的决策单元往前沿面上投影。 A,C,D,F,G均为减小投入而保持产出不变； 而E给出了三种投影方式 (减小投入产出不变； 保持投入不变增大产出；或者同时减小投入和增大产出 )。例2:五个先进制造技术的甄别 ,数据来源于 Wang和Chin(2009) [23]。表5—3 五个先进制造技术的数据及其乐观、悲观以及几何平均值投入产出决策单元X1X2Y乐观效率值悲观效率值几何平均值A4072101.00001.60001.2469B32121050.56391.00000.7509C52203041.00001.73711.3180D35132000.98381.74151.3089E3281500.85801.42861.1071对于每一个决策单元而言，可通过求解模型 (2)和(12)获得全体 DMUs的乐观和悲观效率，结果如上表所示。下面简单介绍一下求解过程和技术实现。以第一个决策单元的 CCR效率(即乐观效率 )为例，将数据代入模型(2)即得模型(16)，显然这是个较为复杂的线性规划模型，需要借助软件计算才会更为简便。因此本书分别给出了 Lingo以及Matlab下的CCR模型的编程。 Lingo的编程一次也只能计算一个 (见下面程序后的计算说明 )，而Maltab程序相对而言更为简便，其可以很快地计算出所有决策单元的效率。此例中通过软件计算所得，在乐观效率下，所得效率为 表 5—3的第五列所示，全体单元的优序关系为：C=A>D>E>B 。，决策单元 A与C均为DEA有效，而 B,D,E为非DEA有效。在悲观模型下，所得的效率值为表5—3的第六列所示，决策单元 B为DEA 无效，而其他单元均为非 DEA 无效，其优序顺序为：D>C>A>E>B 。由此可见，在乐观前沿面和在悲观前沿面下的排序存在着一定的差异。

表5—3的最后一列的值为乐观和悲观效率的几何平均值，显然 Wang等人(2007)提出的该几何平均值较好的综合了乐观和悲观前面的两部分信息，从而五个单元合理的排序为： C>D>A>E>B。Maximize 1 21011+72=1,2101(401+72)0,1051(321+122)0,(16)subjectto304200

11

(52 1+20 2) 0,(35 1+132)0,1501 (321+8 2) 0,1 0, 1 0, 2 0.下面给出 LINGO与Matlab的程序：例2的LINGO程序实现： (以计算第一个决策单元的乐观效率为例 )MODEL:sets:DMU/1..5/:S,T,P; !Decisionmakingunits;II/1..2/:w; !inputindex;OI/1/:u; !outputindex;IV(II,DMU):X; !inputvariable;OV(OI,DMU):Y; !outputvariable;endsetsdata:P=10000;X=403252353271220138;Y=210105304200150;enddatamax=@sum(DMU:P*T);@for(DMU(j):S(j)=@sum(II(i):w(i)*X(i,j));T(j)=@sum(OI(i):u(i)*Y(i,j));S(j)>=T(j));@sum(DMU:P*S)=1;END可得5个决策单P的值(10000)分别替换为(01000),(00100),(00010),(00001),在上述程序中，元的最优效率值依次为1.0000，0.5639，1.0000，0.9838，0.8580。例2的Matlab程序实现：clearall;X=[403252353271220138];Y=[210105304200150];n=size(X',1);m=size(X,1);s=size(Y,1);A=[-X'Y'];b=zeros(n,1);LB=zeros(m+s,1);UB=[];fori=1:n;F=[zeros(1,m)-Y(:,i)'];Aeq=[X(:,i)'zeros(1,s)];beq=1;w(:,i)=linprog(F,A,b,Aeq,beq,LB,UB);E(i,i)=Y(:,i)'*w(m+1:m+s,i);endwEomega=w(1:m,:)mu=w(m+1:m+s,:)EE=diag(E)运行上述 Matlab 程序，即可得全体 DMUs 的CCR效率值。例3.现有14家国际航空公司，数据来源于Tofallis(1997)[24]。已知投入有三项，产出有两项，分别为：x1:飞机容量吨公里x2:营业费用x3:其他资产（预定系统，便利性以及流动资产）y1:每公里乘客数y2:非客运收益表5—4 14家航空公司的数据投入产出DMUx1x2x3y1y215723323920032667769725895422545573081539324099956062671240551266413565749932136473415635518318807832360451361908080323272950115727460334572360221129698120976779647452363200196587334135812650412971056541878191619277972111255980983310419253398125728248122542775498213471517922485313325431422793987441451225281404表5—5CCR效率及其非有效决策单元的改进DMUjCCR效率.投入产出x1jx2jx3jy1jy2j10.8684-753-916-2640020.3379-3903-2940-40323569030.9475-1265-502-329073140.9581-569-1235-13500510000060.9766-447-402-770810710000080.8588-1709-957-20080090.9477-344-175-1392001010000011100000121000001310000014100000利用CCR模型以及将非有效DMU改进为有效DMU的投影公式，可得表5—5的结果。从表中可知，决策单元 5,7,10,11，12,13,14 为DEA 有效，而其它单元为非 DEA有效。对于非有效决策单元，例如对第一家航空公司而言，它的第一项投入应减少 753，第二项投入应减少 916，第三项投入应减少 264，同时保持产出不变，这时该航空公司可达 DEA 有效。DMU4，DMU8和DMU9与DMU1类似也均需减少该三项投入。而对 DMU2而言，其前三项投入应分别减少 3903，2940和 4032，第一项产出需增加 3569，第二项产出保持不变可达有效。而 DMU3和DMU6在减少三项投入的同时，还需要增加第二项产出才会有效。利用攻击型交叉效率模型，我们可得如下表 (表5—6)所示的 14家航空公司的交叉效率表以及其排序。从表中可以看出第 5家航空公司的相对效率为 0.7983，为所有航空公司中最优，其次是第 11家航空公司，其交叉效率值为 0.7742。而第 2家航空公司的交叉效率值为 0.1652，为 14家航空公司中最差。利用仁慈型交叉效率模型，我们可得如下表 (表5—7)所示的 14家航空公司的交叉效率表以及其排序。从表中可以看出第 11家航空公司的相对效率为 0.9193，为所有航空公司中最优， 其次是第 13家航空公司，其交叉效率值为 0.9190。而第 2家航空公司的交叉效率值为 0.1894，为 14家航空公司中最差。此结果与攻击型交叉效率模型所得的结果又较大的差异，然而至今仍无一个准则可以清晰的告诉决策者何时该选择攻击型模型或者是仁慈型模型。因此均对不同的决策问题，选择的模型的不同，所得结果可能出入较大。为此，学者们提出了一些改进模型， 例如Wang 和Chin(2010a) 的中性交叉效率模型， 以及Liang等人(2008a)的博弈交叉效率模型都可以较好的避免这个问题。表5—6攻击型交叉效率值表目标DMU平均排DMU1234567891011121314交叉序效率10.86840.45010.62250.86840.44180.47260.76790.78810.70310.41580.33900.70430.47110.47260.59901220.17190.33790.04720.17190.02240.02470.27700.27240.28080.24650.11520.27890.04170.02470.16521430.88260.19420.94750.88260.65660.68980.64680.68330.62250.25590.19680.62610.74220.68980.62261140.95810.42590.70340.95810.66830.69730.76290.78500.69910.40270.47390.70160.49370.69730.6734750.96530.36581.00000.96531.00001.00000.70110.73590.77780.52720.63820.78190.71811.00000.7983160.88180.11080.95630.88180.96320.97660.57450.60840.50990.13760.17030.51410.67660.97660.6385970.92110.77810.47730.92110.31080.33821.00001.00000.83950.54160.40000.83830.36580.33820.6478880.78130.61140.51620.78130.26830.29240.84150.85880.82080.57030.30110.81940.44180.29240.58551390.78550.72780.50750.78550.24550.26770.88810.90720.94770.75010.35280.94520.45370.26770.630910100.78210.63540.65200.78210.33370.35640.76500.79441.00001.00000.49421.00000.58710.35640.68136111.00001.00000.42871.00000.42020.44181.00001.00001.00000.81071.00001.00000.29610.44180.77422120.94620.63360.75000.94620.40850.43950.90820.93950.99980.76470.42441.00000.63980.43950.73145131.00000.42561.00001.00000.41830.45550.95111.00001.00000.58550.21291.00001.00000.45550.75033141.00000.22771.00001.00000.98061.00000.69190.72750.64780.27470.32990.65210.70971.00000.73164表5—7 仁慈型交叉效率值表目标DMU平均排DMU交叉序1234567891011121314效率10.86840.45010.62250.86840.84920.47260.81080.78810.70310.75120.86840.77130.86840.86840.75431220.17190.33790.04720.17190.17350.02470.24790.27240.28080.20580.17190.20250.17190.17190.18941430.88260.19420.94750.88260.88440.68980.72320.68330.62250.78460.88260.80720.88260.88260.7678940.95810.42590.70340.95810.94130.69730.82280.78500.69910.81120.95810.83410.95810.95810.8222650.96530.36581.00000.96531.00001.00000.77040.73590.77781.00000.96531.00000.96530.96530.8912360.88180.11080.95630.88180.87800.97660.66150.60840.50990.71760.88180.74780.88180.88180.75541170.92110.77810.47730.92110.87950.33821.00001.00000.83950.78080.92110.80120.92110.92110.8214780.78130.61140.51620.78130.77020.29240.84580.85880.82080.75320.78130.76310.78130.78130.72421390.78550.72780.50750.78550.78890.26770.87820.90720.94770.83750.78550.83690.78550.78550.759010100.78210.63540.65200.78210.82500.35640.77800.79441.00001.00000.78210.97190.78210.78210.78038111.00001.00000.42871.00001.00000.44181.00001.00001.00001.00001.00001.00001.00001.00000.91931120.94620.63360.75000.94620.96020.43950.93620.93950.99981.00000.94621.00000.94620.94620.88504131.00000.42561.00001.00001.00000.45551.00001.00001.00000.98431.00001.00001.00001.00000.91902141.00000.22771.00001.00001.00001.00000.77950.72750.64780.85691.00000.88381.00001.00000.86595应用数据包络分析需注意的事项：DEA作为一种非参数方法，将数学、经济和管理的概念与方法相结合，是处理多目标决策问题，解决在经济和管理中评价具有多个投入、多个产出问题的有力工具。主要适合于同种类型部门或单位间的2.

相对有效性排序和评价， 可以通过在生产前沿面上的投影分析， 发现非 DEA有效的产生原因以及改进方向，调整资源投入量和效益产出量使决策单元达到 DEA有效。应用DEA模型进行评价， 不必事先确定指标权重， 只需假定由决策单元的投入产出指标组成的状态可能集满足凸性、无效性、锥性以及最小性等条件即可。 DEA方法本身包含指标的权重分配过程，在计算不同决策单元的最大有效性时，指标的权重是动态可变的，最后排序的结果是每个决策单元在最有3.

利于自身的权重下的结果。应用DEA方法进行评价分析时， 一般要求决策单元数目应大于投入产出变量数目之和。 根据经验法则，最好是决策单元数目应大于投入产出变量数目之和的两倍，这样计算所得的结果才能具有较好的区分性，否则容易出现多个决策单元有效而无法进一步区分的情形。此时，可以采用超效率模型或者是交叉效率模型进行相对效率分析。也可以将这些有效决策单元再采用其他综合评价方法进行分析。4.投入产出指标的确定，一般是根据资源投入量与效益产出量确定。 DEA模型求解时，一般要求投入产出指标具有非负性。如果遇到负的投入指标，一些学者认为可以考虑将取绝对值后纳入产出指标进行考虑，不过这种方法的合理性以及此方面的研究还尚未取得一致的认识。在实践中，通常有两种导向的模型可以供决策者选择，一类是投入导向模型（即在相同产出水平下，比较投入资源的使用情况） ，一类是产出导向模型， 人们通常只从投入导向或者产出导向的角度去分析决策单元的相对有效性， 不过这两种角度在很多时候计算所得的结果是不一致的， 只有CCR模型计算所得的投入导向与产出导向的效率是一致的。 对于采用其他 DEA模型时所得结果存在的不一致性， 在实际中也可以将两个角度通过加权综合的方式一起考虑，相关研究可以参考最新的国内外文献。DEA方法不仅能对管理效率进行横向对比，也可以进行纵向、动态的分析，即评价样本的数据可以选择截面数据、时间序列数据或者是面板数据。 面板数据常用的方法为 DEA视窗分析法与 DEA-Malmquist 生产力指数法。7. 当前已有较多的现成的 DEA 软件可以用于求解 DEA模型，例如 DEAP，DEAsolver 以及MyDEA 等。不过这些软件只能求解常见的 DEA模型，对于改进型的 DEA模型，通常需要编程，此时可借助于：EXCEL 的线性规划求解器， Lingo软件以及 Matlab软件等编程软件。习题：1.现有10 家医院，每家医院有 2个投入（医生人数以及护士人数）和 2个产出（门诊病人人数以及住院病人人数），投入产出表如下所示，试用 DEA方法分析这 10家医院的相对有效性。表5—8 十家医院投入产出数据表DMU12345678910投入医生（人）20192527225533313050护士（人）151131160168158255235206244268产出门诊病人（人）10015016018094230220152190250住院病人（人）9050557266908880100100第7节案例分析本节以一个实例来说明数据包络分析方法的使用，评估对象为中国台湾的森林经营。宝岛台湾面积为36000平方公里，台湾森林覆盖面积占全省土地总面积的一半以上，在 1989年以前有 13个林区，主要以保护林地和木材为主要任务。森林的经营具有非盈利性质，了解各个林区的经营效率，检讨投入资源的使用是否有效，是一个值得探讨的问题。以下将以台湾省各林区的效率评估，探讨 DEA 使用的各个步骤，来说明 DEA 方法在实际中的应用。案例分析中主要集中探讨三个部分内容： (1)确认投入产出项， (2)选择恰当的 DEA 模型，(3)结果的分析与解释。使用DEA 方法首先需要选择适当的投入产出项目。根据森林经营的多目标性及其非营利性等特点，其目标包括实质产出（如：木材、野生动物）以及森林效用（如：净化空气，调节气温，美化环境，洁净水源，保持水土，旅游等） ，为了筛选投入产出项以衡量森林经营目标，参考 Kao和Yang(1992) [25]的研究与高强等人 (2003)[26]编写的书籍。筛选的步骤如下：访问林业局的管理层，确立其组织目标及其管理目标。要求受访者确认投入产出项目。进行过程中间研究者从文献以及经验得知的各种投入产出种类列出，以供受访者参考要求受访者确定投入产出的衡量指标。进行过程中间研究者从文献以及经验得知的各种投入产出衡量指标列出，以供受访者参考收集并整理投入产出数据确认投入产出项目及衡量指标并完成数据收集整理后，进一步与受访者探讨，分析其含义。产出项目根据森林一般多目标经营的想法， 例如美国 1960年颁布的 MultipleUse-SustainedYield 法案，森林的功能包括木材生产、野生动物、水源涵养、放牧、森林游乐等五项。通过与林业局的深访，决定台湾省林业经营目标有林产、水源涵养、游乐三项，而采用以下三个指标来衡量：(1) 主产物：木材的产量，度量单位为立方公尺；(2) 平均蓄积：用以衡量水土保持、调节气温、洁净水源等森林功效，以林地林木蓄积量表示，单位为千立方公尺；(3) 游乐：用以衡量森林之游乐功效，单位为游客人数。投入项依据生产经济理论，基本上是土地、资本与劳力，在森林经营中，林地上林木蓄积量亦可视为资本。因此，林业局的相关人员认为因采用如下四项指标来衡量投入，即：(1) 预算：以千元为单位；(2) 原始蓄积：林地之林木蓄积量越大，其林木的生产量也越大，水源涵养效果也越佳，因此以评估期开始的林地为林地蓄积为指标，单位为千立方公尺；(3) 劳动力：投入的人力，单位为人；(4) 土地面积：单位为千公顷。投入产出项经过专家筛选出来以后， 下一步就是投入产出项的数据的收集与整理。 DEA使用不受指标单位的限制，为了避免各年资料一时的波动造成绩效评估的误判，本案例中选取十年的数据(1978-1987)，各指标取十年的平均值，并将货币价值以 1975年的物价销售指数为基底进行平减。因为台湾省在此段时间的林区总数为 13个，未能满足投入产出项目个数两倍之最少受评单位的要求，因此加入四个与十三个林区大致相当的林区经营单位：林区开发处、林业试验所、台大实验林以及中兴实验林。因此获得了 17个受评单位，其数据如下表所示。表5—9 原始数据表投入产出林区面积劳力预算原始蓄积主产物平均蓄积游乐文山60.852704.115.0415.855.1714.57竹东108.46597.99.313.4547.1918.867大甲79.06421.46.358.2721.5710.4833.73大雪山59.66860.112.2810.958.4111.719.64埔里84.52714.339.9339.0412.250峦山127.2859210.4513.3657.1113.810玉山98.886312.158.1442.8112.43399.83楠濃123.148528.8410.866恒春86.372855.358.6239.246.881081.89关山227.2216.15.8724.0444.0827.280玉里146.432054.0815.7637.319.30木瓜173.48774.912.623.039.6323.5341.86兰阳171.112722.714.5117.8419.7318.8684开发处93.651399128.9417.5842.1117.30林试所13.65350.90.911.4219.071.580台大33.521651.730.3813.570.51061.48中兴8.23490.31.593.861.5767.73通过相应的计算，可得CCR模型，BCC模型，超效率模型，原始交叉效率模型，攻击型交叉效率模型，仁慈型交叉效率模型的结果如下表所示。表5—10 各种DEA模型下的效率值及其排序结果林区CCRBCC超效率排序原始交叉效率排序攻击型交叉效率排序仁慈型交叉效率排序文山0.73380.77630.7338170.4779130.4453120.620515竹东1.00001.00001.152670.730250.675750.95381大甲0.89610.91130.8961120.5734110.5366110.78338大雪山1.00001.00001.0418100.4383140.4123150.670612埔里1.00001.00001.196650.796810.739530.91482峦山0.88141.00000.8814130.603790.5491100.734310玉山1.00001.00001.091790.6032100.557690.85514楠浓0.76571.00000.7657150.4784120.4413130.568616恒春1.00001.00001.173860.637480.571980.707311关山1.00001.00001.340940.738240.726540.78777玉里1.00001.00001.117680.789420.766110.86003木瓜0.81931.00000.8193140.4342150.4162140.636314兰阳0.74390.78490.7439160.3444170.3407170.525717开发处0.97461.00000.9746110.4045160.3586160.651713林试所1.00001.00003.062720.654260.638160.73459台大1.00001.000022.256010.639970.599970.82126中兴1.00001.00001.634230.761230.744120.83465从表5—10中可知，CCR模型和BCC模型中存在多个有效决策单元，因此如果分析的目的是实现所有林区的全排序，显然选择这两个模型难以达到要求。而后四种方法：超效率模型，原始交叉效率模型，攻击型交叉效率模型以及仁慈型交叉效率模型可以较好的实现对所有林区的全排序。原始交叉效率模型与攻击型交叉效率模型所得的排序结果最为接近，如玉里林区在原始交叉效率模型中排在第二位，而在攻击型交叉效率模型中排在第一位；而与之类似的林区还有中兴和玉山。差异最大的是埔里林区的排序，在原始交叉效率模型中排第一位，而在攻击型交叉效率模型中排在第三位。而攻击型交叉效率与仁慈型交叉效率所得的排序结果的差异较大，在仁慈型交叉效率中排在第一位的是竹东林区，而竹东林区在攻击型交叉效率模型中排在第五位。 在超效率模型中， 台大林区以相对有效值 22.256遥遥领先于其他林区排在第一位，第二位的林试所以 3.0627排在第二位，中兴林区以 1.6342排在第三位，也就是说后面加入的四个决策单元与林区具备类似功能的单元占据了所有林区排序的前三位，超效率的结果与其他三种可以实现全排序的排序结果均存在较大差异。本部分基本概念 ：数据包络分析 (dataenvelopmentanalysis,DEA); 决策单元 (Decisionmakingunit ，DMU)七 夕，只因有你，总有一些人牵肠挂肚难以忘记，总有一些日子温暖甜蜜最为珍惜从春夏到秋冬，从陌生到熟悉，虽不能时时联系，却总在特别的日子想起你，七夕快乐，我的朋友。七夕，只因有你，因为有你，再苦生活也不觉得累，再大的险阻也无所畏，再大的波折也不担忧，再痛的经历也会忘记，因为有你，我就拥有了整个世界，谢谢你出现在我的生命里。七夕快乐，我的朋友。七夕，只因有你 ，相识，是最珍贵的缘分，牵挂，是最真挚的心动，思念，是最美丽的心情，问候，是最动听的语言，在这七夕到来之际，最美的祝福送给你，七夕快乐，我的朋友。七夕，只因有你 ，雨点轻敲窗，风吹散了梦想，唯有你的模样依旧在脑海里徜徉，夜深人静时，你占满了心房，舍半生轻狂，半世时光，只为拥有一段和你相处的珍贵情缘，七夕快乐，我的朋友。七夕，只因有你，虽然相距很远，但两颗心却紧紧相连虽然不常见面，音容笑貌犹如眼前，悄悄的挟一缕情丝，放飞在炎炎夏日默默的拽一丝牵挂，悬挂在无垠宇宙静静的捎一声问候，盛开在七夕佳节七夕快乐，我的朋友。七夕，只因有你，祝福，是一种真实的心意，是一种甘甜的快乐，是一种浪漫的味道，是一种温馨的记忆，是一种美丽的幸福，更是我们情谊永远不变的纽带，七夕快乐，我的朋友。七夕，只因有你，爱是种体会，即使心碎也觉得甜蜜，爱是种感受，即使痛苦也觉得幸福，爱是种缘分，即使分离也觉得快乐，七夕到了，最真诚的祝福送给你，七夕快乐，我的朋友。七夕，只因有你，愿天下有情人终成眷属，愿单身人士找到爱的方向，愿情侣们找到幸福的天堂，愿夫妻找到温暖的避风巷，愿岁月抚平生活的忧伤，愿爱的花瓣轻舞飞扬，CharnesA,CooperWW,RhodesE.Measuringtheefficiencyofdecisionmakingunits[J].EuropeanJournalofOperationalResearch,1978,2:429-444.[2]CookWD,SeifordLM.Dataenvelopmentanalysis(DEA)-Thirtyyearson[J].EuropeanJournalofOperationResearch,2009,192(1):1-17.[3]LiuJS,LuLYY,LuWM,LinBJY.Dataenvelopmentanalysis1978 –2010:Acitation-basedliteraturesurvey[4]CharnesA,CooperWW.Programmingwithlinearfractionalfunctional[J].NavalResearchLogisticsQuarterly,1962,9:181-185.BankerRD,CharnesA,CooperWW.Somemodelsforestimatingtechnicalandscaleefficienciesindataenvelopmentanalysis[J].ManagementScience,1984,30:1078-1092.BankerRD,ThrallRM.Estimationofreturnstoscaleusingdataenvelopmentanalysis[J].EuropeanJournal[7]AndersenP,PetersenNC.Aprocedureforrankingefficientunitsindataenvelopmentanalysis[J].ManagementScience,1993,39:1261-1264.[8]BankerRD,ChangH.Thesuper-efficiencyprocedureforoutlieridentification,notforrankingefficientunits[J].EuropeanJournalofOperationalResearch,2006,175:1311-1320.[9]SextonTR,SilkmanRH,HoganAJ.Dataenvelopmentanalysis:Critiqueandextensions[J].In:R.H.Silkman,Editor,MeasuringEfficiency:AnAssessmentofDataEnvelopmentAnalysis,Jossey