1 计算SFA的软件FRONT41.exe
有多种软件支持SFA方法，这里采用FRONT软件来计算我国的生产效率！
FRONT软件是绿色软件，不用安装，直接解压FRONT软件压缩包后，通过设置数据文件和指令文件后，双击front41.exe即可运行计算生产率。
当然front软件的缺点是不够智能，很多数据转换步骤需要在其他软件中完成，且测算后的效率数据存放在txt文件中，不方便效率数据的后续分析！除此之外，front软件只支持无效率项的半正态分布形式和截断正态分布形式！
2 中国生产效率的测算：半正态分布形式
2.1 数据说明
资本存量采用本专题第一部分测度出的数据！
产出增加值以GDP的数据！
劳动采用全国就业人数的数据。在有些文献中会将就业人数乘以劳动者平均受教育年限来表示劳动质量，这里不再演示。
2.2 数据预处理
GDP的数据为现价数据，需要进行平减，消除价格变动的影响。这里同时下载了以1078年为基期的GDP指数。由于资本存量的基期是1990年，所以这里也将gdp指数的基期改为1990年（GDP指数序列除以1990年的指数值），然后将现价GDP除以GDP指数，得到平减后的GDP。
最终得到用于计算全要素生产率的产出、资本、劳动数据
2.3 数据取对数
将整理好的产出、资本、劳动数据取对数处理
2.4 整理成front数据排列要求
front中数据前两列要求生产者序号和时间序号两列，即使是截面数据，也需要这两列。
本例中是截面数据（时间序列），所以时间序列所有数据都为1.
2.5 复制粘贴到fron41.exe目录
在front41.exe文件所在目录下，新建名为eff1-dta.txt的文件。
将Excel中个体序列、时间序列、取对数后的lny、lnk和lnl这5列数据复制粘贴到eff1-dta.txt文件。
2.6 设置指令文件
在deap.exe文件所在目录下，新建名为eff1-ins.txt的文件。
双击打开eff1-ins.txt文件，进入如下设置：
第1行：1是选择SFA模型形式，1表示一般模型（无效率项不会受到外界因素影响）；2表示无效率项受到外界因素影响。这里选择1.
第2行：eff1-dta.txt是刚才粘贴到数据的文件名
第3行：eff1-out.txt是DEAp计算出生产率后需要保持到的文件，文件名可以自己起
第4行：1是函数类型，1表示生产函数，2表示成本函数！这里应该选择1.
第5行：y表示样本是否取对数，y表示已取对数，n表示没有取对数。本例中数据在Excel中已经取对数处理，所以设置为y
第6行：31是样本数量，本例中数据是1990年至2020年共31个样本
第7行：1是时期数，本例中把时间序列当成截面数据来使用，相应的时间数应该设为1
第8行：31是样本总数，本例中为31
第9行：2是是投入的序列数（解释变量），本例中有2列投入，即资本和劳动序列，所以设置为2
第10行：n是SFA模型参数设置。如果第一行选择了“1”，那么第10行设置y表示截断正态分布（mu不等于0）；设置n表示半正态分布（mu等于0）。如果第一行选择了“2”（此时无效率项受到外界影响，即mu不等于0），那么第10行设置y表示影响mu的方程中不存在截距项；设置n表示影响mu的方程中包含截距项。
本例第一行选择了1，这里第10行设置n表示选择半正态分布形式！
第11行：n也是SFA模型参数设置。如果第一行选择了“1”，那么第11行设置y表示效率虽时间变化（eta不等于0，适用于面板数据）；设置n表示效率不随时间变化（eta等于0）。如果第一行选择了“2”，那么第10行应该数据影响mu的方程中解释变量的个数。
本例第一行选择了1，这里第11行设置n表示效率不随时间变化！
第12行：n表示估计SFA模型时是否设置初始值，一般选择n。
2.7 运行fron41.exe
双击运行front41.exe，在出现的新窗口中，首先输入f，点击回车，然后接着输入eff1-ins.txt指定需要运行的指令文件，点击回车，即可开始SFA分析！得到生产效率的计算结果。
2.8 生产率结果
在运行fron41.exe程序后，在该目录下会生成eff1-out.txt文件，即生产率估计结果文件。
注意！由于front41.exe程序不提供报错机制，因而如果运行front41.exe后没有出现eff1-out.txt文件，那么需要仔细检查eff1-dta.txt和eff1-ins.txt文件，查看是否设置正确！
双击打开eff1-out.txt文件，文件中给出了似然估计法得到的生产函数的回归系数，以及地区的技术效率，注意不是全要素生产率！
3 中国生产效率的测算：截断正态分布形式
这里进一步计算截断正态分布形式的技术效率
3.1 复制数据和指令文件
由于原始数据不变，所以这里复制front41.exe目录中的eff1-dta.txt和eff1-ins.txt文件，并命名为eff2-dta.txt和eff2-ins.txt.
3.2 修改指令文件
双击打开eff2-ins.txt文件，将第10行的n改为y，同时将第2行和第3行修改为eff2-dta.txt和eff2-out.txt
3.3 重新运行front41.exe
双击运行front41.exe，在出现的新窗口中，首先输入f，点击回车，然后接着输入eff2-ins.txt指定需要运行的指令文件，点击回车，即可开始SFA分析！得到生产效率的计算结果。
3.4 生产率结果
在运行fron41.exe程序后，在该目录下会生成eff2-out.txt文件，即生产率估计结果文件。
4 中国生产效率的测算：无效率项受到外界因素影响
这里进一步研究当无效率项受到外界影响状态下的生产效率估计。
4.1 数据说明
纯粹为了操作演示，我们这里假定无效率项会受到通货膨胀的影响。由于本例中产出数据和资本存量数据已经进行过价格平减，所以并不包含价格变动因素。然后以cpi来反应通货膨胀的程度！
4.2 复制粘贴到fron41.exe目录
在front41.exe文件所在目录下，新建名为eff3-dta.txt的文件。
将Excel中个体序列、时间序列、取对数后的lny、lnk、lnl和cpi这6列数据复制粘贴到eff3-dta.txt文件。
4.3 设置指令文件
在front41.exe文件所在目录下，新建名为eff3-ins.txt的文件；或者从eff1-ins.txt复制过来。双击打开eff3-ins.txt文件，进入如下设置：
第一行修改为2，第二行和第三行修改为eff3-dta.txt和eff3-out.txt，第10行n表示无效率项中不包含截距项（也可以选择y，包含截距项），第11行修改为1，即本例中只有通货膨胀1个影响无效率项的变量。
4.4 重新运行front41.exe
双击运行front41.exe，在出现的新窗口中，首先输入f，点击回车，然后接着输入eff3-ins.txt指定需要运行的指令文件，点击回车，即可开始SFA分析！得到生产效率的计算结果。
4.5 生产率结果
在运行fron41.exe程序后，在该目录下会生成eff3-out.txt文件，即生产率估计结果文件。
回归结果中delta1是无效率项的解释变量（通货膨胀）回归系数！
5 中国生产效率的测算：超越对数生产函数形式
前面3个例子均采用了柯布道格拉斯生产函数，这里我们进一步采用超越对数生产函数形式来测算生产率！
5.1 数据说明
由于超越对数生产函数要求投入变量的2次项，所以在Excel中整理了lnk2、lnk、lnl2、lnl和lnklnl等变量
5.2 复制粘贴到fron41.exe目录
在front41.exe文件所在目录下，新建名为eff4-dta.txt的文件。
将Excel中个体序列、时间序列、lny、lnk2、lnk、lnl2、lnl和lnklnl这8列数据复制粘贴到eff4-dta.txt文件。
5.3 设置指令文件
在front41.exe文件所在目录下，新建名为eff4-ins.txt的文件；或者从eff1-ins.txt复制过来。双击打开eff4-ins.txt文件，进入如下设置：
第一行修改为1，第二行和第三行修改为eff4-dta.txt和eff4-out.txt，第9行改为5，因为有lnk2、lnk、lnl2、lnl和lnklnl这5个解释变量。
5.4 重新运行front41.exe
双击运行front41.exe，在出现的新窗口中，首先输入f，点击回车，然后接着输入eff4-ins.txt指定需要运行的指令文件，点击回车，即可开始SFA分析！得到生产效率的计算结果。
5.5 生产率结果
在运行fron41.exe程序后，在该目录下会生成eff4-out.txt文件，即生产率估计结果文件。
由于超越对数生产函数有5个变量，所以加上截距项，回归系数总共有6个！