SCM合成控制法.docx
- 文档编号:6869424
- 上传时间:2023-01-11
- 格式:DOCX
- 页数:15
- 大小:847.65KB
SCM合成控制法.docx
《SCM合成控制法.docx》由会员分享,可在线阅读,更多相关《SCM合成控制法.docx(15页珍藏版)》请在冰豆网上搜索。
SCM合成控制法
合成控制法
(一)
经济学家为何热衷反事实
经济学家常要评估某政策或事件的效应。
此政策可能实施于某国家或地区(省、州或城市)。
最简单(天真)的方法是考察政策实施前后的时间序列,看所关心的结果(outcomeofinterest)如何变化。
但此结果还可能受其原有变化趋势的影响,或其他同时发生的混淆性事件(confounder)的作用。
为此,常使用“鲁宾的反事实框架”(Rubin'scounterfactualframework),即假想该地区如未受政策干预将会怎样,并与事实上受到干预的实际数据进行对比,二者之差即为“处理效应”(treatmenteffect,借用医学术语)。
困难之处在于,我们无法观测到“该地区如未受政策干预将会怎样”(反事实)。
选择控制组是门艺术
常用解决方法是,寻找适当的控制组(controlgroup),即在各方面都与受干预地区相似却未受干预的其他地区,以作为处理组(treatedgroup,即受到干预的地区)的反事实替身(counterfactuals)。
但通常不易找到最理想的控制地区(controlregion),在各方面都接近于处理地区(treatedregion)。
比如,要考察仅在北京实施的某政策效果,自然会想到以上海作为控制地区;但上海毕竟与北京不完全相同。
或可用其他一线城市(上海、广州、深圳)构成北京的控制组,比较上海、广州、深圳与北京在政策实施前后的差别,此方法也称“比较案例研究”(comparativecasestudies)。
但如何选择控制组通常存在主观随意性(ambiguity),而上海、广州、深圳与北京的相似度也不尽相同。
为此,AbadieandGardeazabal(2003)提出“合成控制法”(SyntheticControlMethod)。
其基本思想是,虽然无法找到北京的最佳控制地区,但通常可对中国的若干大城市进行适当的线性组合,以构造一个更为优秀的“合成控制地区”(syntheticcontrolregion),并将“真实北京”与“合成北京”进行对比,故名“合成控制法”。
合成控制法的一大优势是,可以根据数据(data-driven)来选择线性组合的最优权重,避免了研究者主观选择控制组的随意性。
西班牙恐怖活动引发的计量方法
AbadieandGardeazabal(2003)的初衷是以合成控制法研究西班牙巴斯克地区(Basquecountry)恐怖活动的经济成本。
MIT经济系教授AlbertoAbadie(此前长期任教于哈佛大学肯尼迪学院),正是来自于巴斯克地区,一个毗邻法国的西班牙自治地区。
巴斯克人长期居住于巴斯克地区,拥有独特的语言与文化,在历史上多次成功对抗强敌入侵。
在1970年代初,巴斯克地区的人均GDP在西班牙17个地区中排第三。
之后,由于民族独立的诉求未获满足,从1975年开始,巴斯克地区陷入有组织的恐怖活动之中。
恐怖活动重创巴斯克经济,至1990年代末,巴斯克地区的人均GDP在西班牙排名降为第六。
然而,70年代末至80年代初,西班牙整体经济也下行,故不易区分恐怖活动的单独效应。
而且,巴斯克地区在恐怖活动之前的经济增长潜力显然与西班牙其他地区也不尽相同。
为此,AbadieandGardeazabal(2003)使用西班牙其他地区的线性组合来构造合成的控制地区,并使得合成控制地区的经济特征与60年代末恐怖活动爆发前的巴斯克地区尽可能相似,然后把此后“合成巴斯克地区”(syntheticBasquescountry)的人均GDP演化与“真实巴斯克地区”(actualBasquecountry)进行对比。
如何构造合成控制
具体而言,假设共有(1+J)个地区,其中第1个地区为受到恐怖活动冲击的巴斯克地区,而其余J个西班牙地区未受冲击(在此J=16),构成潜在的控制组,称为“donorpool”(原意为“器官捐献库”,再次借用医学术语)。
一个潜在假定是,恐怖活动仅影响巴斯克地区,而未波及西班牙的其他地区(事实上恐怖活动也主要集中于巴斯克地区)。
将合成控制地区的权重记为以下J维列向量:
其中,w2表示第2个地区在合成巴斯克地区所占的权重,以此类推;所有权重皆非负,且权重之和为1。
w的不同取值即构成不同的合成控制地区,简称“合成控制”(syntheticcontrol)。
在此研究中,被解释变量为人均GDP,记为y。
影响y的解释变量或预测变量(predictors)包括投资率、人口密度、产业结构、人力资本等,详见下表。
在巴斯克地区爆发恐怖活动之前,记其各预测变量的平均值为向量x1(K×1维列向量,下标1表示“treatedregion”),即上表第
(1)列的数值(除了人均GDP)。
将西班牙其他地区相应预测变量的平均值记为矩阵X0(K×J维矩阵,下标0表示“controlregion”),其中第j列为第j个地区的相应取值。
显然,我们希望选择权重w,使得X0w尽可能地接近于x1,即经过加权之后,合成控制地区的经济特征应尽量接近处理地区。
为度量此距离,可使用二次型(类似于欧几里得空间中两点之间的距离)。
由于x1中的每个预测变量对于y的预测能力有大小之别,应在距离函数中享有不同的权重,故考虑以下有约束的最小化问题:
其中,V为(K×K)维对角矩阵,其对角线元素均为非负权重,反映相应的预测变量对于人均GDP的相对重要性。
此最小化问题的目标函数是二次函数,为“二次规划”(quadraticprogramming)问题,一般进行数值求解。
记此约束最小化问题的最优解为w*(V);显然,它依赖于对角矩阵V。
进一步,选择最优的V,使得在恐怖活动全面爆发之前,合成巴斯克地区的人均GDP与真实巴斯克地区尽量接近。
具体而言,记z1为(10×1)维列向量,包含巴斯克地区在1960-1969年间的人均GDP;记Z0为(10×J)维矩阵,其中每列为相应控制地区在1960-1969年间的人均GDP。
用Z0w*(V)来预测z1,然后选择V,以最小化“均方预测误差”(MeanSquaredPredictionError,简记MSPE),即将每期的预测误差平方后再求各期的平均:
求解此最小化问题,可得构成合成巴斯克地区的最优权重,w*=w*(V*)。
经过计算,AbadieandGardeazabal(2003)发现,只有两个地区的权重为正,即加泰罗尼亚(Catalonia,权重0.8508)与马德里(Madrid,权重0.1492),而其他地区的权重均为0。
直观上,Catalonia与Madrid的经济特征也与巴斯克地区最为相似。
合成控制法的“效果图
得到合成巴斯克地区的权重之后,即可计算其人均GDP在样本期间的演化过程。
记巴斯克地区在样本期间(假设为T期)的人均GDP为向量y1(T×1维列向量)。
记其他地区在样本期间的人均GDP为矩阵Y0(T×J维矩阵),其中每列为相应地区的人均GDP。
由此可得合成巴斯克地区的人均GDP序列y1*=Y0w*。
最直观的方法是将y1与合成控制的y1*画时间趋势图,参见下图。
从上图可知,在1975年大规模恐怖活动爆发之前,真实巴斯克(实线)与合成巴斯克(虚线)的人均GDP十分接近。
二者在1975年后即开始分岔;而在1980与1990年代,真实巴斯克的人均GDP比合成巴斯克低约10%。
换言之,巴斯克恐怖活动的经济成本是损失了约10%的人均GDP
反事实的分析框架
假设共有(1+J)个地区,其中第1个地区受到政策干预(如有多个地区受到干预,可合并为一个大地区;或分别进行估计),而其余J个地区未受冲击(构成donorpool)。
记yit为地区i在第t期实际观测到的结果变量,其中i=1,...,J+1,而t=1,...,T。
记yitN为地区i在第t期如果未受政策干预的结果变量(上标N表示未受干预)。
记T0为政策干预开始之前(preintervention)的时期数,且1≤T0 记yitI为地区i在第t期的结果变量(上标I表示Intervention),如果地区i在第(T0+1)至第T期持续地受到政策干预。 假设政策在前T0期对于结果变量没有影响,即对于所有i与t≤T0,都有yit=yitN=yitI。 如果政策在实施之前即产生影响(比如,通过预期效应),则可重新定义T0为政策实际开始产生影响之前的那个时期。 一个潜在假定是各地区之间不会互相影响(nointerferencebetweenunits);特别地,控制地区的结果变量不受处理地区政策冲击的影响。 我们关心当i=1而t>T0时的处理效应: 在上式中,只要估计y1tN即可。 引入因子模型 假设yitN由以下“因子模型”(factormodel)所决定: 其中,上式右边第 (1)项δt为时间固定效应(timefixedeffects)。 第 (2)项的zi为可观测的向量(不受政策干预影响,也不随时间而变;比如,干预之前的预测变量之平均值)。 zi对于yitN的作用随时间而变,故zi的系数θt(未知参数)带时间下标t。 第(3)项为不可观测的“互动固定效应”(InteractiveFixedEffects),即个体固定效应ui与时间固定效应λt的乘积(Bai,2009)。 第(4)项ԑit为随机扰动项。 根据“因子分析”(factoranalysis)的术语,称第(3)项中不可观测的λt为“共同因子”(commonfactors),可理解为不同地区所面临的共同冲击(commonshocks),比如它有两个分量,分别表示技术冲击(technologicalshocks)与金融危机(financialcrises);而各地区对于共同冲击λt的反应并不相同,以ui来表示,称为“因子载荷”(factorloading)。 如果λt是一维且为常数,则上式简化为“双向固定效应模型”(two-wayfixedeffectsmodel),包含个体固定效应ui与时间固定效应δt。 由此可知,上式是双向固定效应模型的推广,它允许不同个体对于共同冲击的异质性反应(heterogeneousimpacts)。 合成控制的误差 记构造合成控制的权重向量为: 其中,w2表示第2个地区在合成控制所占的权重,以此类推;所有权重皆为非负,且权重之和为1。 对于任意给定的w,可将合成控制地区的结果变量写为: 将y1tN减去上式可得: 显然,如果能找到w,使得上式右边的 (1)式与 (2)式均为0,则上式的期望为0,故合成控制为y1tN的无偏估计。 但 (2)式中的ui不可观测,故不可行。 Abadieetal.(2010)证明,如果能找到w使得 即根据可观测的经济特征与干预前结果变量所选择的合成控制w,也会使得合成控制的不可观测特征接近于处理地区。 反之,如果无法找到w,使得合成控制能很好地复制(reproduce)处理地区的经济特征以及干预之前的结果变量,则不建议使用合成控制法。 合成控制法的基本性质 如果合成控制w*能很好地复制处理地区的经济特征与干预前的结果变量(关于如何计算w*,参见首篇《合成控制法 (一)》),则可定义如下合成控制估计量(SyntheticControlEstimator): Abadieetal.(2010)证明,在一定的正则条件下,如果合成控制w*能完全地复制(perfectlyreproduce)处理地区的经济特征与干预前的结果变量,则当干预前期数T0趋向无穷大时,合成控制估计量是渐近无偏(asymptoticallyunbiased)。 反之,如果合成控制w*只是不完全地复制(imperfectlyreproduce)处理地区的经济特征与干预前的结果变量(更为常见的情形),则合成控制估计量将是渐近有偏的。 此偏差被称为“内插偏差”(interpolationbias),因为使用内插法构造的合成控制并不是处理地区的完美的反事实替身。 因此,在使用合成控制法时,一个重要的模型设定检验就是,考察合成控制的经济特征与干预前的结果变量是否与处理地区足够接近。 合成控制法的Stata命令 Abadieetal.(2010)还提供了合成控制法的Stata程序synth,使得估计十分方便。 打开Stata后,可输入以下命令: .sscinstallsynth,replace (下载并安装synth程序) 其中,选择项“replace”表示如有此命令更新版本,可以新命令覆盖旧命令。 命令synth的基本句型为: synthyx1x2x3,trunit(#)trperiod(#)counit(numlist)xperiod(numlist)mspeperiod()figureresultsperiod()nestedalloptkeep(filename) 其中,“y”为结果变量(outcomevariable),“x1x2x3”为预测变量(predictors)。 必选项“trunit(#)”用于指定处理地区(trunit表示treatedunit)。 必选项“trperiod(#)”用于指定政策干预开始的时期(trperiod表示treatedperiod)。 选择项“counit(numlist)”用于指定潜在的控制地区(即donorpool,其中counit表示controlunits),默认为数据集中的除处理地区以外的所有地区。 选择项“xperiod(numlist)”用于指定将预测变量(predictors)进行平均的期间,默认为政策干预开始之前的所有时期(theentirepre-interventionperiod)。 选择项“mspeperiod()”用于指定最小化均方预测误差(MSPE)的时期,默认为政策干预开始之前的所有时期。 选择项“figure”表示将处理地区与合成控制的结果变量画时间趋势图,而选择项“resultsperiod()”用于指定此图的时间范围(默认为整个样本期间)。 选择项“nested”表示使用嵌套的数值方法寻找最优的合成控制(推荐使用此选项),这比默认方法更费时间,但可能更精确。 在使用选择项“nested”时,如果再加上选择项“allopt”(即“nestedallopt”),则比单独使用“nested”还要费时间,但精确度可能更高。 选择项“keep(filename)”将估计结果(比如,合成控制的权重、结果变量)存为另一Stata数据集(filename.dta),以便进行后续计算。 更多选择项,详见helpsynth。 加州控烟法有效吗? Abadieetal.(2010)随之将合成控制法应用于研究美国加州1988年第99号控烟法(Proposition99)的效果。 1988年11月美国加州通过了当代美国最大规模的控烟法(anti-tobaccolegislation),并于1989年1月开始生效。 该法将加州的香烟消费税(cigaretteexcisetax)提高了每包25美分,将所得收入专项用于控烟的教育与媒体宣传,并引发了一系列关于室内清洁空气的地方立法(localcleanindoor-airordinances),比如在餐馆、封闭工作场所等禁烟。 面板数据集 Abadieetal.(2010)使用的数据为美国1970-2000年的州际面板数据。 由于Massachusetts,Arizona,Oregon与Florida的州政府在1989-2000期间也引入了正式的控烟立法,无法作为控制地区,故将此四州从donorpool中去掉。 外,还去掉了在1989-2000期间将香烟消费税提高了每包50美分及以上的州(Alaska,Hawaii,Maryland,Michigan,NewJersey,NewYork,Washington),以及首都所在地华盛顿特区。 最后剩下38个州作为潜在的控制地区(即使加上去掉的州,结果也依然稳健)。 变量介绍 该研究的结果变量为cigsale(人均香烟消费量,包/年),预测变量包括retprice(平均香烟零售价格)、lnincome(人均收入对数)、age15to24(15-24岁人口所占总人口比重,年轻人为吸烟主力)、beer(人均啤酒消费量,烟酒不分家)。 这些预测变量均为1980-1988年的州平均值,另外还使用1975、1980与1988年的人均香烟消费量作为三个额外的预测变量。 另外,面板变量为state(州),而时间变量为year(年)。 Stata操作 Stata程序synth自带了加州控烟法的数据集smoking.dta。 安装命令synth后,输入以下命令即可打开此数据集,设为面板数据,然后进行合成控制法的估计: .sysusesmoking(打开数据集) .xtsetstateyear(设为面板数据) .synthcigsaleretpricelnincomeage15to24beercigsale(1975)cigsale(1980)cigsale(1988),trunit(3)trperiod(1989)xperiod(1980 (1)1988)figurenestedkeep(smoking_synth) 其中,“cigsale(1975)cigsale(1980)cigsale(1988)”分别表示人均香烟消费在1975、1980与1988年的取值。 必选项“trunit(3)”表示第3个州(即加州)为处理地区;必选项“trperiod(1989)”表示控烟法在1989年开始实施。 选择项“xperiod(1980 (1)1988)”表示将预测变量在1980-1988年期间进行平均,其中“1980 (1)1988”表示始于1980年,以1年为间隔,而止于1988年;其效果等价于“198019811982198319841985198619871988”,而前者的写法显然更为简洁。 选择项“keep(smoking_synth)”将估计结果存为Stata数据集smoking_synth.dta(放在Stata的当前工作路径)。 估计结果如下。 上表显示,大多数州的权重为0,而只有以下五个州的权重为正,即Colorado(0.161),Connecticut(0.068),Montana(0.201),Nevada(0.235)与Utah(0.335),此结果与Abadieetal.(2010)汇报的结果非常接近(细微差别或由于计算误差)。 考察加州与合成加州的预测变量是否接近: 从上表可知,加州与合成加州的预测变量均十分接近,故合成加州可以很好地复制加州的经济特征。 然后比较二者的人均香烟消费量在1989年前后的表现: 从上图可知,在1989年控烟法之前,合成加州的人均香烟消费与真实加州几乎如影相随,表明合成加州可以很好地作为加州如未控烟的反事实替身。 在控烟法实施之后,加州与合成加州的人均香烟消费量即开始分岔,而且此效应越来越大。 更直观地,可打开另一Stata程序,调用已存的数据集smoking_synth.dta,计算加州与合成加州人均香烟消费之差(即处理效应),然后画图。 .usesmoking_synth.dta,clear (如不打开另一Stata程序,则此数据集将覆盖原有的数据集smoking.dta) .geneffect=_Y_treated-_Y_synthetic (定义处理效应为变量effect,其中“_Y_treated”与“_Y_synthetic”分别表示处理地区与合成控制的结果变量) .labelvariable_time"year" .labelvariableeffect"gapinper-capitacigarettesales(inpacks)" (为了画图更漂亮,加上时间变量与处理效应的标签,可使用变量管理器(variablemanager)来方便地加标签) .lineeffect_time,xline(1989,lp(dash))yline(0,lp(dash)) (画处理效应的时间趋势图,并在横轴1989年处与纵轴0处分别画虚线,结果见下图) 上图显示,加州控烟法对于人均香烟消费量有很大的负效应,而且此效应随着时间推移而变大。 具体来说,在1989-2000年期间,加州的人均年香烟消费减少了20多包,大约下降了25%之多,故其经济效应十分显著(economicallysignificant)。 在使用合成控制法时,如何进行稳健性检验与统计推断? 合成控制法与回归法相比有何优缺点? 使用合成控制法有何注意事项? 敬请期待本介绍的完结篇——合成控制法(三)。 合成控制法的稳健性检验 为了检验上述合成控制估计结果的稳健性(参见《合成控制法 (二)》),Abadieetal.(2010)加入了更多的预测变量,比如失业率、收入不平等、贫困率、福利转移、犯罪率、毒品相关的逮捕率、香烟税、人口密度等;发现结果依然稳健。 另外一个担心是,地区之间无互相影响(nointerferencebetweenunits)的假定可能不满足,比如加州的反烟运动可能波及其他州,烟草行业或将其他州的香烟广告预算投入到加州,甚至从其他州走私便宜香烟到加州。 Abadieetal.(2010)根据史实对此进行了探讨,认为这些效应均不大,至少不可能导致上文图中如此大的处理效应。 安慰剂检验 上述结果为对控烟法处理效应的点估计。 此点估计是否在统计上显著(statisticallysignificant)? Abadieetal.(2010)认为,在比较案例研究中,由于潜在的控制地区数目通常并不多,故不适合使用大样本理论进行统计推断。 为此,Abadieetal.(2010)提出使用“安慰剂检验”(placebotest)来进行统计检验,这种方法类似于统计学中的“排列检验”(permutationtest),适用于任何样本容量。 “安慰剂”(placebo)一词来自医学上的随机实验,比如要检验某种新药的疗效。 此时,可将参加实验的人群随机分为两组,其中一组为实验组,服用真药;而另一组为控制组,服用安慰剂(比如,无用的糖丸),并且不让参与者知道自己服用的究竟是真药还是安慰剂,以避免由于主观心理作用而影响实验效果,称为“安慰剂效应”(placeboeffect)。 安慰剂检验借用了安慰剂的思想。 具体到加州控烟法的案例,我们想知道,使用上述合成控制法所估计的控烟效应,是否完全由偶然因素所驱动? 换言之,如果从donorpool随机抽取一个州(而不是加州)进行合成控制估计,能否得到类似的效应? 为此,Abadieetal.(2010)进行了一系列的安慰剂检验
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- SCM 合成 控制
![提示](https://static.bdocx.com/images/bang_tan.gif)