清朝末期，清政府與西方列強簽訂了壹系列不平等條約，開放沿江沿海等城市作為通商口岸即是不平等條約的主要內容之壹。賈瑞雪老師（2014）將近代通商口岸的設置作為壹項準自然實驗，評估了通商口岸對中國近現代人口和經濟發展的長期影響。
從1840年至1910年，中國壹***被迫開放了40多個通商口岸。不同通商口岸開放的時間（政策時點）是不同的，例如，廣州、福州、廈門、寧波和上海是在1842年《南京條約》後開放的，漢口、九江、南京、鎮江等城市是在1858年《天津條約》後開放的，而蘇州、杭州等城市是在1895年《馬關條約》後開放的......
如上表所示，蘇州府是在1896年被開放為通商口岸的，所以蘇州府的政策分期變量period在1896年之後取值為1，之前取值為0；而同處蘇杭地區的嘉興府則壹直沒有被開放為通商口岸，所以嘉興府的政策分組變量treat和政策分期變量period壹直取值為0。
事實上，我們可以發現交互項treat×period的取值和政策分期變量period的取值是壹毛壹樣的，所以 在多期DID中，我們其實就沒有必要去生成什麽交互項，只需用壹個政策虛擬變量予以替代就可以了 ，用以表示地區i在t期是否實施政策。當然，如果為了便於理解的話，可以嘗試去生成交互項，結果都壹樣滴！
Stata操作
1設置環境
cd?D:\experiments（位置壹般為數據所在文件夾）
多期雙重差分法（DID）的Stata操作可以分為如下兩步：
2導入數據
use "data_for_DID.dta", clear
3設置面板
xtset FID time（FID為妳的研究對象，time為時間。例如有12個城市，2013~2019年，那麽FID的個數就是12個，time的個數就是7個）
（1）我們需要生成壹個政策虛擬變量post_cmc，用以表示地區i在t期是否被開放為通商口岸。我們只需比較樣本各期時間與開放時間（政策時點）即可，如果是在政策時點後，則取值為1，否則為0。
gen post_cmc=(year>cmcyear)
對於像嘉興府這樣未被開放為通商口岸的控制組個體，其開放時間是缺失的，Stata中將缺失值定義為無窮大，所以政策虛擬變量post_cmc會壹直取值為0。
（2）使用被解釋變量人口增長率AnnualGrowth對政策虛擬變量post_cmc進行回歸（加入個體固定效應和時間固定效應，用以更為精確地反映個體特征和時間特征），政策虛擬變量post_cmc反映的就是通商口岸的設置對近現代人口增長的影響。
DID模型與固定效應模型有著千絲萬縷的關系，和之前壹樣，多期DID的Stata命令主要有三個，分別是reg命令、xtreg命令和reghdfe命令。reg命令使用的估計方法最小二乘虛擬變量方法（LSDV），通過在回歸方程中引入虛擬變量來代表不同的個體，可以起到和固定效應組內估計方法（FE）同樣的效果（已經被證明），賈瑞雪老師使用的就是reg命令。
xtset id year
reg AnnualGrowth post_cmc $control_fe i.id i.year, vce(cluster id)
（id 就是上文導入的FID，目的是研究對象進行聚類；i.year是上文導入的time，這裏我的理解是控制時間；i.id 我的理解是固定個體）
xtreg,fe是固定效應模型的官方命令，使用這壹命令估計出來的系數是最為純正的固定效應估計量（組內估計量），所以對於面板數據的DID模型，我們使用更多的是xtreg,fe命令。
xtset id year
xtreg AnnualGrowth post_cmc $control_fe i.year,fe vce(cluster id)
第三個命令是reghdfe命令，也是壹直以來我最推薦的固定效應命令。reghdfe命令支持多維固定效應，運算速度快，並且不會匯報壹大長串虛擬變量回歸結果。
reghdfe AnnualGrowth post_cmc $control_fe, absorb(id year) vce(cluster id)