1.背景
一般来讲美国民主党更倾向于更多的联邦支出，案例研究是否民主党获选对于联邦支出的影响。通常情况下如果得票率大于50%即会获选，反之小于50%则会落选。因而得票率则为驱动变量X，此处0.5则可作为断点cutoff值（研究中为了方便使用，因而将得票率 – 0.5）作为驱动变量，即最终cutoff值为0，大于0则应该获选，小于0则应该落选）。结果变量Y为联邦支出。而且还包括另外2个控制变量。除此之外，还包括另外一个变量‘是否获选’作为判断是否模糊断点。本案例数据使用Stata软件的votex.sta数据，各数据的定义如下：
项 | 变量简称 | 变量说明 |
结果变量Y | y | Log fed expenditure in district |
处理变量或驱动变量X | x | Dem vote share minus .5 |
模糊变量fuzzy项 | fuzzy | Dem Won Race |
控制变量 | control1 | Incumbent |
控制变量 | control2 | Voting Age Population Share |
2 理论
RDD断点回归的分析知识点相对较多，从分析步骤包括，具体可精确断点或模糊断点的选择，模型选择，模型基本假定分析，模型分析，模型稳健性检验等。具体分为以下5步。
第一步、精确断点和模糊断点判断
判断精确断点或模糊断点的思路在于处理变量X被cutoff分为左右两侧后，是否真正决定‘实验走向’，比如本案例中cutoff值分成两组后即认为‘民主党是否当选’（命名为new_x），而fuzzy模糊项即真实是否当选项，如果与new_x与fuzzy项没有特别明显的不一致，甚至完全一样，此时则应该使用精确断点。反之如果new_x与fuzzy项有着明显的差异，此时使用模糊断点较为适合。
第二步、模型选择，通常指模型阶数的判断
研究X对于Y的影响时，二者的关系是线性关系（一阶），还是曲线二阶关系，也或者三阶关系。可首先通过图示直观查看，并且得出结论。待定模型阶数后，后续分析基于该阶数进行分析使用。至于‘带宽值’或者‘核函数’，通常默认即可，SPSSAU会自动找出最优带宽值，默认使用triangular三角核函数。
第三步、模型基本假定分析
RDD模型通常包括着一定的假设，通常包括‘断点适用性检验’和‘局部平滑性检验’。如下所述：
假定 | 说明 | SPSSAU是否提供 |
断点适用性检验 | 驱动变量不应存在人为操控，正常情况下断点cutoff值附近样本量基本一致才会具有随机性 | 通过图示直接查看断点左右两侧样本量密集性情况，如果出现一侧明显更多样本点则说明不适合 |
局部平滑性检验 | 如果有控制变量，那么控制变量在断点处不应该存在明显的跳跃现象 | 将控制变量看成‘驱动变量X’进行分析，图示查看即可 |
第四步、模型分析
在上述确认好精确或模糊断点，并且确定好模型阶数，并且模型适合时，则开始分析X对于Y的影响关系情况。
第五步、模型稳健性检验
模型分析后，还需要对模型稳健性进行检验。模型稳健性检验有多种方式，包括更换核函数法、更换断点值法，是否加入控制变量法，更换带宽值法，更换阶数，改变样本选择法，如下表格所述。
稳健性检验方式 | 说明 | SPSSAU是否提供 |
更换核函数 | 更换不同核函数时结论是否稳定 | 研究者更换核函数 |
更换断点 | 更改不同断点值时，结论是否稳定 | 研究者更换断点 |
是否加入控制变量 | 是否加入控制变量时，结论是否稳定（断点安慰剂检验Placebo Tests） | 研究者更换模型，即是否加入控制变量 |
换带宽 | 更换带宽值时，结论是否稳定 | 自动提供，且提供coefplot图 |
换阶数 | 更换阶数时，结论是否稳定 | 研究者更换阶数 |
样本选择 | 更换筛选样本时，结论是否稳定 | 研究者自行处理 |
3 .操作
第一步、精确断点和模糊断点判断
将X按cutoff值0分为两组，并且与fuzzy项进行交叉卡方分析。操作截图分别如下：
使用SPSSAU数据处理->数据编码功能，并且选择‘范围编码’，将x按cutoff值0分为两组（可通过描述分析得到x的最小值为-0.276，最大值0.470）。系统会自动生成一个新标题‘New_x’，将该项与fuzzy项进行交叉卡方，得到如下结果：
New_x代表x分为两组后的新变量，0代表cutoff值左侧（落选），1代表cutoff值右侧（当选）。而fuzzy项里面的0和1代表真实情况下‘是否当选’（0为落选，1为当选）。从上表格可以看到：二者数据完全一致，按cutoff值得到的131个‘落选’样本真实情况下也是‘落选’，按cutoff值得到的218个‘当选’样本真实情况下也是‘当选’。即意味着应该使用精确断点。
提示：
实际研究中，如果数据的gap较小，此时也可直接使用精确断点回归。
确认好为精确断点模型之后，接着进行第二步。
第二步、模型阶数判断
模型阶数判断时使用直观图示法。即首先进行模型分析，通过图示查看模型应该是一阶、二阶或三阶更加适合。首先操作如下图：
首先放入结果变量y，驱动变量x，2个控制变量。以及设置好断点值为0（默认不设置即为0），选中‘绘图’复选框。至于另外4个参数（带宽值、核函数、阶数和稳健性检验）默认即可。此步骤主要查看绘图，用于确认‘阶数’。得到图形如下：
从上图可以看到，最左侧‘线性拟合’即一阶时或者中间‘二次型拟合’即二阶时，模型拟合相对较好。可能‘二次型拟合’相对更适合。因而确定模型为二阶。并且后续以二阶为准进行分析。另外从上面三个图可以看到，断点值左右两侧附近的样本量基本均匀，即说明断点值选择适合没有受到人为操纵。
第三步、模型基本假定分析
模型基本假定分析时，通常包括‘断点适用性检验’和‘局部平滑性检验’。关于‘断点适用性检验’如果说cutoff值两侧附近的断点样本量基本均匀则说明断点选择适合，不受人为操纵。从第二步中得到的图形也可以看出，断点值附近两侧的点基本差不多，说明当前案例设置的断点值准确，并没有受到人为操纵干扰。
除此之外，还需要查看‘局部平滑性检验’，即分别将控制变量作为驱动变量X进行断点回归，通过图示法查看断点值是否在控制变量身上也起效果，即‘同样的断点值不应该在控制变量身上也起作用’，此检验通常并不完全需要。并且有时候控制变量并不能被当前断点cutoff值区分为两侧因而不能进行分析，本案例即是此类情况，本案例不进行‘局部平滑性检验’。
第四步、模型分析
第一步确认好模型为精确断点，并且第二步确定为二阶模型最优，而且满足基本模型假定。因而进行操作，准备得到最终结果。操作如下图所示：
分别设置结果变量y，驱动变量x，2个控制变量。以及设置好断点值为0（默认不设置即为0），设置为‘二阶’，选中‘稳健性检验’复选框。至于另外2个参数（带宽值、核函数）默认即可，并且将‘绘图’复选框取消（因为已经不再需要通过图示查看阶数）。此步骤为了得到最终结果。见‘SPSSAU输出结果’部分说明。
第五步、模型稳健性检验
在得到模型最终结果时，选中‘稳健性检验’复选框，系统默认提供不同带宽值（0.25倍、0.5倍、0.75倍、1倍、1.25倍、1.5倍、1.75倍和2倍共8个不同带宽值）时的结果，便于进行稳健性检验查看，实际研究中，可能并不需要8个不同带宽值情况下的结果对比，通常只需要1倍带宽值附近（比如0.75倍、1倍和1.25倍）共3项带宽值时结果对比，如果结论基本稳定即说明模型具有稳健性。
模型稳健性检验并没有固定的做法，只要可以证明模型具有稳健性（不同情况下模型结论基本一致则说明具有稳健性），具体稳健性方式上有很多种，一般使用1种或2种即可并没有固定标准。
至于其它的方式，比如‘更换核函数’法，‘更换断点’法，‘是否加入控制变量’法，‘更换阶数’法和‘样本选择法’。研究者可自行更换模型进行结果对比研究。尤其是‘更换核函数’法，‘是否加入控制变量’法和‘更换阶数’这3种方式，其操作简单方便，只需要在SPSSAU系统中下拉选择下参数更换即可进行，建议研究者尝试使用查看对比即可。比如‘更换核函数法’，操作截图如下所示：
SPSSAU默认是使用‘triangular三角核函数’，可选为‘Epanechnikov核函数’和‘Uniform核函数’，来回切换另外两个核函数，将结果进行汇总对比即可，如果结论基本一致则说明模型具有稳健性。
4 .SPSSAU输出结果
针对本案例结果，即‘精确断点’且‘二阶模型’时结果，SPSSAU共输出表格和图形，具体说明如下：
名称 | 说明 |
RDD基本情况 | 展示结果变量、处理变量、控制变量或fuzzy项对应的项名称 |
RDD参数情况 | 展示模型的参数值情况，包括断点值，带宽值，核函数，阶数和是否进行稳健性检验 |
RDD样本数据情况 | 展示RDD断点模型，断点值及断点两侧的样本量情况 |
RDD断点回归结果汇总 | 展示模型最终结果 |
稳健性检验-不同带宽时 | 展示不同带宽值时，模型汇总结果 |
稳健性检验（基于不同带宽） | 不同带宽值时模型回归系数95%置信区间coefplot图 |
5 文字分析
本案例得到最终结果，包括RDD基本情况、RDD参数情况、RDD样本数据情况、RDD断点回归结果汇总，稳健性检验结果和稳健性检验coefplot图，分别说明如下：
RDD基本情况 | |
项 | 名称 |
结果变量 | y |
处理变量或驱动变量 | x |
控制变量 | control1 |
control2 | |
Fuzzy模糊项 | - |
从上表格可以看到结果变量、处理变量、控制变量或fuzzy项对应的项名称，本案例中有两个控制变量，另本案例最终为精确断点，因而没有设置fuzzy模糊项。
RDD参数情况 | |
参数 | 名称 |
断点值 | 0 |
带宽值 | 0.096 |
核函数 | triangular |
阶数 | 2阶 |
稳健性检验 | 是 |
从上表可以看到，断点值为0，并且没有设置过带宽值，模型自动计算出‘最优带宽值’为0.096，并且默认使用三角triangular核函数，模型设定为2阶。以及选中‘稳健性检验’，SPSSAU自动会提供不同带宽值时模型汇总结果。
RDD样本数据情况 | |
项 | 数据 |
断点值 | 0 |
断点值左侧样本量 | 131 |
断点值右侧样本量 | 218 |
分析样本量 | 349 |
上表格可以看到，断点值为0，断点值左侧样本量为131个，右侧为218个，总共分析样本量为349个。
上表格可模型结果核心表格，从上表格可以看到，p值均大于0.05，但是小于0.1，也即意味意味着如果以0.1作为显著性水平，那么回归系数呈现出显著性，如果以0.05作为标准，则说明没有显著性。无论是Conventional法，也或者校正bias法（Bias-Corrected），也或者稳健法检验robust时。
提示：断点回归时对于回归系数的检验共提供3种方式，分别是Conventional法、Bias- Corrected和Robust法，三种方式并没有优劣之分。通常使用其中一种即可，比如Conventional法。下述基于不同带宽稳健性检验时默认汇总Conventional法。
由于本案例时选择‘稳健性’检验，因而SPSSAU提供上表格展示不同带宽值（0.25倍、0.5倍、0.75倍、1倍、1.25倍、1.5倍、1.75倍和2倍共8个不同带宽值）时回归系数显著性检验结果，默认汇总Conventional法的显著性检验结果。从上表格可以看到，8种情况下时，只有其中4种带宽下呈现出0.1水平显著性，另外远离1倍较远的带宽时并没有呈现出显著性。
整体上看，如果模型以0.1作为显著性水平，那么模型具有一定的稳健性（如果是0.05作为标准，则稳健性非常强，因为全部p值均大于0.05）。建议还可进一步通过其它方式，比如‘更换核函数’法，‘是否加入控制变量’法和‘更换阶数’等进一步查看。本案例中如果使用‘更换核函数’，‘是否加入控制变量’或者‘更换阶数’，也会有出现0.1水平显著的结论（但并不完全是）,但全部均会出现0.05水平不显著的结论，即意味着如果模型以0.1作为显著性标准，此时模型稳健性较弱，而模型以0.05水平作为标准，此时模型稳健性非常强，无论如何显著性值均大于0.05，最终模型以0.05作为显著性水平，即意味着模型并不显著，即‘民主党当选对于联邦支出并没有实际性影响’，并且此结论非常稳健。
上图为基于不同带宽时，模型回归系数95%置信区间进行展示的coefplot图，从图中可以看到，8种不同带宽时，95%置信区间均包括数字0，即明显的可以看到，模型在0.05水平上不显著，此结论稳定。
6 剖析
涉及以下几个关键点，分别如下：
ü RDD断点回归时，分析步骤较多，建议逐步进行，且在判断时尽量多的对比综合分析，而不能只查看某一个模型基础上得到最终结论；
（1）多数情况下使用精确断点，如果确实有非常强的理由证明应该使用模糊断点，也可使用模糊断点；
（2）模型阶数判断上结合图示进行，但带有一定的主观性，建议对比选择；
（3）模型假定分析时，图示直观上满足即可，不太可能模型完美的满足；
（4）模型分析时有3种显著性检验方式，选择其中一种即可；
（5） 模型稳健性检验有很多种方式，通常选中一个或者两个即可。
