PSM法匹配后如何操作

27683

收藏 2015-12-31

悬赏 100 个论坛币已解决

情况是这样的：目的是要分析处理组和对照组患病率的比较，现已通过PSM法完成了匹配。想问：下一步该怎么做才能回到目的分析患病率的变化？难道要把处理组和与之匹配的对照组抽取出来成为一个新的样本吗？如果是，该怎么做？悬赏求助，价格虽不高，但已是我的极限，拜托各路大神非大神都进来看看、讨论讨论，跪谢！

最佳答案

gingerlinger 查看完整内容

%bias在10%以内为佳，多数match后的bias小于unmatch的

扫码加我拉你入群

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

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

全部回复

gingerlinger

2015-12-31 16:07:11

%bias在10%以内为佳，多数match后的bias小于unmatch的

扫码加我拉你入群

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

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

hzl619

2015-12-31 16:21:58

不知可否让你满意，https://bbs.pinggu.org/thread-1490982-1-1.html

扫码加我拉你入群

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

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

wenwubei659

2016-1-1 10:24:44

hzl619 发表于 2015-12-31 16:21
不知可否让你满意，https://bbs.pinggu.org/thread-1490982-1-1.html

这个我看了，但是还是不明白，谢谢

扫码加我拉你入群

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

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

gingerlinger

2016-1-1 21:16:03

直接用psmatch2

. use 你的样本链接

. set seed 10101
. gen ranorder = runiform()
. sort ranorder
//这是把顺序打乱，便于最邻近随机匹配
//接下来匹配
. psmatch2 （是否处理变量）（所有协变量）, out ( 患病率变量) n(1) ties common
//n(1)是一对一匹配，ties 是如果得分相同，随机选择一个匹配，common是只对common support的部分进行匹配，psm的匹配方法很多，不知道你是用一对一，还是近邻匹配还是kernel匹配还是局部匹配还是卡尺匹配，根据我经验kernel和卡尺效果比较好，我给你的是一对一匹配
接下来结果里面你就找ATT对应的difference就是你要的处理效应，（处理会增加患病率或者减少患病率，以及增减的多少）
然后还可以用-pstest-测试匹配效果如何
psgragh把图画出来
看你具体需要

扫码加我拉你入群

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

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

wenwubei659

2016-1-3 16:18:57

gingerlinger 发表于 2016-1-1 21:16
直接用psmatch2

. use 你的样本链接

你好，以下是我运行的命令，其中V8表示是否处理变量（取值为0，1）；V32 V33 V34 V35 V36 V747 d1为协变量；V110为就诊变量（取值为1，2。1代表就诊，2代表未就诊）。这样的话我是不是就可以得到就诊率的比较？
set seed 10101
gen ranorder = runiform()
sort ranorder
logistic V8 V32 V33 V34 V35 V36 V747 d1
predict pscore, p
psmatch2 V8,out(V32 V33 V34 V35 V36 V747 d1) pscore(pscore) neighbor(1) caliper(0.05) common odds index logit ties nowarnings quietly ate
pstest V32 V33 V34 V35 V36 V747 d1,both
attnd V110 V8 V32 V33 V34 V35 V36 V747 d1,comsup boot reps(2) dots logit

最后一个命令得出的结果中，ATT值是否可理解为 “处理组就诊率 — 对照组就诊率”？

扫码加我拉你入群

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

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

点击查看更多内容…

wenwubei659

2016-1-4 20:05:46

gingerlinger 发表于 2016-1-1 21:16
直接用psmatch2

. use 你的样本链接

你好，我想问一下，pstest检验的结果怎么解释呢？谢谢，这个回答完后我就设为最佳答案

扫码加我拉你入群

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

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

wenwubei659

2016-1-4 20:05:49

gingerlinger 发表于 2016-1-1 21:16
直接用psmatch2

. use 你的样本链接

你好，我想问一下，pstest检验的结果怎么解释呢？谢谢，这个回答完后我就设为最佳答案

扫码加我拉你入群

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

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

gingerlinger

2016-2-15 20:44:20

wenwubei659 发表于 2016-1-4 20:05
你好，我想问一下，pstest检验的结果怎么解释呢？谢谢，这个回答完后我就设为最佳答案

%bias在10%以内为佳，多数match后的bias小于unmatch的

扫码加我拉你入群

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

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

wenwubei659

2016-3-6 18:24:20

gingerlinger 发表于 2016-1-1 21:16
直接用psmatch2

. use 你的样本链接

请问画密度图的命令是什么？

扫码加我拉你入群

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

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

wenwubei659

2016-3-15 19:32:35

gingerlinger 发表于 2016-2-15 20:44
%bias在10%以内为佳，多数match后的bias小于unmatch的

你好！请教一下，为什么运行卡尺匹配（0.05）等了好半天都不出结果？样本量在2～3万。
希望回答一下，毕竟悬赏都给你了，谢谢～！

扫码加我拉你入群

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

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

gingerlinger

2016-4-13 01:25:59

wenwubei659 发表于 2016-3-15 19:32
你好！请教一下，为什么运行卡尺匹配（0.05）等了好半天都不出结果？样本量在2～3万。
希望回答一下， ...

天哪我才看到额，没做过这么大的样本量，不过我有同学要等一晚上的也有

扫码加我拉你入群

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

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

gingerlinger

2016-4-13 01:26:49

wenwubei659 发表于 2016-3-6 18:24
请问画密度图的命令是什么？

psgraph?这么久了都快忘光了。。。。你试试这个

扫码加我拉你入群

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

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

gingerlinger

2016-4-13 01:37:03

wenwubei659 发表于 2016-1-3 16:18
你好，以下是我运行的命令，其中V8表示是否处理变量（取值为0，1）；V32 V33 V34 V35 V36 V747 d1为协变量 ...

ATT是average effect of treatment on the treated
如果是随机试验的话，有ATE=ATT
记不大清了。。。要看你处理和对照是不是随机抽取的。。。。

扫码加我拉你入群

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

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

恋上香草味

2016-8-31 15:47:13

gingerlinger 发表于 2016-1-1 21:16
直接用psmatch2

. use 你的样本链接

你好，你的答案让我们受益匪浅，想再咨询一下，如果做完psmatch2，但是本文想把样本分为三个：总样本，A样本和B样本，分别看它的ATT差异。应该如何加入条件变量？非常感谢

扫码加我拉你入群

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

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

TeresaTLH

2016-11-29 08:40:07

恋上香草味发表于 2016-8-31 15:47
你好，你的答案让我们受益匪浅，想再咨询一下，如果做完psmatch2，但是本文想把样本分为三个：总样本，A样 ...

同问，请问应该怎么做呀？

扫码加我拉你入群

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

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

TeresaTLH

2016-11-29 08:41:25

gingerlinger 发表于 2016-4-13 01:25
天哪我才看到额，没做过这么大的样本量，不过我有同学要等一晚上的也有

请问如果处理组样本量是2000，控制组是600，这样的样本可以做吗？

扫码加我拉你入群

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

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

gingerlinger

2017-1-16 03:44:03

恋上香草味发表于 2016-8-31 15:47
你好，你的答案让我们受益匪浅，想再咨询一下，如果做完psmatch2，但是本文想把样本分为三个：总样本，A样 ...

太久远了。。。学完就忘。。。。其实没大看懂你的问题意思。。。

扫码加我拉你入群

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

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

gingerlinger

2017-1-16 03:44:40

TeresaTLH 发表于 2016-11-29 08:41
请问如果处理组样本量是2000，控制组是600，这样的样本可以做吗？

为啥不可以

扫码加我拉你入群

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

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

xiefengqi

2017-2-11 16:25:25

gingerlinger 发表于 2017-1-16 03:44
太久远了。。。学完就忘。。。。其实没大看懂你的问题意思。。。

你之前提问的匹配后的核密度指令解决了吗。连老师给的那个指令没有体现权重的思想啊，如果是用1对3匹配，或者核匹配，匹配后画出来的图基本上一样

扫码加我拉你入群

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

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

兰郡月光

2017-2-25 20:19:11

gingerlinger 发表于 2016-1-1 21:16
直接用psmatch2

. use 你的样本链接

请问 pstest出来的表格怎么看平稳性呢？是看p值由小变大就好吗？

扫码加我拉你入群

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

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

最佳答案

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群