因子模型与合成控制法

作者:Butts,原文地址:https://www.kylebutts.com/blog/factor-model/
译者:许文立,CIMERS/IMIM,xuweny87@hotmail.com.
注:并不是完全按照原文翻译,有一些译者自己的理解和注释,并增加了一些相关的知识介绍。

当我第一次看到合成控制法时,我真的很难理解它什么时候有效,什么时候无效。当然,你需要很多处理前时期。但还需要什么呢?"较好的处理前拟合程度"。如果我在思考一个研究问题,我如何判断是否能事先获得较好å的处理前拟合程度呢?当然,通过观察结果来判断并不是一个好主意(请参考Roth(2022,AER:Insights))。

我是在一个完全不同的环境中偶然发现因子模型的。但是,在我学会了这个模型之后,我开始在有关合成控制的论文中一次又一次地看到它的出现。它们对我的理解产生了巨大的影响,因此我想与大家分享我的心得。

首先,我将与大家分享什么是因子模型,然后再回过头来谈谈它与合成控制的联系。作为一个重要的免责声明:合成控制估计量在各种数据生成过程中都是一致的。然而,因子模型是一个主要案例,在教学上非常有用。

一、因子模型

最简单的因子模型可以表示如下:
y_{it}=f_t' \gamma_i+u_{it}
其中,f_tp \times 1的不可观测“因子”向量,\gamma_ip \times 1的不可观测“因子载荷”向量。注意,这两者看起来像时间固定效应和个体固定效应,不过是以交互项的形式出现。

我们可以将因子f_t理解成宏观经济冲击(注:也就是对所有地区,所有行业,所有个体都产生影响的冲击),而\gamma_i是个体暴露于冲击的程度(注:例如,有些行业受到宏观经济冲击较小,而有些行业影响较大)。二者乘积描述的是个体如何受到宏观经济冲击的影响。例如,考虑y是GDP。因子载荷\gamma_i是一个地区制造业就业份额。而因子f_t可能是对所有制造业就业的国家层面的共同冲击(例如,对国家经济的进口竞争或其他变化,Aurto et al,2013)。两项的乘积表示国家层面冲击对一个地区GDP的效应。

另一种理解是\gamma_i表示时间不变的特征,而其对结果的效应f_t随时间变化。例如,y是工资,因子载荷\gamma_i是工人的教育,因子f_t可能是对教育回报的国家层面共同冲击。两者乘积比哦啊是国家冲击对一个工人工资的影响(注,最近的研究生扩招对教育-工资溢价的影响)。

上述例子仅仅只是因子模型的一些例子而已,它的真实能力其实在于我们不需要测度\gamma_if_t来控制两类因素。就像个体和时间固定效应一样,我们可以思考它们是重要的混淆因子因子,但我们实际上并不用测度它们。

注意,TWFE实际上是因子模型的一种特例。我们取\gamma_i=(1,\mu_i)f_t=(\lambda_t,1)。然后,模型变成y_{it}=\mu+i+\lambda_t+u_{it}。然后,因子模型允许个体有差异化的趋势(主要是因为不同个体暴露在共同冲击中的程度不同)。

因为因子模型也包含了TWFE,那么,我们为什么不在DID中仅仅使用因子模型?最大的问题在于,使用因子模型需要长面板(长时期)来估计因子和因子载荷。而我们使用DID的大部分情形都是较短的时期(短面板)。

二、合成控制法

合成控制法的基本思想:识别出能较好地近似处理前时期处理个体结果的控制组。合成控制估计量加权这些控制组。Ababie等指出,当结果模型是因子模型的时候,当处理前时期非常大,合成控制估计量的偏误会趋向于0(需要注意的是,在其它DGP下,合成控制估计量也是一致的)。

在因子模型下,处理个体的未处理潜在结果的合成控制估计量事控制组结果的加权平均\hat{y}_{it}^{sc}=\sum_j w_j y_{jt},其中,w_j是合成控制算法得到最优权重。我们将i=0作为处理个体。忽略掉误差项,我们可以写出下式:
\hat{y}_{it}^{sc}=\sum_j w_j y_{jt}=\sum_j w_j f_t'\gamma_j=f_t'\sum_j w_j \gamma_j=f_t'\sum_j \gamma^{sc}
最后一项显示:合成控制个体的结果受到相同因子的影响,且因子载荷是控制组个体因子载荷的加权平均值。

那么,合成控制组什么时候是有效的?从上式可以明确地看出,当合成控制个体的因子载荷看出来与处理个体的因子载荷尽可能的相近时(\gamma^{sc}=\gamma_0)。也就是说,处理组和挑选出的控制组受到国家层面的共同冲击相近时,合成控制法才比较适用。如果满足这个条件,我们预测的反事实结果才比较好。直觉上来说,合成控制个体和处理个体受到相同的因子冲击影响,且暴露在冲击中的程度相同才行。

我发现这种解释非常有用。作为一名研究人员,我可以事先开始考虑这些条件。首先,我会考虑在我的研究环境中很重要的因素(例如,在 GDP 的例子中,我认为一个重要的混杂冲击可能是产业随时间变化的冲击)。其次,我会考虑所处理的个体的因子载荷(例如其产业组合)是否可以很好地近似于备选控制组。在上述的例子中,如果我所处理的地区在其产业组合方面是一个异常值,那么合成控制法就会失效。

我不确定我是否在任何地方看到过这样的建议,但我认为计算重要特征的加权平均值是一个很好的检验方法,你认为这些重要特征可能会混淆处理选择(由基准协变量X替代),看看合成控制加权值是否接近处理个体加权,\sum w_j X_j \approx X_0

最后一个联系是关于长面板对估计合成控制的重要性。在上述工作中,请注意我忽略了误差项。这是因为,在长面板中,误差项为零。但在短面板中却不是这样,因为合成控制过程最终会与噪声匹配,而不是与因子载荷匹配,这使得 \gamma^{sc}\neq \gamma_0

三、合成控制论文中的因子模型

下面,我们来看看合成控制文献中,因子模型多么常见。

We consider the properties of the SC and related 
estimators, in a linear factor model setting, when the 
pre-treatment fit is imperfect. We show that, in this 
framework, the SC estimator is generally biased if 
treatment assignment is correlated with the 
unobserved heterogeneity, and that such bias does 
not converge to zero even when the number of pre-
treatment periods is large.

-3、Ben-Michael, Feller, and Rothstein(2021,JASA):"The Augmented Synthetic Control Method"

Considers how to use regularization on the weights 
to prevent over-fitting in shorter panels. This helps 
better approximate the true treated unit's factor-
loadings in a factor-model.
They impute untreated potential outcomes by 
estimating a low-rank factor model using nuclear-
norm regularization to prevent over-fitting.

四、短面板

上面提到的论文都依赖于长面板,才能恰当识别出因子结构。Butts与他的合作者研究了短面板中也保持合成控制法有效的估计量。

他们的研究表明,只要有多个处理个体,处理效应估计量仅仅要求因子f_t的一致估计,而不要求\gamma_i的估计。例如,IV(准长差分quasi-long differencing)估计出f_t的一致估计量,或者用相同因子结构(共同相关效应CCE,参见Kyle Butts and Nicholas Brown(2023))
所影响的协变量估计出f_t的一致估计量。这两种方法都可以得到事件研究图。参见 Kyle Butts and Nicholas Brown(2023):Generalized Imputation Estimators for Factor Models

image.png
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 204,189评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 85,577评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 150,857评论 0 337
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,703评论 1 276
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,705评论 5 366
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,620评论 1 281
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 37,995评论 3 396
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,656评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 40,898评论 1 298
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,639评论 2 321
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,720评论 1 330
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,395评论 4 319
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 38,982评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,953评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,195评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 44,907评论 2 349
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,472评论 2 342

推荐阅读更多精彩内容