上图为数据结构，其中year 为年份；region代表地区；t（若实施此政策为1）和treated（若处在该区域为1）为虚拟变量；jyl为被解释变量；lnrgdp、edu、gz、ins 为控制变量；gd 为t*treated 的交互项。
（1）双重差分模型：
Stata命令：reg jyl gd treated t edu ins lnrgdp gz,r
（2）加入个体和时间固定效应的模型
Stata命令：reg jyl gd edu ins lnrgdp gz i.year i.region,r
问题一：请问老师，加入个体和时间固定效应后的模型命令正确吗？如果再在此基础上加入个体与时间交乘项该如何操作命令？
（3）双重差分倾向得分匹配
Stata命令：
diff jyl, t(treated) p(t) cov(edu ins lnrgdp gz) kernel id(region) support logit report robust
diff jyl, t(treated) p(t) kernel id(region) logit cov(edu ins lnrgdp gz) support test
回复一
上述两个模型的命令都是正确的。在此基础上，无法加入个体（region）与时间（year）的彻底的交互项，因为这样的交互项们将与双重差分的关键变量gd
完全共线性。但如果你想反映个体的效应可能随时间变化，可以尝试加入个体的参数形态时间趋势，最简单的是加入个体的线性时间趋势。
命令如下：
reg jyl gd edu ins lnrgdp gz i.year i.region i.region#c.year, r