Stata18异质性DID之Bacon分解新操作案例

来源：2023-06-01-Stata18异质性DID之Bacon分解新操作案例（四）

2023年4月25日，Stata如约而至发布了18.0新版本。里面带来了异质性DID稳健估计量等最新功能！为了让大家更好的了解Stata新版本更新的重要功能，计量经济学服务中心&数量经济学已于2023年4月26日推出了系列文章资源+视频讲解，这些专题资源基于Stata18.0软件基本操作讲解的，这估计是当时国内*早附带Stata18.0软件操作的文章了。

Stata18新功能一览表

---

• 2023年4月24日推文：重磅！倒计时！Stata18.0本周会发布么？？？
• 2023年4月25日推文：重磅！预言成真，Stata18.0新版本已经发布，异质性DID新官方命令发布了

---

关于Stata18.0中异质性DID里面更新的培根分解新操作，推文如下：

• 4月28日推文：4大Bacon分解操作技能，你都Get到了吗 (qq.com)
• 5月3日推文：Stata18新版本Bacon分解+bacondecomp新+旧版本操作及ddtiming命令（四合一）操作应用 (qq.com)

---

今天我们继续介绍关于Bacon分解的那些知识!

下面我们将讨论在estat bdecomp中实现的处理效果分解方法。

处理时间变化的设计的 DID 处理效果分解

我们使用一个（虚构的）数据集，其中单位是犬种，在美国养犬俱乐部的注册号是在 10 年内观察到的。我们希望了解注册号是否因广泛放映的电影中的品种而变化。以下是数据摘录：

use https://www.stata-press.com/data/r18/akc

我们列出了两个品种的数据，Old English Sheepdog and the Otterhound。我们可以看到，古英国牧羊犬Old English Sheepdog出现在 2037 年上映的电影中，因此是一个经过处理的单位。the Otterhound水獭犬品种没有这样的运气，属于从未接受过处理的群体。我们的结果变量注册量衡量美国养犬俱乐部的注册数量。我们使用 xtdidregress 来拟合双向固定效应 DID 模型：

xtset breed
xtdidregress (registered) (movie), group(breed) time(year) 

注意：ATET估计值根据面板效应和时间效应进行了调整。注意：处理发生在不同的时间。结果为：

我们可以看到，我们的整体 ATET  估计约为 2130 只，也就是说，如果该品种的成员作为电影主角，我们预计狗的注册数量将增加 2，130 只。我们现在可以使用 estat bdecomp 将此数字分解为其组成 2×2 部分：

estat bdecomp 

结果为：

. estat bdecomp

DID treatment-effect decomposition

ATET = 2129.655                                      Number of obs     = 1,410
                                                     Number of groups  =   141
                                                     Number of cohorts =     4

ATET decomposition summary                ATET component                Weight
------------------------------------------------------------------------------
Treated vs never treated                        2166.021              0.970810
Treated earlier vs later                       936.70729              0.013157
Treated later vs earlier                       906.60588              0.016033
------------------------------------------------------------------------------

Full ATET decomposition                  2x2 coefficient                Weight
------------------------------------------------------------------------------
Treated vs never treated
  2034 vs never treated                        1673.5688              0.157125
  2036 vs never treated                        1999.2515              0.140290
  2037 vs never treated                        2315.6702              0.673394

Treated earlier vs later
  2034 vs 2036                                 580.59722              0.001132
  2034 vs 2037                                 836.17778              0.008488
  2036 vs 2037                                 1291.9333              0.003537

Treated later vs earlier
  2036 vs 2034                                   542.825              0.001886
  2037 vs 2034                                 831.70833              0.011318
  2037 vs 2036                                 1448.7167              0.002829
------------------------------------------------------------------------------
Note: Number of cohorts includes never treated.
Note: The ATET reported by xtdidregress is a weighted average of the ATET components. If any
      component is substantially different from the ATET reported by xtdidregress and the
      weight is large, consider accounting for treatment-effect heterogeneity by using
      xthdidregress.

. 

注意：xtdidregress 报告的 ATET  是 ATET 组件的加权平均值。如果任何成分与 xtdidregress 报告的 ATET 有很大差异，并且权重很大，请考虑使用 xthdidregress 来考虑处理效果异质性。

在此示例中，处理次数相对较小。但是，如果该数字很大，则根据其权重绘制 2×2 系数可能会有所启发。我们通过使用estat bdecomp的graph选项来做到这一点。因为我们不需要再次看到数值结果，所以我们还指定了值得注意和无标题选项：

estat bdecomp， graph notable noheader

结果为：

蓝色圆圈是来自处理与从未处理的组件的 2×2 系数; 红色菱形是较早处理与较晚处理的组件的 2×2 系数; 绿色三角形是较晚处理与较早处理的 2×2 系数。三角形和菱形虽然与整体 ATET（由虚线参考线表示）有很大不同，但权重都非常接近于零。这表明这些系数对整体 ATET  的贡献可以忽略不计。