数据全部来源于国泰安：
mainincome 为企业主营业务收入
labor 为企业员工数量
median 购买商品接受劳务支付的现金（也称中间投入）
netassest 为净固定资产

计算代码如下：
generate lnY=ln(mainincome)
generate lnL=ln(labor)
gen lnM=ln(median)
gen lnK=ln(netassest)

*去除缺失值
foreach i in lnY lnL lnM lnK {
    drop if `i'==.
}
*缩尾处理
winsor2 lnY lnL lnM lnK, cut(1 99) replace by(year)

*国泰安上面下载的企业代码不是数字格式，这里是把字符转化为数字
destring stkcd,replace

*计算全要素生产率
levpet lnY, free(lnL) proxy(lnM) capital(lnK) i(stkcd) t(year)
predict tfp_lp,omega
gen tfp=ln(tfp_lp)

*注意：最后得到的tfp是已经取过对数的

参考资料包括：https://zhuanlan.zhihu.com/p/264359691?utm_id=0 知乎帖子
连享会： https://www.lianxh.cn/details/13.html
b站：https://www.bilibili.com/read/cv24660865/
论文：中国工业企业全要素生产率估计：1999-2007 鲁晓东、连玉君

帖子的方法如果存在错误欢迎批评指正，谢谢！