1 计算DEA的软件DEAP.exe
对于省市层面的生产率，只要要本量不是特别大，依然可以用DEAP.exe来计算。
2 省市生产率测算
2.1 数据说明及处理
这里仍然采用OLS方法中地区全要素生产率中的数据，即用1990年至2017年中国31个地区的gdp、资本和劳动力数据来估计各地区的生产率！
本例的数据属于面板数据，如果仅仅想估算技术效率、或者配置效率、规模效率，则可以将面板数据看成截面数据，直接用截面数据的DEA估计方法来测算，如前面DEA方法中“中国生产率测算”的例子。**需要注意的是如果想得到所有年份所有地区对应的生产率效率，则设置指令文件时第4行个体数应填868（31×28），第4行时期数应填1；如果第3行个体数填31，第4行时期数填28，则deap只会计算31个地区28年的平均生产效率！**具体操作这里不再演示！
对于面板数据，可以采用MALMQUIST-DEA方法来估计全要素生产率！我们主要演示采用DEA来测算省市全要素生产率的操作过程！
2.2 数据转换
deap软件在处理面板数据时，需要将面板数据先按照时间顺序，在按照个体顺序进行排序，如下图所示：
在Excel中，将数据按照年份和地区进行排序，然后将gdp、资本存量和劳动力数据取对数处理！
2.3 制粘贴到DEAP.exe目录
在deap.exe文件所在目录下，新建名为eff1-dta.txt的文件。
将Excel中取对数后的gdp、资本存量和就业人口数据复制粘贴到eff1-dta.txt文件。注意只用复制粘贴数据列，不用复制粘贴年份和地区这两列！
2.4 设置指令文件
在deap.exe文件所在目录下，新建名为eff1-ins.txt的文件。
双击打开eff1-ins.txt文件，进入如下设置：
第1行：eff1-dta.txt是刚才粘贴到数据的文件名
第2行：eff1-out.txt是DEAp计算出全要素生产率后需要保存到的文件，文件名可以自己起
第3行：31是截面个体的数量，本例中截面为31个省市，共31个样本。
第4行：28是时期数，本例时间范围为1990-2017年，共28年数据。
第5行：1是产出的序列数，本例中只有1列产出，即gdp序列，所以设置为1
第6行：2是投入的序列数，本例中有2列投入，即资本和劳动序列，所以设置为2
第7行：0是计算全要素的前沿方向，0表示基于投入来计算生产前沿；1表示基于产出来计算生产前沿。一般来说，基于投入或产出来计算生产前沿，最终计算的生产率相差不大。经验的作用是首先考虑基于投入来计算生产前沿；而如果企业的投入资源不受限制，尽可能增加产出，此时应该使用基于产出来计算生产前沿。本例中设置为0，采用基于投入来计算生产前沿！
第8行：0是规模报酬的形式，0表示规模报酬不变，1表示规模报酬可变。一般情况下采用规模报酬不变的方式来计算生产率，需要计算规模效率时，再采用规模报酬可变的方式来计算生产率。但是当计算全要素生产率时，dea程序默认会同时计算规模报酬不变和可变两种情况下的效率，进而计算全要素生产率，所以选择0和1的结果是一样的！
第9行：2是计算生产率的方法。0表示采用多阶段方法来计算生产率，3和4分别表示采用1阶段和2阶段方法来计算生产率。在1阶段、2阶段和多阶段方法中，多阶段方法优势明显，故0，3，4中优先选择0。2表示MALQUIST-DEA方法，适用于面板数据，1表示成本DEA方法，适用于有投入价格的数据。只有MALQUIST-DEA方法能够计算全要素生产率，所以本例选择2.
2.5 运行DEAP.exe
双击运行deap.exe，在出现的新窗口中，输入eff1-ins.txt指定需要运行的指令文件，点击回车，即可开始DEA分析！得到生产效率的计算结果。
2.6 生产率结果
在运行deap.exe程序后，在该目录下会生成eff1-out.txt文件，即生产率估计结果文件。
注意！由于deap.exe程序不提供报错机制，因而如果运行deap.exe后没有出现eff1-out.txt文件，那么需要仔细检查eff1-dta.txt和eff1-ins.txt文件，查看是否设置正确！
另外，由于面板数据计算时间较长，需要等一小会儿才能出现eff1-out.txt文件！
双击eff1-out.txt打开，MALMQUIST INDEX SUMMARY处列出了从第2年开始的生产率测算结果，即分别列出了各年份31个地区的生产率，firm的次序与原始地区的排序一致。效率包括技术效率变动（effch）、技术进步变动（techch）、纯技术效率变动（pech）、规模效率变动（sech）和全要素生产率变动（tfpch），其中有如下关系：
全要素生产率变动=技术效率变动×技术进步变动
技术效率变动=纯技术效率变动×规模效率变动