DEA各种模型原理及stata代码实现

DEA各种模型原理及stata代码实现

一、CCR和BCC

1.原理

CCR模型产出导向下的效率通过求解以下规划得出：
$CCR\_TE=max\theta$

$s.t.{\sum }_{k=1}^K{z}_k{y}_{km}\ge y_{km}{\theta }_m,m=1,...,M$

${\sum }_{k=1}^K{z}_k{x}_{kn}\le x_{kn},,n=1,...,N$

$z_k\ge 0$

其中，$({\sum }_{k=1}^K{z}_k{y}_{km},{\sum }_{k=1}^K{z}_k{x}_{kn})$可以理解为前沿面，$(x_{km},y_{km})$为每个决策单元（dmu）的值。

BCC模型在上述规划的约束条件中加入${\sum }_{k=1}^K{z}_k=1$.

2.效率测算stata代码

代码格式如下：

CCR模型对应规模报酬不变crs

dea inputvars = outputvars ,rts(crs) 

BCC模型对应规模报酬可变vrs

dea inputvars = outputvars ,rts(vrs) 

3.Malmquist指数

3.1M指数

Malmquist指数是效率的变化率，简单地想，如果以t期为基期，那么公式为：
$Malmquis{t}_t=\frac{D^t(x^{t+1},y^{t+1})}{D^t(x^t,y^t)}$
其中,$D^t(x^{t+1},y^{t+1})$某一个决策单元在t+1期的生产情况基于t期的前沿面计算的效率（决策单元在t+1期的生产情况可能超出t期的前沿面，因此可能无解）；

$D^t(x^t,y^t)$某一个决策单元在t期的生产情况基于t期的前沿面计算的效率，也就是正常来说的$T{E}_t$

同样的，如果以t+1期为基期，那么公式为：
$Malmquis{t}_{t+1}=\frac{D^{t+1}(x^{t+1},y^{t+1})}{D^{t+1}(x^t,y^t)}$
基期不同值不同，为了解决这个问题，定义：

$Malmquis{t}_{t+1,t}=(Malmquis{t}_t\times Malmquis{t}_{t+1}{)}^{0.5}=M(x^{t+1},y^{t+1},x^t,y^t)$

$=[\frac{D^t(x^{t+1},y^{t+1})}{D^t(x^t,y^t)}\times \frac{D^{t+1}(x^{t+1},y^{t+1})}{D^{t+1}(x^t,y^t)}{]}^{0.5}$

= D t + 1 ( x t + 1 , y t + 1 ) D t ( x t , y t ) × [ D t ( x t + 1 , y t + 1 ) D t + 1 ( x t + 1 , y t + 1 ) D t ( x t , y t ) D t + 1 ( x t , y t ) ] 0.5 =\dfrac{D^{t+1}(x^{t+1},y^{t+1})}{D^t(x^t,y^t)}\times [\dfrac{D^{t}(x^{t+1},y^{t+1})}{D^{t+1}(x^{t+1},y^{t+1})}\dfrac{D^{t}(x^{t},y^{t})}{D^{t+1}(x^t,y^t)}]^{0.5} =Dt(xt,yt)D