查看原文
其他

合成控制法(SCM)安慰剂检验及可视化操作






来源:SYNTH2: Stata module to implement synthetic control method (SCM) with placebo tests, robustness test and visualization


合成控制法简介

合成控制法可视为DID方法的延伸,即考察某一政策冲击对个体影响的效应分析,此时DID方法便失效了,经济学研究中评估某项政策或事件的效应时,使政策可能实施于某个国家,地区或城市。较为简单的方法是考察政策实施前后的时间序列,看所关心的结果如何变化。但结果还可能受其自身时间趋势的影响和其他事件带来的综合影响,为此,常使用“鲁宾的反事实框架”。即假设该地区如果未受到政策干预将会如何。这一方法的困难之处就在于,这一假设是无法观测的。

鉴于反事实无法观测,通常的解决办法是寻找适当的控制组,即在各方面都与受干预地区相似却未受干预的其他地区,以作为处理组的反事实替身。但现实情况中这种控制组的寻找也是相当困难的。为了解决控制组选择的难题,Abadie和Gardeazabal(2003)出了合成控制法”,基本思想为,假设考察仅在北京实施的政策的效果评估,可能可以选择上海,广州,深圳等地,但是实际上他们与北京也不尽相同。虽然无法找到最佳控制地区,但通常可以对中国的若千大城市进行适当的线性拟合,以构造一个“合成控制地区”,并将真实北京”和合成北京”进行对比。

详情阅读:

Stata&一文读懂合成控制法操作及应用(附:安慰剂检验完整代码)


synth2

作为synth实现合成控制方法(SCM)的包装程序,synth2提供了方便的实用程序来自动化使用假治疗单元的空间安慰剂测试、使用假治疗时间的即时安慰剂测试和一次排除一个非零重量控制单元的遗漏鲁棒性测试。Synth2生成一系列用于可视化的图形。命令synth(可从统计软件组件获得)是必需的。

语法格式为:

  synth2 depvar indepvars, trunit(#) trperiod(#) [options]


选项含义为:

  • ctrlunit(numlist)用于捐助池的Ctrlunit (numlist)控制单元

  • preperiod(numlist)  干预发生前的预处理阶段

  • postperiod(numlist)   干预发生时及之后的治疗后时期

  • Xperiod (numlist) indepvar中指定的预测器的平均周期

  • mspeperiod(numlist)均方预测误差(MSPE)应最小化的周期

  • customV(numlist)提供自定义V-Weights,确定变量在预处理期间对结果的预测能力

  • nested 在所有(对角)正半定v矩阵和w权集之间进行搜索的嵌套全嵌套优化过程

  • allopt 如果指定了嵌套,将获得完全健壮的结果

  • placebo([{unit|unit(numlist)} period(numlist) cutoff(#_c)]) 使用假治疗单位的空间安慰剂试验和/或使用假治疗时间的时间安慰剂试验 

操作应用



. use "C:\Users\Metrics\Desktop\smoking.dta"
(Tobacco Sales in 39 US States)

.  xtset state year

Panel variable: state (strongly balanced)
 Time variable: year, 1970 to 2000
         Delta: 1 unit


使用假治疗时间1985和删除协变量cigsale(1988)进行及时安慰剂试验

.  synth2 cigsale lnincome age15to24 retprice beer cigsale(1980) cigs
> ale(1975), trunit(3) trperiod(1989) xperiod(1980(1)1984) nested pla
> cebo(period(1985))
Fitting results in the pre-treatment periods:
--------------------------------------------------------------------------------
 Treated Unit             : California     Treatment Time           :       1989
--------------------------------------------------------------------------------
 Mean Absolute Error      =    1.51510     Number of Control Units  =         38
 Mean Squared Error       =    4.86334     Number of Covariates     =          6
 Root Mean Squared Error  =    2.20530     R-squared                =    0.95253
--------------------------------------------------------------------------------

Predictor balance in the pre-treatment periods:
--------------------------------------------------------------------------------
   Covariate   |  V.weight    Treated    Synthetic Control     Average Control  
               |                        Value          Bias   Value        Bias 
---------------+----------------------------------------------------------------
      lnincome |   0.0000     10.0372      9.8639    -1.73%     9.7892    -2.47%
     age15to24 |   0.0000      0.1815      0.1825     0.55%     0.1814    -0.06%
      retprice |   0.0384     76.2200     76.1523    -0.09%    71.8353    -5.75%
          beer |   0.0000     25.0000     23.0089    -7.96%    23.6947    -5.22%
 cigsale(1980) |   0.9597    120.2000    120.0846    -0.10%   138.0895    14.88%
 cigsale(1975) |   0.0019    127.1000    126.8324    -0.21%   136.9316     7.74%
--------------------------------------------------------------------------------
Note: "V.weight" is the optimal covariate weight in the diagonal of V matrix.
      "Synthetic Control" is the weighted average of control units in the donor pool with
      optimal weights.
      "Average Control" is the simple average of control units in the donor pool with equal
      weights.

Optimal Unit Weights:
--------------------------
     Unit    |    U.weight
-------------+------------
        Utah |     0.3600 
      Nevada |     0.2880 
 Connecticut |     0.1990 
    Colorado |     0.1020 
   NewMexico |     0.0500 
--------------------------
Note: The unit Alabama Arkansas Delaware Georgia Idaho Illinois Indiana Iowa Kansas Kentucky
      Louisiana Maine Minnesota Mississippi Missouri Montana Nebraska NewHampshire
      NorthCarolina NorthDakota Ohio Oklahoma Pennsylvania RhodeIsland SouthCarolina
      SouthDakota Tennessee Texas Vermont Virginia WestVirginia Wisconsin Wyoming in the
      donor pool get a weight of 0.

Prediction results in the post-treatment periods:
-----------------------------------------------------------
 Time | Actual Outcome  Predicted Outcome  Treatment Effect
------+----------------------------------------------------
 1989 |       82.4000            93.0322          -10.6322 
 1990 |       77.8000            89.4297          -11.6297 
 1991 |       68.7000            82.4727          -13.7727 
 1992 |       67.5000            80.6731          -13.1731 
 1993 |       63.4000            79.5929          -16.1929 
 1994 |       58.6000            78.1272          -19.5272 
 1995 |       56.4000            75.6207          -19.2207 
 1996 |       54.5000            74.8372          -20.3372 
 1997 |       53.8000            74.5395          -20.7395 
 1998 |       52.3000            71.1561          -18.8561 
 1999 |       47.2000            71.4380          -24.2380 
 2000 |       41.6000            65.8382          -24.2382 
------+----------------------------------------------------
 Mean |       60.3500            78.0631          -17.7131 
-----------------------------------------------------------
Note: The average treatment effect over the post-treatment periods is -17.7131.

Implementing placebo test using fake treatment time 1985...

Placebo test results using fake treatment time 1985:
-----------------------------------------------------------
 Time | Actual Outcome  Predicted Outcome  Treatment Effect
------+----------------------------------------------------
 1985 |      102.8000           106.1262           -3.3262 
 1986 |       99.7000           103.2850           -3.5850 
 1987 |       97.5000           106.1524           -8.6524 
 1988 |       90.1000            98.4873           -8.3873 
 1989 |       82.4000            96.5237          -14.1237 
 1990 |       77.8000            91.9127          -14.1127 
 1991 |       68.7000            83.7156          -15.0156 
 1992 |       67.5000            81.4730          -13.9730 
 1993 |       63.4000            79.7911          -16.3911 
 1994 |       58.6000            77.9078          -19.3078 
 1995 |       56.4000            76.2193          -19.8193 
 1996 |       54.5000            75.2010          -20.7010 
 1997 |       53.8000            75.1958          -21.3958 
 1998 |       52.3000            71.9437          -19.6437 
 1999 |       47.2000            72.2260          -25.0260 
 2000 |       41.6000            67.1861          -25.5861 
------+----------------------------------------------------
 Mean |       69.6437            85.2092          -15.5654 
-----------------------------------------------------------
Note: The average treatment effect over the post-treatment periods is -15.5654.

Finished.

结果图为:


*实现leave-one-out稳健性测试,并创建一个Stata数据框“california”存储生成的变量


未完待续!


您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存