Stata用diff命令做基于kernel psm的did，每次做出来的结果都不一样

llishaa 发表于 2018-12-18 17:06
黄老师，请教一下，我用PSM做匹配，处理组有70个样本，而待配对的控制组有2300多个，我如果使用1:1不放回 ...

应该是你做错了！

黃河泉 发表于 2018-12-18 17:16
应该是你做错了！

黄老师，我的命令没有问题啊，而且匹配的平衡性检验结果也能出来，那是什么地方错了呢，可否把数据发给您，帮忙给看看呢

llishaa 发表于 2018-12-18 17:36
黄老师，我的命令没有问题啊，而且匹配的平衡性检验结果也能出来，那是什么地方错了呢，可否把数据发给您 ...

你怎确认你的命令没错？请发出来！

黃河泉 发表于 2018-12-18 18:16
你怎确认你的命令没错？请发出来！

*\PSM相关的命令操作\*
*第一步，首先生成一个随机数种子
set seed 10101
gen ranorder=runiform()
sort ranorder

*\第二步，进行核匹配，默认Probit回归,一般使用logit回归，Rosenbaum and Rubin（1985）建议使用形式灵活的logit模型
psmatch2 treat 匹配变量1 匹配变量2 匹配变量3 匹配变量4 匹配变量5 匹配变量6 匹配变量7,outcome(y) kernel  ate ties common
    1:1不放回
psmatch2 treat 匹配变量1 匹配变量2 匹配变量3 匹配变量4 匹配变量5 匹配变量6 匹配变量7,outcome(y) n (1) norepl  ate ties common        1:1可放回
psmatch2 treat 匹配变量1 匹配变量2 匹配变量3 匹配变量4 匹配变量5 匹配变量6 匹配变量7,outcome(y) n (1)  ate ties common
    1:2
psmatch2 treat 匹配变量1 匹配变量2 匹配变量3 匹配变量4 匹配变量5 匹配变量6 匹配变量7,outcome(y) n (2)  ate ties common
    1:3
psmatch2 treat 匹配变量1 匹配变量2 匹配变量3 匹配变量4 匹配变量5 匹配变量6 匹配变量7,outcome(y) n (3)  ate ties common
    我分别用了这几种方法，匹配出来的控制组都是1000多，匹配偏差都在20%以下
*自助法得到标准误，即ATU和ATE
bootstrap r(att) r(atu) r(ate),reps(500):psmatch2 treat 匹配变量1 匹配变量2 匹配变量3 匹配变量4 匹配变量5 匹配变量6 匹配变量7,outcome(y) kernel  ate ties common

*检验匹配后数据是否平衡以及画图得出倾向得分的共同取值范围
pstest  匹配变量1 匹配变量2 匹配变量3 匹配变量4 匹配变量5 匹配变量6 匹配变量7,both graph  \\ 画出各变量的标准化偏差图示
psgraph  \\画图得出倾向得分的共同取值范围的直方图

llishaa 发表于 2018-12-18 19:18
*\PSM相关的命令操作\*
*第一步，首先生成一个随机数种子
set seed 10101

1. 1:1不放回
   
2. psmatch2 treat 匹配变量1 匹配变量2 匹配变量3 匹配变量4 匹配变量5 匹配变量6 匹配变量7,outcome(y) n (1) norepl  ate ties common        
   
3. 1:1可放回
   
4. psmatch2 treat 匹配变量1 匹配变量2 匹配变量3 匹配变量4 匹配变量5 匹配变量6 匹配变量7,outcome(y) n (1)  ate ties common

复制代码

不可能有 1000多 控制组，请确定。

黃河泉 发表于 2018-12-19 06:33
不可能有 1000多 控制组，请确定。

黃河泉 发表于 2018-12-19 06:33
不可能有 1000多 控制组，请确定。

黄老师，我把运行的结果发给您了，您帮我看看

llishaa 发表于 2018-12-19 10:23
黄老师，我把运行的结果发给您了，您帮我看看

哪一个？ 不就是 70 个吗？

黃河泉 发表于 2018-12-19 10:32
哪一个？ 不就是 70 个吗？

llishaa 发表于 2018-12-19 10:36

平衡检验结果出不来，r（2000）