模型系列-PSM原理介绍

2019-04-12 本文已影响106人凡有言说

这篇文章是和小伙伴LF协作完成的

倾向值匹配(propensity score matching,PSM)

一、为什么需要匹配？

随机实验很好，但成本高，做起来有时无法落地，不切实际。比如为了研究大学教育对收入的影响，不可能找一批已经考上大学的学生，随机分为两组，一组是上大学，一组不去上大学，然后等毕业再来看两组的收入差异吧。

在多数情况下，研究者只有观测数据，没有实验数据，处理组和控制组的成员不是由随机分配产生，而是可以自由选择参与其中。正是由于选择效应的存在，并非随机分组，倾向值匹配(propensity score matching,PSM)就有用武之地啦。

针对非随机分组的个体，匹配估计量的基本思路是：找到属于控制组的某个个体j，使其与属于处理组的个体i在可测变量的取值上尽可能的相似（“匹配”），即X_i = X_j

又由于给定X_i，Y_0i，Y_1i对立于D_i ,因此个体i与j进入处理组的概率详尽，具有可比性。

因此可以将Y_j作为Y_0i的估计量；进一步，可以将Y_1i-Y_j作为对个体i的处理效应。

什么是匹配？

匹配是一种非实验方法，是对于一些没有采用或不方便采用实验方法区分处理组和控制组的数据采用的一种近似实验的方法。

匹配方法假定，控制协变量之后，具有相同特征的个体对政策具有相同的反应。换句话说，不可观测因素不影响个体是否接受政策干预的决策，选择仅仅发生在可观测变量上。

因此，对每一个实验组个体而言，可以根据可观测特征为其选择一个控制组个体构成反事实。

匹配的目的在于确保干预效应估计，是建立在可比个体之间的不同结果基础上。最简单的匹配方式是将处理组和控制组中协变量值相同的两个个体进行配对分析。

举个栗子

例子

对于控制组的个体1，由于X₁ = X₅ =2 ，因此个体1与处理组的个体5匹配。Y₀₁的估计量为7，Y₁₁的估计量为8；

对于控制组的个体2，没有相同的匹配，只有近似的匹配X₄ = X₆ =3。Y₀₂的估计量为8，而Y₁₂的估计量为(Y₄ + Y₆)/2 = 7.5

对于整个样本的匹配结果，此时需要计算平均处理效应（又称“平均因果效应ACE”）：

ATE = E(Y_1i-Y_0i) =((8-7)+(7.5-8)+(7.5-6)+(9-7.5)+(8-7)+(6-7.5)+(7-5))/7 = 0.143

但我们更关心参加培训的匹配结果，即只需要计算参与者平均处理效应（又叫“参与者处理效应TOT”）
ATT = E(Y_1i-Y_0i|D_i=1)=((9-7.5)+(8-7)+(6-7.5)+(7-5))/4=0.25

但如果协变量不是某一个变量，而是一组变量时，这种简单的匹配方式就不再适用。

在现实生活中，数据偏差和混杂变量比较多。由于存在很多其他变量混淆自变量和因变量之间的关系，研究者很难直接探索二者间的净效果。这些混淆变量的影响通常被称为选择性误差，需要通过倾向值匹配的方法来控制和消除。

倾向值指的是，被研究个体在控制可观察到的混杂变量的情况下，接受某种干预的条件概率。

给定X_i，个体i的倾向值（倾向得分）为个体i进入处理组的条件概率，即P(X_i) = P(Di=1|X=X_i)

以D_i={0,1}表示个体i是否参加培训

PSM的实质就是把许多可观察到的混杂变量整合成一个变量：倾向值。

由于具有相同或相近倾向值的个体的其他变量在分布上具有相同的特征，故将处理组和控制组的个体根据倾向值进行匹配，从而平衡两组样本的基线数据，达到类似随机分组的效果。

由于混杂变量在倾向值匹配的过程中被控制起来，两组个体在结果上的差异就只能归因于有无干预措施。

回到大学教育对收入影响的例子，PSM想要解决的问题就是，由于A同学已经读了大学，如何估计出A要是不读大学，A的收入会是多少？**

倾向值匹配能从样本中，对每个人读大学的概率进行估计。然后选出与 A 同学有相似念大学的概率，却没有去读的同学 B 作为 A 同学的对照，然后再来看他们的区别。当样本中的每位“大学生 A” 都找到了匹配的“非大学生 B” ，就能对这两组样本进行比较研究了。

关于stata操作放在下一篇推文

参考资料：

公众号推荐

阿虎定量笔记