1 数据说明
工业企业数据库每年包含了数十万条企业样本数据，是研究中国工业企业发展的重要微观数据！但是工业企业数据库样本、变量非常多，且不同年份指标还存在差异，给分析带了较大困难！这里我们已2011-2013年3年的工业企业数据库为例来演示工业企业全要素生产率的计算！
参考相关文献，企业投入产出的指标做如下近似处理
-
产出增加值-用主营业务收入近似表示
-
资本-用固定资产合计近似表示
-
总产出-用工业总产值近似表示
-
劳动投入-用就业人口近似表示
-
中间投入-用总产出减去产出增加值近似表示
这里采用Stata16软件来分析。2011-2013年工业企业数据库已经匹配并保存成stata数据格式！（数据及变量中文对照文件放置于文末！）
注意:在计算全要素生产率时，产出数据有两种形式，一种是增加值，一种是总产值。如果用阐述产出用增加值，那么生产函数中不需要加入中间投入变量，但如果产出用总产值，那么生产函数中需要加入中间投入变量！
2 数据预处理
///第一部分：加载工业企业数据.dta文件
use D:\Desktop\工业企业数据.dta
///第二部分：删除缺失值及不符合常理的样本
//2.1 删除工业总产值、销售产值、资产总计、固定资产合计、从业人数为负、小于零或缺失的企业
drop if (indo==0 | indo<0 )
drop if (isov==0 | isov<0 )
drop if (toasA==0 | toasA<0 )
drop if (tofa==0 | tofa<0 )
drop if (fias==0 | fias<0 )
drop if (nvoe==0 | nvoe<0 )
//2.2 删除企业成立时间在1949年之前的企业
drop if opye<1949
//2.3 删除企业销售利润率超过100%的企业
drop if (oppr/isov)>1
//2.4 删除资产总计小于固定资产的企业
drop if toasA<tofa
//2.5 剔除最大和最小的1%变量
egen p01=pctile(nvoe) ,p(1)
g temp=0
replace temp=1 if nvoe<=p01
bys id : egen x=total(temp)
drop if x>=1
drop p01 temp x
egen p99=pctile(nvoe) ,p(99)
g temp=0
replace temp=1 if nvoe>=p99
bys id : egen x=total(temp)
drop if x>=1
drop p99 temp x
///第三部分：计算工业企业投入和产出变量
//3.1 删除id和year重复的样本
duplicates drop id year,force
//3.2 对企业以现价表示的投入和产出数据进行价格平减
g ppi_index=100 if year==2011 /*用ppi指数对工业总产值、工业销售产值进行平减*/
replace ppi_index=98.3 if year==2012
replace ppi_index=96.4 if year==2013
g assets_index=100 if year==2011 /*用固定资产价格指数对固定资产进行平减*/
replace assets_index=101.1 if year==2012
replace assets_index=101.4 if year==2013
g real_add=mabi/ppi_index*100 /*用主营业务收入近似表示产出增加值*/
g real_k=tofa/assets_index*100 /*用固定资产合计近似表示固定资产*/
g real_total=indo/ppi_index*100 /*用工业总产值近似表示总产出*/
//3.3 投入产出变量取对数处理
g lny=ln(real_add)
g lnl=ln(nvoe) /*就业人口近似表示劳动投入*/
g lnk=ln(real_k)
g lnm=ln(real_total-real_add) /*总产出减去产出增加值近似表示中间投入*/
g lnt=ln(real_total)
//3.4 设置面板数据格式
xtset id year
3 LP法计算全要素生产率
运行levpet命令，估计增加值形式下的企业全要素生产率：
levpet lny, free(lnl) proxy(lnm) capital(lnk) valueadded reps(5)
predict tfp_lp1, omega
该命令采用了产出为增加值的形式，其中free() proxy() capital()为必填字段；valueadded表示增加值形式。
运行levpet命令，估计总产值形式下的企业全要素生产率：
levpet lnt, free(lnl) proxy(lnm) capital(lnk) revenue reps(5)
predict tfp_lp2, omega
该命令采用了产出为总产值的形式，revenue表示增加值形式。