識別策略
清朝末期变泄，清政府與西方列強簽訂了一系列不平等條約，開放沿江沿海等城市作為通商口岸即是不平等條約的主要內(nèi)容之一恼琼。賈瑞雪老師（2014）將近代通商口岸的設(shè)置作為一項準(zhǔn)自然實驗妨蛹，評估了通商口岸對中國近現(xiàn)代人口和經(jīng)濟發(fā)展的長期影響。
從1840年至1910年晴竞，中國一共被迫開放了40多個通商口岸蛙卤。不同通商口岸開放的時間（政策時點）是不同的，例如，廣州颤难、福州神年、廈門、寧波和上海是在1842年《南京條約》后開放的行嗤，漢口瘤袖、九江、南京昂验、鎮(zhèn)江等城市是在1858年《天津條約》后開放的捂敌，而蘇州、杭州等城市是在1895年《馬關(guān)條約》后開放的......
如上表所示既琴，蘇州府是在1896年被開放為通商口岸的占婉，所以蘇州府的政策分期變量period在1896年之后取值為1，之前取值為0甫恩；而同處蘇杭地區(qū)的嘉興府則一直沒有被開放為通商口岸逆济，所以嘉興府的政策分組變量treat和政策分期變量period一直取值為0。
事實上，我們可以發(fā)現(xiàn)交互項treat×period的取值和政策分期變量period的取值是一毛一樣的，所以在多期DID中弦蹂，我們其實就沒有必要去生成什么交互項背犯，只需用一個政策虛擬變量予以替代就可以了翁垂，用以表示地區(qū)i在t期是否實施政策。當(dāng)然，如果為了便于理解的話，可以嘗試去生成交互項岛马，結(jié)果都一樣滴！
Stata操作
【1】設(shè)置環(huán)境
cd?D:\experiments（位置一般為數(shù)據(jù)所在文件夾）
多期雙重差分法（DID）的Stata操作可以分為如下兩步：
【2】導(dǎo)入數(shù)據(jù)
use "data_for_DID.dta", clear
【3】設(shè)置面板
xtset FID time（FID為你的研究對象屠列，time為時間啦逆。例如有12個城市，2013~2019年笛洛，那么FID的個數(shù)就是12個夏志，time的個數(shù)就是7個）
（1）我們需要生成一個政策虛擬變量post_cmc，用以表示地區(qū)i在t期是否被開放為通商口岸苛让。我們只需比較樣本各期時間與開放時間（政策時點）即可沟蔑，如果是在政策時點后，則取值為1蝌诡，否則為0溉贿。
gen post_cmc=(year>cmcyear)
對于像嘉興府這樣未被開放為通商口岸的控制組個體枫吧，其開放時間是缺失的浦旱，Stata中將缺失值定義為無窮大，所以政策虛擬變量post_cmc會一直取值為0九杂。
（2）使用被解釋變量人口增長率AnnualGrowth對政策虛擬變量post_cmc進行回歸（加入個體固定效應(yīng)和時間固定效應(yīng)颁湖，用以更為精確地反映個體特征和時間特征）宣蠕，政策虛擬變量post_cmc反映的就是通商口岸的設(shè)置對近現(xiàn)代人口增長的影響。
DID模型與固定效應(yīng)模型有著千絲萬縷的關(guān)系甥捺，和之前一樣抢蚀，多期DID的Stata命令主要有三個，分別是reg命令镰禾、xtreg命令和reghdfe命令皿曲。reg命令使用的估計方法最小二乘虛擬變量方法（LSDV），通過在回歸方程中引入虛擬變量來代表不同的個體吴侦，可以起到和固定效應(yīng)組內(nèi)估計方法（FE）同樣的效果（已經(jīng)被證明）屋休，賈瑞雪老師使用的就是reg命令。
xtset id year
reg AnnualGrowth post_cmc $control_fe i.id i.year, vce(cluster id)
（id 就是上文導(dǎo)入的FID备韧，目的是研究對象進行聚類劫樟；i.year是上文導(dǎo)入的time，這里我的理解是控制時間织堂；i.id 我的理解是固定個體）
xtreg,fe是固定效應(yīng)模型的官方命令叠艳，使用這一命令估計出來的系數(shù)是最為純正的固定效應(yīng)估計量（組內(nèi)估計量），所以對于面板數(shù)據(jù)的DID模型易阳，我們使用更多的是xtreg,fe命令附较。
xtset id year
xtreg AnnualGrowth post_cmc $control_fe i.year,fe vce(cluster id)
第三個命令是reghdfe命令，也是一直以來我最推薦的固定效應(yīng)命令潦俺。reghdfe命令支持多維固定效應(yīng)翅睛，運算速度快，并且不會匯報一大長串虛擬變量回歸結(jié)果黑竞。
reghdfe AnnualGrowth post_cmc $control_fe, absorb(id year) vce(cluster id)
