DID边际分析：让政策评价结果更加丰满

👇 连享会 · 推文导航 | www.lianxh.cn

• 🍎 Stata：Stata基础 | Stata绘图 | Stata程序 | Stata新命令
• 📘 论文：数据处理 | 结果输出 | 论文写作 | 数据分享
• 💹 计量：回归分析 | 交乘项-调节  | IV-GMM | 时间序列 | 面板数据  | 空间计量 | Probit-Logit | 分位数回归
• ⛳ 专题：SFA-DEA | 生存分析  | 爬虫 | 机器学习 | 文本分析
• 🔃 因果：DID | RDD | 因果推断 | 合成控制法 | PSM-Matching
• 🔨 工具：工具软件 | Markdown | Python-R-Stata
• 🎧 课程：公开课-直播 |  计量专题 | 关于连享会

作者：卢梅

目录

• 1. 前言

• 2. margins 命令介绍

• 3. DID 的边际分析

• 4. 对 margin 结果的修正

• 5. 文中使用的主要代码

温馨提示： 文中链接在微信中无法生效。请点击底部「阅读原文」。或直接长按/扫描如下二维码，直达原文：

1. 前言

本文在翻译下文的基础上进行了一些修改和补充：

Wheeler, Andrew. P., Some Stata notes – Difference-in-Difference models and postestimation commands, -Link-

倍差法 (DID) 是政策评价中被广泛应用的方法，DID 的结果表明某项政策的实施将会对因变量平均产生多少影响，但是很多时候我们都想知道，在某一个自变量处于某一特定水平，或者是某几个自变量处于特定水平时，政策效应到底有多大。

比如说，某项政策对于居民健康有正的影响，但是我们想知道这项政策对于 20 岁、 30 岁、 60 岁的居民影响有多大。将边际分析 (margin) 应用到 DID 中，可以观察到在任何一个想要的位置上政策冲击的所带来的边际效应是多少。margins 命令可以实现这一功能。接下来让我们近距离观察一下 margin。

2. margins 命令介绍

有关 margins 命令的详情，参见 Stata: 边际效应分析 和 Stata：因子变量全攻略 (微信版)。

在 Stata 命令窗口中输入 help margins 可以查看完整的帮助文件。事实上，它与 Stata 中的 predict，lincom 等命令有密切的关系。

• predict 的功能是用现有的估计结果 (主要是估计所得的系数) 来估计观测对象的统计指标。
• lincom、nlcom的功能在于估计已有系数的线性或非线性表达式。
• margin 命令综合了 predict 与 lincom、nlcom 的功能，来估计边际效果。margin 是指在基于数据集计算的拟合模型基础上，看在不同的协变量值的情况下，边际效应是多少。

简单来说，margin 其实就是，先利用已有的数据估计出模型系数，然后再用这些系数算出线性或非线性表达式的值，接着预测协变量处于一定水平上时的边际效果。

margin 本身就是一个可以实现众多功能的命令，关于它此处不做过多介绍，大家有兴趣的话可以 help 一下。

3. DID 的边际分析

Andrew P Wheeler 教授生成一份数据来具体讲述 margin 运用到 DID 中的示例，我们也使用 Andrew P Wheeler 生成的数据来给大家说明。

数据结构：Y 代表被解释变量，Post 和 Exper 可以理解为政策冲击的虚拟变量和是否为政策实施地区的虚拟变量，它们二者的交乘项系数表示政策效果，还有一个 Ord 序列变量就相当于我们经常用的 id。

很简单的数据，但是要注意，虽然 DID 及后续分析中用到了 i.month，但是这份数据并不是多期 DID，不同的月份只是代表不同的分类， i.month 表示在全年 12 个月中每个月都有一个特有的 DID 回归系数，可以将 month 类比为对于 20 多岁组， 30 多岁组， 40 多岁组这种，一月份组，二月份组，十二月份组。

DID 回归代码如下，之所以用 xtgee、poisson 广义估计泊松模型是由这份数据的性质决定的，这个我们可以先不用管它，就把它当作一般的回归就好。

回归结果为：

接下来对一月份组进行的边际分析，就是说此项政策实施对于一月份组来说效果有多大 (base 默认为一月) ，代码如下：

结果如下：

以图形展示 margin 结果，代码如下：

图形如下：

如果不想要两边边的 end cross hairs (还是不要让我翻译了，哈哈哈，我的第一反应是边边起横长的毛发) 可以输入如下代码，然后看一下效果，就知道差别在哪里了。

去掉 end cross hairs 的图形：

直观来看，对照组因变量上升了 1 个单位，实验组因变量下降了 2 个单位。

但是 Andrew 教授的同事提出，设置对照组的意义在于看到如果没有实行政策时会是怎样的情形，对照组因变量是下降的，那么实验组因变量也应该 是下降的，我们更应该看的是实验组因变量的相对变动。

我们需要借助lincom、nlcom等指令计算已估系数的表达式，然后还原出真正实验组对于对照组的相对变化。

4. 对 margin 结果的修正

lincom、nlcom之前，我们可以先用test检验一下相对下降是否存在，代码如下：

结果为：

结果显示，相对下降显著。当然也可以对各分组系数做联合检验，代码如下：

结果如下：

结果依然显著。接下来我们使用lincom、nlcom来看真实的相对下降是多少。

结果如下：

这一结果与在回归模型中，直接加入i.Exper i.Post的效果是一样的。输入如下代码：

结果如下：

结果与通过 lincom 命令得到的交乘项系数是一样的，因为模型是泊松模型，所以我们再使用 nlcom 还原到更直观的数值变动。代码如下：

结果如下：

因为看的是第一个月组的情况，month 前的系数为 0，所以nlcom命令常数项后面不用加 month 的系数，如果想看 2 月份组、3 月份组的政策效应的话，就应该把 month 的系数加上。

政策实施后因变量是 11.79 接近 12 的值，再想想之前那个长着 end cross hairs 的图，实验组从 8 变到 6 ，其实是不准确的，应该是从 12 变到 6，政策效应是使因变量减少了 8 个单位每月！

所以 DID 时，其实应该使用 xtgee Y i.Exper i.Post Exper#Post i.Month, family(poisson) corr(ar1) 这种形式的回归，得到回归结果后，再用 nlcom 计算出实际变动值。

5. 文中使用的主要代码

New！ Stata 搜索神器：lianxh 和 songbl  GIF 动图介绍
搜： 推文、数据分享、期刊论文、重现代码 ……
👉 安装：
  . ssc install lianxh
  . ssc install songbl
👉  使用：
  . lianxh DID 倍分法
  . songbl all

🍏 关于我们

• 连享会 ( www.lianxh.cn，推文列表) 由中山大学连玉君老师团队创办，定期分享实证分析经验。
• 直通车： 👉【百度一下：连享会】即可直达连享会主页。亦可进一步添加 「知乎」,「b 站」,「面板数据」,「公开课」 等关键词细化搜索。