你想评估一项补贴政策。政策后,居民平均支出几乎没变。听起来像“没效果”。但也可能是另一回事:低支出者花得更少,高支出者花得更多,两边刚好在平均数上抵消了。
差分中的差分(Difference-in-Differences,DiD)常用来回答这类政策问题。它比较处理组与对照组在政策前后的变化,再用对照组扣除共同的时间趋势。关键前提叫“平行趋势”:如果没有政策,两组本来应该以相同方式变化。
问题在于,传统 DiD 主要看均值。若想知道政策怎样改变低位、中位和高位人群,就要研究整个结果分布。现有做法往往需要把平行趋势施加到整套分布变化上,这可能过强。
Kyle Schindl、Joseph Fry、Florian Gunsilius 和 Larry Wasserman 的论文《Parallel Trends as a Spectrum》提出:不要把平行趋势当成一道非真即假的判断题。研究者可以只对自己关心的分布特征作出假设,从均值一路走向完整分布,形成一条连续的限制谱。以下方法与结论均来自该论文;材料虽为全文节选,但仍只有这一个独立信源,论文也未提供可供本文引用的实证效果数字。
均值与完整分布之间,原来还有很多档
先看两个端点。
最弱的一端是均值平行趋势。它假设:没有政策时,处理组和对照组的平均结果会以相同幅度变化。这个条件可以识别 ATT——处理组的平均处理效应,却无法说明分布形状怎样改变。
最强的一端是分布层面的平行趋势。以论文讨论的加性分位数平行趋势为例,它要求两组在每个分位位置上的未处理结果都发生相同位移。分位数可以理解为把人按结果从低到高排列后所处的位置。这个条件能补出处理组在没有政策时的完整反事实分布,但也要求两组几乎整条分布都同步移动。
所谓“反事实”,就是现实中没有发生、却是因果比较所必需的另一种结果:接受政策的人,如果没接受政策,本来会怎样。它无法直接观察,只能依赖假设推断。
论文认为,这两个端点之间缺了一大片实用空间。它用 L-moments 填上这段空白。
像调音台一样选择分布特征
L-moments 是一组概括分布的统计量。它们由次序统计量的线性组合定义,也可以写成分位数函数的加权积分。第一阶对应均值,第二阶描述离散程度,第三阶描述不对称性,更高阶逐步刻画尾部和更细的形状。与普通高阶矩相比,L-moments 通常较少受极端值影响;只要结果存在有限的一阶绝对矩,各阶 L-moment 就有定义。
可以把一条分布想成一段复杂声音。完整分布是整段波形,L-moments 则像调音台上的不同旋钮:一个控制整体位置,一个控制铺得多开,另一个控制左右是否对称。研究者不必为了研究某个旋钮,就承诺整段波形都遵循同一规律。
论文把分位数函数投影到 shifted Legendre polynomials——一组彼此正交的多项式基——上。这样,复杂的分布变化就能拆成一串 L-moment 坐标。只约束第一阶,会回到传统的均值平行趋势;约束所有坐标,则等价于分布层面的平行趋势;选取有限或不连续的一组坐标,就得到中间档位。
关键在于“按坐标识别”。如果研究者只关心离散程度,就只需要相应坐标满足平行趋势,不必同时要求均值、偏度或其他未选特征也平行。这不是检验平行趋势究竟为真还是为假,而是先问清楚研究目标,再选择足以识别该目标的限制。
第二阶为何特别有意思
论文指出,第二阶 L-moment 与 Gini mean difference 成比例。Gini mean difference 是从总体中独立抽取两个人后,两者结果绝对差的期望,用来描述绝对离散程度。
按论文推导,第二阶坐标本身就能识别政策对这种两两差异的影响,而且不要求均值或更高阶特征满足平行趋势。摘要把它简写成“第二阶 L-moment identifies the Gini mean difference”,但两者的具体倍数取决于归一化方式,因此更稳妥的说法是“成比例并可据此识别”。
若同时识别第一、第二阶,并且结果非负、均值为正,还可以研究 Gini coefficient——把绝对离散程度除以平均水平后的相对不平等指标。论文特别区分了这两个概念:给所有人增加同样的正数,不改变两两之间的绝对差,却会因为均值上升而降低 Gini coefficient。
类似地,第二、第三阶可以组合成标准化的 L-skewness;第二、第四阶可以组合成 L-kurtosis,用来描述不对称和尾部形状。这些比率不受整体平移和正比例缩放影响,更接近纯粹的“形状变化”。
选得越多,看到得越完整
选择一组 L-moments 后,研究者可以重建分位数处理效应曲线的一部分。选取前 个坐标,相当于得到一条最高 次多项式对真实曲线的最佳投影;随着坐标增加,近似会更细。
但“近似更细”不等于“免费获得完整分布”。每增加一个坐标,也就增加一条需要相信的识别限制。有限阶重建甚至可能不是单调函数,因而未必能直接作为合法的分位数函数。论文讨论了用单调重排修正的方法,但这会改变部分高阶坐标。
作者还研究了高阶 L-moments 在什么条件下衰减,也就是后续坐标会多快变小。结论取决于分位数函数的矩条件、平滑程度和端点行为。仅仅拥有很多普通矩,并不自动保证很快的近似速度;即使分布有界,分位数函数的不连续也可能拖慢衰减。论文进一步给出有效估计量、非参数效率界,以及针对阶数随样本增长的一组坐标做同时推断的方法,并把框架扩展到重复横截面、面板数据和错期采用等 DiD 场景。
为什么这值得关注
这套思路改变的不是一道技术检验,而是提问顺序。过去,研究者容易先在“只看均值”和“识别完整分布”之间选边;现在可以先问:我真正关心的是平均水平、绝对离散、相对不平等,还是尾部形状?然后只为这个目标承担必要的假设。
论文还指出,非零的高阶坐标意味着两种边际结果分布之间不只是整体平移。它也能给个体处理效应的差异程度提供下界。不过,这仍不能识别个体效应的完整分布,因为那还取决于同一个人在处理与未处理两种状态下结果的联合关系。下界为正可以说明存在异质性;下界为零却不能证明人人受到的影响相同。
作者以 Medicare Part D 推行后的医疗支出分布变化作为应用案例。现有材料没有给出该案例的具体估计值,因此不能据此判断方法在现实数据中改善了多少。
局限与未知
- 这条“谱”没有消除平行趋势假设,而是把假设拆细。选择哪些坐标仍需领域知识,也无法用政策后的反事实直接验证。
- 加性分位数框架与 changes-in-changes(通过非线性映射补出反事实分布的方法)给出的是不同外推规则,并非同一假设的两种写法。论文把如何检验或组合二者留作未来方向。
- 论文建立了估计与推断理论,但现有供稿没有可引用的模拟或实证数值,也没有与既有方法的量化比较。因而“更弱条件”应理解为限制范围更窄,而不是已经得到外部验证的全面性能优势。