对于面板数据，可以采用MALMQUIST-DEA方法来计算全要素生产率变动！
1 数据说明
这里选取了2018-2020年上市工业企业的数据，通过随机前沿模型来计算企业生产效率。以主营业务收入表示企业产出；员工总数表示劳动投入；固定资产净额表示资本投入。
2 录入数据
打开数据编辑窗口，将数据从Excel中复制到Stata，注意粘贴时选择“将第一行作为变量名”
3 数据格式与变量转换
设置面板格式：
reshape long 主营收入 员工总数 固定资产 主营成本 国资背景 是否陆港通 是否出口 资产总计 资产负债率 流动比率 总资产周转率,i(证券代码) j(year)
encode 证券代码,g(id)
xtset id year
变量转换：
g lny=ln( 主营收入)
g lnk=ln( 固定资产 )
g lnl=ln( 员工总数 )
4 估计全要素生产率变动
切换到新的目录：
cd D:\Desktop\dea
由于malmq2命令会将全要素测算结果保存到本地，所以最好提前切换到新的目录，以免找不到结果文件！
安装malmq2命令：
search malmq2
运行malmq2命令：
malmq2 lnk lnl= lny, ort(o) rd seq sav(tfp_result,replace)
ort()参数指定基于投入或产出来计算效率， ort(o)表示基于产出， ort(i)表示基于投入！
rd参数是分解全要素生产率的一方法，也可以指定fgnz来使用另一分解全要素生产率的方法！在malmq2命令中rd和fgnz只能指定1个
seq参数是指定顺序生产技术。除此之外，还有 global、win()、bi参数可以指定不同时期的生产技术。
运行malmq2命令后，在stata目录下面会出现名为tfp_result.dta的文件，将其打开，TFPCH即为所求的全要素生产率变动！存在如下关系：TFPCH=TECH×TECCH×SECH
修改malmq2命令参数：
malmq2 lnk lnl= lny, ort(o) rd global sav(tfp_result1,replace)
运行malmq2命令后，在stata目录下面会出现名为tfp_result1.dta的文件，将其打开，TFPCH即为所求的全要素生产率变动！存在如下关系：TFPCH=TECH×BPC×SECH
修改malmq2命令参数：
malmq2 lnk lnl= lny, ort(i) fgnz win(3) sav(tfp_result2,replace)
运行malmq2命令后，在stata目录下面会出现名为tfp_result2.dta的文件，将其打开，TFPCH即为所求的全要素生产率变动！存在如下关系：TFPCH=TECH×TECCH
malmq2命令还可以估计其他多种全要素生产率的分解形式，这里不再一一演示！