Stata论文复现|环境规制实现了绿色技术创新的“增量提质”吗——来自环保目标责任制的证据

论文原文和相关数据代码可在《中国工业经济》网站中下载。

陶锋,赵锦瑜,周浩.环境规制实现了绿色技术创新的“增量提质”吗——来自环保目标责任制的证据[J].中国工业经济,2021(02):136-154.

网址：中国工业经济 (ajcass.org)

本公众号刊登的代码为《中国工业经济》公布的源代码，仅作学习使用；若使用作者提供的数据和程序等附件内容，请务必在研究成果上注明引文和数据附件出处。

为了继续致力于打造良好的学术生态，提高论文的可重复性，《中国工业经济》编辑部决定，自 2016 年第 11 期开始，在官方网站 www.ciejournal.org 公开所刊发论文的原始数据和程序、处理数据、案例研究材料、因正文篇幅所限而未刊发的附录等附件。出于保护作者对公开的数据、程序等原创内容的知识产权，编辑部在此郑重提醒，如在科研工作中使用了作者提供的数据和程序等附件内容，请务必在研究成果上注明引文和数据附件出处。如因未明确引用而造成违反学术规范、学术道德等的学术不端问题，作者和编辑部有权予以责任追究。——《中国工业经济》编辑部

Part1主回归

image-20221108143920396

image-20221108143950688

代码如下：

use 数据,clear
xtset id appyear
xi:xtreg number did i.appyear,fe cluster (id)
xi:xtreg number did i.appyear*i.ipc3,fe cluster (id)

xi:xtreg quality did i.appyear,fe cluster (id)
xi:xtreg quality did i.appyear*i.ipc3,fe cluster (id)


Part2橙陈关于交互固定效应的思考

交互FE本质上包含了单个FE，交互FE本质上包含了单个FE，交互FE本质上包含了单个FE。

以第一个主回归为例，我尝试了不同的代码，发现结果一致。

xi:xtreg number did i.appyear*i.ipc3,fe cluster (id)                 *个体FE\年份FE\IPC小类FE\年份-IPC交互FE
xtreg number did i.appyear##i.ipc3,fe cluster (id)                   *个体FE\年份FE\IPC小类FE\年份-IPC交互FE
xtreg number did i.appyear i.ipc3 i.appyear#i.ipc3,fe cluster (id)   *个体FE\年份FE\IPC小类FE\年份-IPC交互FE
xtreg number did i.appyear#i.ipc3,fe cluster (id)                    *个体FE\年份-IPC交互FE
reghdfe number did,absorb(id appyear#ipc3) cluster (id)              *个体FE\年份-IPC交互FE
reghdfe number did,absorb(id ipc3 appyear#ipc3) cluster (id)         *个体FE\IPC小类FE\年份-IPC交互FE
reghdfe number did,absorb(id appyear appyear#ipc3) cluster (id)      *个体FE\年份FE\年份-IPC交互FE
reghdfe number did,absorb(id appyear ipc3 appyear#ipc3) cluster (id) *个体FE\年份FE\IPC小类FE\年份-IPC交互FE

也可以采用egen命令产生交互项。

egen appyear_ipc3 = group(appyear ipc3)
xtreg number did i.appyear_ipc3,fe cluster (id)     *一定要i.appyear_ipc3,否则将appyear_ipc3视为连续变量


我发现reghdfe的代码运行结果更加清爽，由于交互项太多了，用xtreg会展示由于多重共线性删除了哪些项目，但reghdfe不显示多重多线性的处理过程，直接给结果，但最终结果xtreg和reghdfe一致。

image-20221108151357041

image-20221108151323521

Part3平行趋势检验

image-20221108152356176

*动态效应
use 数据,clear
xtset id appyear
gen Dyear=appyear-2006
gen before4=(Dyear==-4 & treat==1)
gen before3=(Dyear==-3 & treat==1)
gen before2=(Dyear==-2 & treat==1)
gen before1=(Dyear==-1 & treat==1)
gen current=(Dyear==0 & treat==1)
gen after1=(Dyear==1 & treat==1)
gen after2=(Dyear==2 & treat==1)
gen after3=(Dyear==3 & treat==1)
gen after4=(Dyear==4 & treat==1)
gen after5=(Dyear==5 & treat==1)
gen after6=(Dyear==6 & treat==1)
gen after7=(Dyear==7 & treat==1)
gen after8=(Dyear==8 & treat==1)
gen after9=(Dyear==9 & treat==1)
gen after10=(Dyear==10 & treat==1)

xi:xtreg number before4 before3 before2 current after1 after2 after3 after4 after5 after6 after7 after8 after9 after10 i.appyear*i.ipc3,fe cluster(id) 
xi:xtreg quality  before4 before3 before2 current after1 after2 after3 after4 after5 after6 after7 after8 after9 after10 i.appyear*i.ipc3,fe cluster(id)