倾向得分匹配的stata命令_[论文复制5]倾向得分匹配+双重差分（PSMDID）

最新推荐文章于 2024-08-22 10:15:15 发布

weixin_39885412

最新推荐文章于 2024-08-22 10:15:15 发布

阅读量4.4k

点赞数 4

文章标签：倾向得分匹配的stata命令

本文链接：https://blog.csdn.net/weixin_39885412/article/details/111374888

版权

[论文复制5]自然资源资产离任审计与空气污染防治：“和谐锦标赛”还是“环保资格赛”

PSM部分

上一期我们已经讲完了安慰剂检验，作者为展示其研究结论的稳健性，还使用了PSMDID这一方法。

PSM的思想是：假设个体

属于处理组，匹配估计量的基本思路是，找到数据控制组的某个体

，使得个体

与个体

的可测变量取值尽可能相似（匹配），即

。

关于PSM的详细介绍，可以阅读陈强教授编写的《高级计量经济学与Stata应用》。

使用PSM之前，需要对样本进行随机排序。此时，我们需要使用runiform() 函数生成一列均匀分布的随机数ranorder，再使用sort命令使其从小到大依次排列，此时，各样本即实现了随机排序。

细心的读者可能发现了，我们再生成随机数的时候，前面加了set seed 0001这行命令，这是为了设定一个抽样种子，使每次随机抽样的结果保持一致。因为不设定随机种子，生成的随机是将无法重现。当然，你也可以更改0001这一数字，使用其他数字。

cd D:stata15personal论文复制自然资源离任审计
use datafordid.dta, clear
xtset serial year

global xlist2 "lnpgdp lnpgdp2 popdst age edu gender tenure1 tenure5"
global xlist3 "coast mount temp water lnpgdp"

set seed 0001
gen ranorder = runiform()
sort ranorder

随后，即可开始对样本进行配对了。我们使用psmatch2进行配对（该命令为外部命令，需要使用ssc install psmatch2, replace 进行安装）。下面展示具体的运行命令。treat为处理变量，$xlist3为协变量，包括城市的自然地理特征和经济发展水平。out(aaqi) 用来指定aaqi为结果变量，logit是指用logit模型计算倾向得分。ate表示同时汇报ATE、ATU和ATT，默认仅汇报ATT。neighbor(1)