时间跨度：2012-2022年

数据范围：全国30个省份（不含港澳台、西藏）

数据指标：

参考韩文龙等的做法，收集了全部控制变量与稳定性检验所需变量。

样例数据：

处理代码（do）：

//此处根据文件存放位置来调整路径
cd D:\
//导入数据
use 新质生产力水平+被解释变量+控制变量.dta
global cont_var ind inv gov con edu mar
xtset id year

*基准回归
reghdfe ggdp newp ,noabsorb vce(r)
reghdfe ggdp newp ,absorb(id year) vce(r)
reghdfe ggdp newp $cont_var ,absorb(id year) vce(r)
*分指数回归
reghdfe ggdp 新劳动者 $cont_var ,absorb(id year) vce(r)
reghdfe ggdp 新劳动资料 $cont_var ,absorb(id year) vce(r)
reghdfe ggdp 新劳动对象 $cont_var ,absorb(id year) vce(r)
reghdfe ggdp 新技术 $cont_var ,absorb(id year) vce(r)
reghdfe ggdp 生产组织 $cont_var ,absorb(id year) vce(r)
reghdfe ggdp 数据要素 $cont_var ,absorb(id year) vce(r)
*内生性GMM
xtabond ggdp $cont_var ,lag(1) maxldep(2) endogenous(newp,lag(0,2)) twostep
xtabond2 ggdp l.ggdp newp $cont_var , gmmstyle(newp,lag(0 2)) iv($cont_var) nolevel twostep
*机制分析
reghdfe fac newp $cont_var ,absorb(id year) vce(r)
xtabond fac $cont_var ,lag(1) maxldep(2) endogenous(newp,lag(0,2)) twostep
reghdfe icc newp $cont_var ,absorb(id year) vce(r)
xtabond icc $cont_var ,lag(1) maxldep(2) endogenous(newp,lag(0,2)) twostep
*地区异质性
reghdfe ggdp newp  $cont_var if region==1,absorb(id year) vce(r)
reghdfe ggdp newp  $cont_var if region==2,absorb(id year) vce(r)
reghdfe ggdp newp  $cont_var if region==3,absorb(id year) vce(r)
reghdfe ggdp newp  $cont_var if region==4,absorb(id year) vce(r)
*稳健性
winsor newp,gen(w99newp) p(0.01)
winsor newp,gen(w95newp) p(0.05)
reghdfe ggdp w99newp  $cont_var ,absorb(id year) vce(r)
reghdfe ggdp w95newp  $cont_var ,absorb(id year) vce(r)
reghdfe gdp newp  $cont_var ,absorb(id year) vce(r)
reghdfe ggdp newp  $cont_var if hml<=0.5,absorb(id year) vce(r)
reghdfe ggdp newp  $cont_var if hml>0.5,absorb(id year) vce(r)
*空间杜宾
xsmle ggdp newp  $cont_var ,fe model(sdm) wmat(w) type(both) robust nolog effects

下载链接（stata数据集）：https://download.csdn.net/download/li514006030/89447866