三阶段DEA模型操作步骤笔记

参考B站视频三阶段DEA操作步骤，做个简单的笔记，便于用的时候复习 https://www.bilibili.com/video/BV1F4411n75p?p=1

使用到的软件：DEAP2.1、Frontier4.1、excel

根据数据包络法对DMU的内在要求，即DMU的个数要大于等于所有的投入变量和产出变量的两倍

这里原始数据是截取博主给的数据中的前10个（主要是记录操作步骤） ：

3个投入，2个产出，2个环境变量

运用DEAP2.1软件进行效率分析

如何使用deap软件

（1）打开123.DTA，用于存储数据，产出在前、投入在后

将数据复制到 123.dta

crste = TE vrste = PTE scale = SE

TE表示综合效率，PTE表示纯技术效率，SE表示规模效率，满足TE=SExPTE；drs表示规模报酬递减，irs表示规模报酬递增，- 表示规模报酬不变。

deap 2.1软件分析过程及结果解释

第一阶段分析出来的投入产出松弛变量可能受到外部环境因素、随机误差以及内部管理因素影响，通过随机前沿分析方法（Stochastic Frontier Analysis，简称SFA）对上述因素进行测算并将其影响剔除。

（1）

然后

松弛变量S = 原始投入值 - 目标投入值

关于 这样算出来的松弛变量与DEA结果带的松弛变量结果 SUMMARY OF INPUT SLACKS 部分差距很大的原因

四阶段DEA松弛变量及模型调整问题 [经济学] 请教：DEA分析中投影变量与松弛变量之间的关系 https://bbs.pinggu.org/thread-3347318-1-1.html

其次，对效率为1的DUM，松弛变量可以令为0

注：因原始数据为残缺的，这些结果数值不必在意，主要是操作步骤

理论部分我参考的是：研发投入产出效率的国际比较研究—基于三阶段DEA模型分析 以及 关于三阶段DEA模型的几点研究

实际操作运用Frontier4.1软件+ excel

投入的松弛变量作为因变量，K个环境变量作为自变量

注意：SAF只能做单产出分析，即每次只能调整一个投入变量，若有 m 个投入变量，则需要进行 m 次调整。

首先对投入1进行调整，将excel中整理好的数据复制粘贴进 EG1.DTA 中

如下图

注意顺序

第一列是评价体系的序号；第二列是时期t；.第三列是因变量；.第四列之后是K全自变量

接下来设置 EG1.ins

这里根据自己的实际情况设置参数

我这里选择

第一个 y ：假设U的分布。y表示截断分布，n表示半正态分布. 第二个 y ：表示时变模型，n表示非时变模型。 选择 n

其他设置保持不变.

设置好之后，运行front41.exe 输入 ： F EG1.ins

回车之后，文件夹中会新生产一个eg1.out 文件

打开新生成的 out 文件

在使用Frontier4.1分析投入差额值时，估计回归系数使用的是最大似然法。

检验 参数是否显著（t检验是检验参数是否显著为0，参数是否有统计学意义） 以及 判断LR是否通过检验（单边误差似然检验值，为了判断使用SFA模型是否合理）

对于 t 值（ t-ratio） ，一般t值大于2，通常都显著，但t值并不是最终判定标准，还需要结合 p值 一起看（计算P值的方法在下方t值分布表中有）。

以及 LR 值，单边似然比统计量LR=-2×ln (L (H0) /L (H1) )

t检验的自由度是样本数-待估参数-1，LR检验的自由度见结果：with number of restrictions

环境变量数

注:***、**、*分别表示在1%、5%、10%的显著性水平上通过正态检验

t分布表如下 注：对 α < 0.5 ，有 t α < − t 1 − α \text{注：对}\alpha