倾向得分匹配法（PSM）举例及stata实现

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

      本例选自Cameron&Trivedi《微观计量经济学：方法与应用》（中译本，上海财经大学出版社，2010）pp794-800

      所有数据及程序均来自于本书的配套网站本文是末学向两位善知识Cameron和Trivedi的教材恭敬学习后整理的读书笔记，没有任何创意，只是简化整理一番而已，整理中难免错谬之处，敬请指正，不胜感激。末学联系方式houhh1125@163.com

倾向打分+选用得分进行匹配

     

      其中倾向得分匹配法的 stata 实现过程如下：

      总的讲来分为两大部分：

一是为倾向打分，

二是运用得分进行样本匹配并比较。

      由于 stata11 没有倾向得分匹配法相关的命令，所以要通过网络下载相关的 命令，相关的内容见附录。

一、倾向打分

---

第一步：设定宏变量 breps 表示重复抽样 200 次

1. . global breps 200

复制代码

第二步，设定宏变量 XDW02，

表示变量 AGE AGESQ EDUC EDUCSQ MARR NODEGREE BLACK HISP RE74 RE75 RE74SQ U74 U75 U74HISP

1. . global XDW02 AGE AGESQ EDUC EDUCSQ MARR NODEGREE BLACK HISP RE74 RE75 RE74SQ U74 U75 U74HISP

复制代码

第三步，通过 logit 模型进行倾向打分

1. . pscore TREAT $XDW02, pscore(myscore) comsup blockid(myblock)
   
2. numblo(5) level(0.005) logit

复制代码

二、运用得分进行样本匹配并比较

---

方法一： Nearest neighbor matching (random version) 最邻近方法

首先，产生随机数种子

1. . set seed 10101

复制代码

然后

1. . attnd RE78 TREAT $XDW02, comsup boot reps($breps) dots logit

复制代码

【说明】注意变量排列次序， RE78 正是我们所关心的变量，放在最前面； 其次是处理指示变量，最后是其他变量。关于选项内容，使用 help attnd 命令可 以详细学习。

方法二： Radius matching（半径匹配法）

1. . set seed 10101 . attr RE78 TREAT $XDW02, comsup boot reps($breps) dots logit radius(0.001)

复制代码

方法三： Stratification Matching（分成匹配法）

1. . set seed 10101 . atts RE78 TREAT, pscore(myscore) blockid(myblock) comsup boot reps($breps) dots

复制代码

方法四： Kernel Matching（核匹配方法）

1. . set seed 10101 . attk RE78 TREAT $XDW02, comsup boot reps($breps) dots logit

复制代码

三、附件

---

倾向得分匹配法（PSM）举例及stata实现

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏66 回帖

关键词：倾向得分匹配法 倾向得分匹配 Stata 倾向得分 tata 倾向得分 PSM stata 倾向得分匹配法

感谢分享

好贴，感谢分享。

现在不是有专门的命令了吗啊

谢谢分享

玄一无相 发表于 2015-8-4 14:35
现在不是有专门的命令了吗啊

请问“专门的命令”是什么？？

谢谢楼主

多难兴邦02 发表于 2016-3-25 16:53
请问“专门的命令”是什么？？

可以用psmatch2

谢谢楼主

merlinny96 发表于 2016-4-17 18:41
可以用psmatch2

请问哪个版本的stata直接有这个命令 不需要package么