1 数据说明
这里选取了2020年上市工业企业的数据，通过随机前沿模型来计算企业生产效率。以主营业务收入表示企业收入；主营业务成本表示企业生产成本；员工总数表示劳动投入；固定资产净额表示资本投入。
2 录入数据
打开数据编辑窗口，将数据从Excel中复制到Stata，注意粘贴时选择“将第一行作为变量名”
3 变量转换
将企业投入产出变量取对数：
g lny=ln( 主营收入)
g lnk=ln( 固定资产 )
g lnl=ln( 员工总数 )
g lnc=ln( 主营成本 )
g n=_n
4 生成dea要求的数据形式
切换目录：
cd D:\Desktop\dea
由于dea命令会将全要素测算结果保存到本地，所以最好提前切换到新的目录，以免找不到结果文件！
生成dmu变量：
encode 证券代码,g(dmu)
dea命令需要指定企业标识，且必须命名为“dmu”
生产投入、产出变量：
g o_lny=lny
g i_lnk=lnk
g i_lnl=lnl
dea命令中，投入变量必须以i开头，产出变量必须以o开头
删除确实值：
drop if lny==. | lnk==. | lnl==.
如果存在缺失值，比如某个企业资本投入数值缺失，dea依然会计算该企业的生产效率，且会高估该企业的生产效率！
5 计算生产效率
安装dea命令：
search dea
运行dea命令：
dea i_lnk i_lnl = o_lny, rts(vrs) ort(out) stage(2) saving(dea_result)
dea命令中i_lnk、i_lnl和o_lny分别为投入和产出变量，rts()指定规模报酬是否变化，vrs参数值是报酬可变，crs参数是报酬不变；ort()指定基于投入或产出方向，out参数是产出方向，in参数是投入方向。stage()是dea计算的方法，1是一步法，2是二步法。
dea命令运行速度非常慢，这也是stata中dea命令的缺点！
在运行完dea命令后，会在stata目录下面出现名为dea_result.dta的文件，将其打开。可以看出dea命令分别计算了每个企业规模报酬不变下的技术效率（CRS_TE）、规模报酬可变下的技术效率（VRS_TE），以及规模效率（SCALE）