多时点DID保姆级教程(下)-安慰剂检验

阅读全文：多时点DID保姆级教程(下)-安慰剂检验 (lianxh.cn)

作者：李闯 (新疆大学)
邮箱：chuangli0808@163.com

---

安慰剂检验核心思想在于虚构处理组或者虚构政策时间进行估计，如果不同虚构方式下的估计量的回归结果依然显著，那么就说明说明原来的估计结果很有可能出现了偏误。现有的关于 DID 的安慰剂检验主要有以下几种方式：

1. 改变政策时间

改变政策发生时点，具体包括政策时点提前、滞后两种情况，还包括将政策时间随机化这种更一般化的处理方式 (多用于政策时点随机提前)。

政策时点提前，若回归系数不显著，则说明是该政策起作用，反向验证了结论的稳健性。政策时点滞后，系数依旧显著，一般而言绝对值增大，表明政策效果明显，并且具有一定的持续性。

1.1 政策时点提前、滞后一定期数

数据来自石大千等 (2018)《智慧城市建设能否降低环境污染》这篇文章，严格意义来说本文只使用了 2012 年的一期智慧城市试点，不属于多期DID，但是方法相同。

use smart_city, clear
merge m:1 id using smart_pt //匹配政策变量，生成政策时间变量
drop _merge
replace pt=0 if pt==.
g did_1=0
replace did_1=1 if pt==2012&year>=2011 //政策时点提前1期
g did_2=0
replace did_2=1 if pt==2012&year>=2010 //政策时点提前2期
g did1=0
replace did1=1 if pt==2012&year>=2013 //政策时点滞后1期
g did2=0
replace did2=1 if pt==2012&year>=2014 //政策时点滞后2期
xtreg lnrso did_1 lnrgdp lntgdp lninno lnurb lnopen lnss,fe
xtreg lnrso did_2 lnrgdp lntgdp lninno lnurb lnopen lnss,fe
xtreg lnrso did1 lnrgdp lntgdp lninno lnurb lnopen lnss,fe
xtreg lnrso did2 lnrgdp lntgdp lninno lnurb lnopen lnss,fe //以人均废气排放量为例

数据表明在假定政策时点提前 1 期、2 期后，智慧城市政策对人均废气排放量均无显著抑制作用，从侧面证实了智慧试点政策的环境友好效应。同时在政策分别滞后 1 期、2 期时，系数依旧显著为负，且绝对值增大，表明政策的促进环境友好发展作用在加强。

. lxhuse smart_city, clear
. merge m:1 id using https://file.lianxh.cn/data/s/smart_pt.dta // 匹配政策变量，生成政策时间变量
. drop _merge
. replace pt=0 if pt==.
. g did_1=0
. replace did_1=1 if pt==2012&year>=2011 // 政策时点提前1期
. g did_2=0
. replace did_2=1 if pt==2012&year>=2010 // 政策时点提前2期
. g did1=0
. replace did1=1 if pt==2012&year>=2013  // 政策时点滞后1期
. g did2=0
. replace did2=1 if pt==2012&year>=2014  // 政策时点滞后2期

. xtset id year
. g lntgdp=lnrgdp*lnrgdp
. gen lnww=log(ww)
. gen lnso=log(so)
. g lnrww=log(ww/pop)
. xtreg lnrso did_1 lnrgdp lntgdp lninno lnurb lnopen lnss,fe
. est  store m1
. xtreg lnrso did_2 lnrgdp lntgdp lninno lnurb lnopen lnss,fe
. est  store m2
. xtreg lnrso did1 lnrgdp lntgdp lninno lnurb lnopen lnss,fe
. est  store m3
. xtreg lnrso did2 lnrgdp lntgdp lninno lnurb lnopen lnss,fe
. est  store m4

. esttab m1 m2 m3 m4, star(* 0.1 ** 0.05 *** 0.01) b(%6.3f) ///
>     order(did_1 did_2 did1 did2) compress nogap

--------------------------------------------------------------
                 (1)          (2)          (3)          (4)   
               lnrso        lnrso        lnrso        lnrso   
--------------------------------------------------------------
did_1         -0.051                                          
             (-0.91)                                          
did_2                      -0.091                             
                          (-1.63)                             
did1                                    -0.114*               
                                       (-1.86)                
did2                                                 -0.128*  
                                                    (-1.84)   
lnrgdp         2.695***     2.621***     2.600***     2.619***
              (5.79)       (5.66)       (5.63)       (5.69)   
lntgdp        -0.137***    -0.133***    -0.132***    -0.133***
             (-5.96)      (-5.82)      (-5.79)      (-5.86)   
lninno        -0.037*      -0.037*      -0.037*      -0.038** 
             (-1.90)      (-1.90)      (-1.89)      (-1.97)   
lnurb          0.172***     0.172***     0.174***     0.174***
              (5.06)       (5.06)       (5.11)       (5.11)   
lnopen         0.030        0.030        0.030        0.031   
              (1.23)       (1.25)       (1.27)       (1.28)   
lnss           0.033        0.033        0.026        0.027   
              (0.24)       (0.25)       (0.19)       (0.20)   
_cons        -14.771***   -14.430***   -14.326***   -14.415***
             (-6.14)      (-6.03)      (-6.00)      (-6.06)   
--------------------------------------------------------------
N               2167         2167         2167         2167   
--------------------------------------------------------------
t statistics in parentheses
* p<0.1, ** p<0.05, *** p<0.01

 阅读全文：多时点DID保姆级教程(下)-安慰剂检验 (lianxh.cn)