did模型用什么软件

stata。因为学到的方法是用reg或xtreg命令与treated*t的交互项孙态来丢可以腊和用实现DID，因此did模型用stata软件。DID是一款为创造时尚和追轮凯盯求个性的消费者提供表达自身情感诉求的原创设计平台软件。

识别策略

清朝末期，清政府与西方列强签订了一系列不平等条约，开放沿江沿海等城市作为通商口岸即是不平等条约的主要内容之一。贾瑞雪老师（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)

可以用absdid命令。 absdid depvar [if] [in], tvar(varname) xvar(varlist) [options]

例毁含：

use absdid

absdid dlwage, tvar(union97) xvar(age black hispanic married grade)

absdid dlwage, tvar(union97) xvar(age black hispanic married grade) 宏雹  蔽余帆 sle

csinf(0.01) csup(0.99)

absdid dlwage, tvar(union97) xvar(age black hispanic married grade)     order(4)

等。

---