专栏名称: 连享会
连玉君老师团队分享，主页：lianxh.cn。白话计量，代码实操；学术路上，与君同行。
目录
相关文章推荐
笔吧评测室  ·  小米首款 AI PC ... ·  昨天  
笔吧评测室  ·  戴尔全新命名 14/16 Plus ... ·  2 天前  
笔吧评测室  ·  新一代 ROG XG 显卡扩展坞发布：可选 ... ·  3 天前  
笔吧评测室  ·  惠普将推 RTX 5060/5050 ... ·  3 天前  
笔吧评测室  ·  聊一个“知名老品牌”推出的新型笔记本 ·  3 天前  
51好读  ›  专栏  ›  连享会

Stata: psestimate-倾向得分匹配(PSM)中匹配变量的筛选

连享会  · 公众号  ·  · 2024-12-18 22:00

正文

👇 连享会 · 推文导航 | www.lianxh.cn

🍓 课程推荐： 连享会：2025 寒假班
嘉宾：连玉君（初级|高级）；杨海生（前沿）
时间：2025 年 1 月 13-24 日
咨询：王老师 18903405450（微信）

温馨提示: 文中链接在微信中无法生效。请点击底部 「阅读原文」 。或直接长按/扫描如下二维码，直达原文：



倾向得分匹配分析 (PSM) 已经在诸多领域得到了应用。虽然 PSM 不能完全解决内生性问题，但却能在很大程度上缓解自我选择问题导致的偏差。在前期文献中，Becker & Ichino (2002, Stata Journal, 2(4):358-377) 对 PSM 的分析过程进行了详细的介绍，Stata 中也有多个命令可以执行 PSM 分析，如 pscore , psmatch2 , treatrew (Stata Journal, 14(3): 541-561), gpscore (SJ 8(3):354--373), kmatch 

  1.  net describe st0328, from(http://www.stata-journal.com/software/sj14-1)

平衡性假设

在 PSM 匹配时，用treat变量对控制变量进行Logit回归，得到倾向得分值。倾向得分值最接近的控制组个体即为实验组的配对样本，通过这种方法可以最大程度减少实验组与控制组个体存在的系统性差异，从而减少估计偏误。在进行PSM匹配后的其他估计前比如PSM-DID 估计前，还需进行协变量的平衡性假设检验，即匹配后各变量在实验组和控制组之间是否变得平衡，也就是说实验组和控制组协变量的均值在匹配后是否具有显著差异。如果不存在显著差异，则支持进一步的模型估计。

在平衡性检验之前，我们先使用 psmatch2 命令进行PSM匹配，处理变量为 train ，协变量为 age 、 educ 、 black ，结果变量为 re78 ，采用一对一近邻匹配，具体操作如下： 

  1.  use ldw_exper.dta,clear
  2.  psmatch2 train age educ black, out(re78) logit ate neighbor(1) common caliper(.05) ties

PSM 匹配完成之后，我们需要检验匹配后的样本是否满足平衡性假设，即实验组与控制组的匹配协变量是否没有显著性差异，在这里可以使用 pstest 命令进行检验，具体如下： 

  1.  pstest age educ black hisp married , t(train)

平衡性假设检验结果如下： 

  1. ------------------------------------------------------------------------------
  2.                         |       Mean               |     t-test    |  V(T)/
  3. Variable                | Treated Control    %bias |    t    p>|t| |  V(C)
  4. ------------------------+--------------------------+---------------+----------
  5. age                     | 25.527   24.714     11.4 |   1.19  0.234 |  1.24
  6. educ                    | 10.291   10.401     -6.0 |  -0.59  0.557 |  1.60*
  7. black                   | .84066   .87363     -8.9 |  -0.90  0.370 |     




    
.
  8. hisp                    | .06044   .09066    -10.9 |  -1.09  0.277 |     .
  9. married                 | .18681    .1522      9.2 |   0.88  0.380 |     .
  10. ------------------------------------------------------------------------------

根据t检验结果发现，以上5个协变量在实验组与控制组之间不存在显著性差异。

那么，在进行 PSM 分析之前，应当如何选择匹配协变量，使模型实现最佳的拟合效果呢？今天介绍的 psestimate 命令可以通过比较不同模型的极大似然值，帮助我们选择能实现最佳拟合效果的协变量的一阶和二阶形式。

The psestimate command estimates the propensity score proposed by Imbens and Rubin (2015).  The main purpose of the program is to select a linear or quadratic function of covariates to include in the estimation function of the propensity score.

1. 命令的安装与示例数据导入

在Stata命令窗口执行第一行代码即可完成对 psestimate 命令的下载，然后输入第二行命令下载网上示例数据 nswre74.dta（LaLonde, 1986），并执行第三行命令导入数据。 

  1. ssc install psestimate, replace //安装命令
  2. net get psestimate  //下载命令附带的数据到当前工作路径下
  3. use "psestimate.dta", replace

2. 命令的语法

该命令的语法如下： 

  1.  psestimate depvar [indepvars] [if] [in] [, options]
  2.  options：             




    

  3.       totry(indepvars)     
  4.       notry(varlist)       
  5.       nolin
  6.       noquad
  7.       clinear(real)       
  8.       cquadratic(real)    
  9.       iterate(#)          
  10.       genpscore(newvar)   
  11.       genlor(newvar)

各个主要选项的含义如下：

  • depvar ，必选项，填入处理变量（如 treat），即标记是否参与实验的虚拟变量

  • indepvars ，可选项，指定基准模型中的协变量

  • totry(indepvars) ，可选项，放入供选择的协变量列表，默认为全部

  • notry(varlist) ，可选项，指定不包括的协变量列表，默认为没有

  • nolin ，可选项，指定不进行一阶多项式的选择

  • noquad ，可选项，指定不进行二阶多项式的选择

  • clinear(real) ，可选项，指定一阶协变量似然比检验的门槛值，默认值为 1

  • cquadratic(real) ，可选项，指定二阶协变量似然比检验的门槛值，默认值是 2.71

  • iterate(#) ，可选项，指定在每个 logit 中执行循环的最大值，默认值是 16000

  • genpscore(newvar) ，可选项，由于指定程序自动生成的用于记录倾向得分值的新变量的名称

  • genlor(newvar) ，可选项，生成对数似然比的新变量的名称

3. 命令操作

3.1 命令基本操作

下面本文将基于 psestimate 命令的作者提供的数据集 nswre74.dta 来简要说明如何使用 psestimate 这一命令来选择能最好拟合处理变量 (treat) 的协变量的一阶及二阶形式。

在这里，我们事先选定教育变量 ed 作为基准模型中的一个协变量，意味着 Stata 自动将 ed 放入基准模型中。另外，我们将指定 age 、 black 、 hisp 、 nodeg 四个变量作为待选协变量。代码如下： 

  1. use "nswre74.dta", clear
  2. psestimate treat ed, totry(age black hisp nodeg)

运行结果如下： 

  1. Selecting first order covariates... (10)
  2. ----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
  3. ...s..s..
  4. Selected first order covariates are: nodeg hisp
  5. Selecting second order covariates... (21)
  6. ----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
  7. .....s.....
  8. Selected second order covariates are: c.nodeg#c.ed
  9. Final model is: ed nodeg hisp c.nodeg#c.ed

根据以上结果，可以确定在倾向得分匹配中，我们应该选取的一阶协变量为 nodeg 、 hisp ，二阶协变量为 c.nodeg#c.ed 。综上，根据 psestimate 命令的运算结果，我们应该选取 ed 、 nodeg 、 hisp 、 c.nodeg#c.ed 等四个变量作为倾向得分匹配的协变量。

3.2 提升运算速度

psestimate 命令在运算中会耗费较长时间，而通常来说，该命令在选择协变量的一阶形式时要快于二阶形式的选择，因此，为了加快运算速度，我们可以首先通过加入 noquad 选项，只对协变量的一阶形式进行筛选，当一阶形式选定后，将其作为解释变量放入基准模型中，然后加入 nolin 选项跳过一阶形式筛选步骤，只对协变量的二阶形式进行筛选。具体操作如下。

首先，加入入 noquad 选项，只筛选协变量的一阶形式，如下： 

  1. use "nswre74.dta", clear
  2. psestimate treat ed, totry(age black hisp nodeg) noquad

一阶协变量的筛选结果如下： 

  1. Selecting first order covariates... (10)
  2. ----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
  3. ...s..s..
  4. Selected first order covariates are: nodeg hisp
  5. Final model is: ed nodeg hisp

然后，将选定的 ed 、 nodeg 、 hisp 作为解释变量放入基准模型中，加入 nolin 选项值进行二阶形式的筛选，操作如下：。 

  1. psestimate treat ed nodeg hisp , totry(age black hisp nodeg) nolin

二阶协变量的筛选结果如下： 

  1. Selecting second order covariates... (21)
  2. ----+---




    
 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
  3. .....s.....
  4. Selected second order covariates are: c.nodeg#c.ed
  5. Final model is: ed nodeg hisp  c.nodeg#c.ed

4. psestimate 的核心思想

4.1 协变量一阶形式的选择

第一步，该程序首先在基准模型（ logit treat ed ）基础上通过循环分别加入 totry() 中指定的四个变量 age 、 black 、 hisp 、 nodeg ，进行四次模型估计，如下所示： 

  1. logit treat ed age
  2. logit treat ed black
  3. logit treat ed hisp
  4. logit treat ed nodeg

每次估计完成后，它将得到的新的极大似然值与基准模型比较，选择上述四个模型中对数极大似然值 (Log-Likelihood, 简称 LL 值) 最大的模型中的协变量放入基准模型中，除非上述四个模型的极大似然值都低于 clinear(real) 中指定的门槛值。若此处假设为 nodeg ，则基准模型扩展为 logit treat ed nodeg , 然后第二步，它将估计如下模型： 

  1. logit treat ed nodeg age
  2. logit treat ed nodeg black
  3. logit treat ed nodeg hisp

这一步的协变量筛选原则与第一步相同。可以看出，当供选择的协变量数量为 时，在确定协变量的一阶形式时，该程序理论上会估计 个 Logit 模型。本例中有 4 个供选择的协变量，则需要估计 10 次（如下括号中所示），该命令选择的协变量一阶形式结果如下： 

  1. Selecting first order covariates... (10)
  2. ----+---




    
 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
  3. ...s..s..
  4. Selected first order covariates are: nodeg hisp

4.2 协变量二阶形式的选择

在协变量二阶形式的选择上，主要分为协变量平方项以及协变量间的交乘项。
如果在一阶形式中只选择了 a 这一个协变量，则二阶形式的选择只需要检验 a^2 这一变量。但是如果有 a 、 b 两个一阶协变量被选择，则二阶形式的选择需要检验 a^2 、 b^2 、 ab 三个二阶协变量形式。具体到本例，确定的一阶协变量有 ed 、 nodeg 、 hisp 三个，则需要检验的二阶协变量有六个，即 ed^2 、 nodeg^2 、 hisp^2 、 c.ed#c.nodeg 、 c.ed#c.hisp 、 c.nodeg#c.hisp ，筛选过程与选择协变量一阶形式的方法一致。因此本例中共需估计 即 21 次（如下括号中所示），结果如下所示： 

  1. Selecting second order covariates... (21)
  2. ----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
  3. .....s.....
  4. Selected second order covariates are: c.nodeg#c.ed
  5. Final model is: ed nodeg hisp c.nodeg#c.ed

4.3 流程图展示

如下流程图可以更加直观地展现 psestimate 筛选协变量一阶及二阶形式的过程，为简化分析，我们可供选择的协变量为 a 、 b 两个变量，假设各模型的对数极大似然值存在如下大小关系，LL1>LL2> clinear() >LL3,LL4>LL5>LL6> cquadratic() >LL7>LL8。

图1.流程图

5. PSM估计的完整流程示例

5.1 psestimate 筛选匹配变量的一阶、二阶形式

第一步，使用 psestimate 筛选匹配变量 

  1. use "nswre74.dta", clear
  2. psestimate treat ed, totry(age black hisp nodeg)

匹配变量选择如下： 

  1. Selecting first order covariates... (10)
  2. ----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
  3. ...s..s..
  4. Selected first order covariates are: nodeg hisp
  5. Selecting second order covariates... (21)
  6. ----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5 
  7. .....s.....
  8. Selected second order covariates are: c.nodeg#c.ed
  9. Final model is: ed nodeg hisp c.nodeg#c.ed

最终选择的匹配变量为 ed 、 nodeg 、 hisp 、 c.nodeg#c.ed

5.2 psmatch2 基于筛选出的匹配变量进行PSM匹配

基于上述匹配变量进行PSM匹配： 

  1.  psmatch2 treat ed nodeg hisp c.nodeg#c.ed, logit ate neighbor(1) common caliper(.05) ties

结果如下： 

  1. Logistic regression                             Number of obs     =        445
  2.                                                 LR chi2(4)        =      17.03
  3.                                                 Prob > chi2       =     0.0019
  4. Log






请到「今天看啥」查看全文


推荐文章
笔吧评测室  ·  小米首款 AI PC 官宣，卢伟冰透露新笔记本搭载 99Wh 超大电池、刚好能上飞机
昨天
笔吧评测室  ·  戴尔全新命名 14/16 Plus 轻薄本发布：酷睿 Ultra 5 226V 型号 6999/7499 元起
2 天前
笔吧评测室  ·  新一代 ROG XG 显卡扩展坞发布：可选 RTX 5070Ti / 5090，雷电 5 接口，1999 美元起
3 天前
笔吧评测室  ·  惠普将推 RTX 5060/5050 显卡光影精灵游戏本，基于 Raptor Lake-H 处理器
3 天前
笔吧评测室  ·  聊一个“知名老品牌”推出的新型笔记本
3 天前
创意铺子  ·  这样的神兽你见过？
8 年前
吴晓波频道  ·  月薪3万的白领与扫地阿姨，谁是真正的中产？｜荐读
7 年前
上海陆家嘴并购联盟  ·  周小川陆家嘴论坛主旨演讲全文：保护导致懒惰 中国金融服务业要进一步扩大开放
7 年前
哈佛商业评论  ·  推广 ｜tutorabc CTO汤峥嵘：人工智能时代，技术能够真正改变教育
7 年前
南宁财神哥  ·  自动赚钱机器，零成本创业项目，已经有人实现屌丝逆袭了，你就看看吗？
7 年前
Sov5搜索   ·   小百科   ·   今天看啥   ·   移动版
51好读 - 好文章就要读起来!