关于AER上大运河论文，2023-08-15推文DID前沿应用_DID模型时间、空间、混合安慰剂检验新命令didplacebo里面介绍到的命令就用到了大运河数据，在2023年9月DID前沿专题讲座上，也专门讲过该论文，详见【DID异质稳健估计量+合成DID+安慰剂检验汇总】2023年9月DID前沿专题讲座第3课（附板书截图）：可以点击链接查看板书截图
复现AER大运河--SCM方法论文推荐：synth_runner操作应用
本文资源推荐：AER量化历史论文--大论文论文--Table3--笔记
🍓复现AER大运河--更换实证方法--SCM方法：synth_runner操作应用
🌹AER量化历史论文--大运河论文Fig3+附录图1--笔记
🌈复现AER大运河--SCM方法论文推荐：synth_runner操作应用
👉复现AER大运河_交叠DID时间、空间、混合安慰剂检验新命令didplacebo
💎2024Stata暑假班_北京现场/远程_0基础掌握Stata方法+操作+论文+前沿
该文源于：Cao Y, Chen S. Rebel on the Canal: Disrupted Trade Access and Social Conflict in China, 1650–1911[J].American Economic Review, 2022, 112(5): 1555-90.
文章数据资料来源于：https://doi.org/10.1257/aer.20201283
下面小编将带大家回顾一下该文主要内容。
Cao and Chen (2022) 在 American Economic Review (AER) 上发表的关于大运河的论文，主要探讨了1650年至1911年间中国大运河的废弃对社会冲突和贸易准入的影响。
简介
中国大运河作为世界上最大和最古老的人工水道，对区域贸易准入起到了重要作用。然而，随着时间的推移，大运河的废弃对沿线地区的社会和经济产生了深远影响。
Cao and Chen (2022) 的研究旨在分析大运河废弃后，贸易准入中断如何导致社会冲突的增加。
研究方法与数据
数据来源：该研究使用了涵盖262年、575个县的原始数据集。这些数据为大运河废弃前后沿运河县的叛乱程度提供了重要信息。
模型设定：
研究基于标准的双重差分模型（DID），以1826年作为大运河废弃的时点，比较了运河县与临近非运河县在大运河废弃前后叛乱水平的变化。
详见【DID异质稳健估计量+合成DID+安慰剂检验汇总】2023年9月DID前沿专题讲座第3课（附板书截图） (qq.com)：
ssc new---didplacebo命令
AER上大运河文章
时间安慰剂检验
时间安慰剂检验
DID前沿_DID模型时间、空间、混合安慰剂检验新命令didplacebo
1、简介
DID安慰剂检验的Stata新命令didplacebo，可以自动进行DID模型的时间、空间及时空混合安慰剂检验，并提供可视化展示。
didplacebo实现了用于估计差异中的差异(DID)模型的安慰剂检验，其中政策采用可能是同步的或交错的。特别是，didplacebo使用假的治疗时间进行时间安慰剂检验，使用假的处理单位(随机选择)进行空间安慰剂检验，以及使用假的处理单位和时间(都是随机选择)进行混合安慰剂测试。此外，didplacebo为可视化提供了方便的图表。
要执行安慰剂检验，必须指定估计名称estimatename和处理变量名称treatvarname，作为didplacebo命令的输入。didplacebo自动执行estimatename来获取用于DID估计的命令(cmdline)。然后用相同的命令再次估计模型，使用假的处理变量来产生安慰剂效应。Treatvar (treatvarname)指定一个虚拟变量，指示一个单元在特定时期是否被处理。
除了传统的双向固定效应估计(TWFE)之外，要使用异质性稳健估计(例如，Callaway和Sant'Anna, 2021)进行安慰剂检验，可以使用我们的命令xtshuffle, xtrantreat和tofirsttreat。有关详细说明，请参阅Chen, Qi和Yan(2023)。
didplacebo已被确认与官方命令regress、xtreg、areg和xtdidreress以及第三方命令reghdfe兼容。其他命令虽然没有经过专门测试，但也可能与didplacebo兼容。
2、语法格式
didplacebo estimatename , treatvar(treatvarname) [options]
选项含义：
Pbotime (numlist)时间安慰剂测试 Pbounit空间安慰剂试验 Pbomix([1 2 3])混合安慰剂试验 Repeat (int)复制的个数 Seed (int)随机数所使用的种子 frame(framename)存储估计安慰剂的数据框结果 Nofigure不显示图表 Savegraph(前缀，[asis replace])将生成的所有图形保存到当前路径。
注意事项：
必须使用Xtset panelvar timevar来声明面板数据集，参考[XT] xtset。 Estimatename应该是由估计存储保存的估计;参考[R] estimate储存
3、 案例应用
3.1 案例1：标准DID: Cao and Chen:2022
1、导入数据
use cao_chen.dta, clear
2、进行TWFE模型估计
. xtset county year
panel variable: county (strongly balanced)
time variable: year, 1650 to 1911
delta: 1 unit
. reghdfe rebel canal_post, absorb(i.county i.year) cluster(county)
(MWFE estimator converged in 2 iterations)
HDFE Linear regression Number of obs = 140,432
Absorbing 2 HDFE groups F( 1, 535) = 5.23
Statistics robust to heteroskedasticity Prob > F = 0.0226
R-squared = 0.0308
Adj R-squared = 0.0253
Within R-sq. = 0.0002
Number of clusters (county) = 536 Root MSE = 0.3848
(Std. Err. adjusted for 536 clusters in county)
------------------------------------------------------------------------------
| Robust
rebel | Coef. Std. Err. t P>|t| [95% Conf. Interval]
-------------+----------------------------------------------------------------
canal_post | .0380143 .016621 2.29 0.023 .0053639 .0706647
_cons | .0313251 .0007227 43.35 0.000 .0299054 .0327447
------------------------------------------------------------------------------
Absorbed degrees of freedom:
-----------------------------------------------------+
Absorbed FE | Categories - Redundant = Num. Coefs |
-------------+---------------------------------------|
county | 536 536 0 *|
year | 262 0 262 |
-----------------------------------------------------+
* = FE nested within cluster; treated as redundant for DoF computation
. estimates store did_cao_chen
.
end of do-file
结论
漕粮海运导致运河县的叛乱显著增加，DID系数为 0.038, 且在 5% 水平上显著
3、时间安慰剂检验
结论：安慰剂效应的 p 值均大于 0.1，95% 区间均包括 0，故认为满足平行趋势
4、空间安慰剂检验
5、混合安慰剂检验
. didplacebo did_cao_chen, treatvar(canal_post) pbomix(1)
Implementing mixed placebo test for standard DID (version 1) using both fake treatment units
> and times:
-------------------------------------------------------------------------
The number of units randomly | The range within which fake
selected as fake treatment units | treatment times are randomly selected
----------------------------------+--------------------------------------
73 | [1650, 1911]
-------------------------------------------------------------------------
Simulations (500):.........10.........20.........30.........40.........50.........60........
> .70.........80.........90.........100.........110.........120.........130.........140.....
> ....150.........160.........170.........180.........190.........200.........210.........22
> 0.........230.........240.........250.........260.........270.........280.........290.....
> ....300.........310.........320.........330.........340.........350.........360.........37
> 0.........380.........390.........400.........410.........420.........430.........440.....
> ....450.........460.........470.........480.........490.........500
Results of mixed placebo test for standard DID (version 1) using both fake treatment units a
> nd times:
--------------------------------------------------------------
| | P-value
| Coefficient | Two-sided Left-sided right-sided
-----------+-------------+------------------------------------
canal_post | 0.038014 | 0.0100 0.9960 0.0040
--------------------------------------------------------------
Note: (1) The two-sided p-value is the frequency that the absolute values of the placebo
effects are greater than or equal to the absolute value of estimated treatment
effect.
(2) The left-sided (right-sided) p-value is the frequency that the placebo effects
are smaller (greater) than or equal to the estimated treatment effect.
Finished.
.
处理效应估计值位于安慰剂效应分布的右侧尾部
原文稳健性检验使用合成控制法（SCM）进行检验，结果如下图
使用Stata软件进行论文复刻结果如下：
附录对应图为：
