PSM

set seed 10000 ///定义种子

gen tmp = runiform() ///生成随机数

sort tmp ///将数据库随机整理

gen treated = (idcode > 2000) &! missing(idcode) ///设置实验组和控制组

*用二值选择模型 logit 回归估计倾向值，并基于近邻匹配（默认 K=1）实现一对一匹配；

psmatch2 treated $xlist,out(y) logit ate neighbor(1) common caliper(.05) ties 

pstest $xlist, both graph

gen common = _support ///检验协变量在处理组与控制组之间是否平衡

drop if _weight==. ///去掉不满足共同区域假定的观测值

用原回归模型回归

©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容