分享

SPSS实现1:M倾向性匹配

 松哥精鼎统计 2020-10-23
缘起
      倾向性评分火的不能再火,学习难度最低,使用人数最多的统计软件SPSS在22.0以上的版本,通过加载Python插件的方式,引入了1:1的匹配,虽然很好了,但是1:M始终是个缺憾。其实SPSS专门有个R插件,可以进行1:1和1:M,效果非常好,今天就给大家做期教程.
    松哥统计制作了35期倾向性评分视频课程,感兴趣可以订阅,新课早鸟,优惠进行中!

1.案例数据






    这是倾向性评分的经典哪里,stata和R软件都采用这个案例进行演示。不过多解释,只要知道红色框是本次干预措施,绿色为协变量,黄色为结局变量即可。

图1案例数据

2.SPSS-PSMatching






2.1打开PSMatching:

图2 SPSS-PSMatching

2.2PSMatching界面操作:

      图3中A位倾向性评分匹配算法,B为变量放入框。C为卡钳值设定框,默认为0.2。

图3 PSMatching主界面

2.3各种平衡图设置框

    图4中,可以选择PS直方图,散点图,标准化差异直方图,标准化均数差异点图,个体差异线图。

图4 各种平衡图复选框

2.4输出与检验设置

    图5A为匹配结果输出设置,可以成组式输出,也可以直接输出配对(配比)形式数据集,便于后续任性的分析。
    图5B为均衡性检验统计表,有基础检验和详细检验2方式,不用说,详细更加嗨哦!

图5输出数据与检验设置

2.5 1:M设置

      眼睛瞪大!眼睛瞪大!图6A处,就是设置1:M的哦!图6B处是选择从大到小,还是从小到大进行匹配的顺序的。

图6 1:M设置处

3.结果





3.1如何分析真实世界研究数据

    图7A告诉我们到底匹配了多少对。我们发现221-113,并不是2:1呀,是不是软件不靠谱呀!其实不是,答案是在给定的条件下,113个干预组只能找到221个对照。剩下的实在配不上。所以不能怪软件。怪松哥,找的这个案例不恰当,呵呵!

    图7B:匹配效果检验,可是没有值,为啥?因为这个Hansen检验只能针对1:1做,人家原则性杠杠滴,所以这里就没检验了.

图7结果(1)

 图7 C:不平衡检验,关键看前后对比,本例匹配前为0.914,匹配后为0.814,只下降一点点,不好不好!(该值取值0-1,越接近0,效果越好!)

图7D:既然效果不好,是那个因素再捣蛋呢?结果发现Propensity(咱们的评分)和married因素不好。他们的Std Mean Diff绝对值都大于0.25了。差评,差评!

3.2详细均衡性检验结果

    重点看图8最后一列,After绝对值越小越好,你也可以发现,确实是Propensity和Married出格了!

图8详细均衡性检验

3.3抖动散点图

图9可见两组匹配数据的概况。

图9 抖动散点图

3.4绝对标准化均数差异线图

图10中可见数据的收敛趋势!

图10 绝对标准化均数差异线图

3.5匹配前后效果直方图

    图11重点看右边红色框中,上下两组的形状越趋于一致,说明匹配效果越好!左边绿色框为匹配前的直方图。

图11前后匹配效果直方图

3.6Love Plot

    图12非常重要的一张Love Plot,发文章频频见到它!可以反映匹配前后均衡性改善情况!

图12 Love Plot

3.7配对后输出数据集

    配对再好,也要输出匹配后结果的数据集,然后才可以对该数据集进行后续干预因素效果分析。SPSS数据的数据集再Stata和R中,最简单干净利索的!

图13输出的配对后的数据集


整理不易,欢迎点亮再看!

---SPSS实战与疑难杂症系列---

【2013.】SPSS中如何设置哑变量

【2012.】分层回归之SPSS实践

【2011.】SPSS直接物理拆分为两个文件

【2010.】倾向性评分(PSM)的SPSS24.0实战

【2009.】如何SPSS实现随机分组

【2008.】趋势性卡方检验专题讨论

【2007.】2个因素均为重复测量数据方差分析

【2006.】为什么配对样本t检验,SPSS结果默认输出相关性

【2005.】SPSS软件作图如何去除顶部与右边边框

【2004.】如何让SPSS自动按照四分位数间距分组某变量【技巧】

【2003.】SPSS两独立样本t检验的简易算法(无需原始数据)

【2002.】太伤自尊了,自以为学会了,现在用SPSS连两独立样本t检验都做不好

【2001.】SPSS作图坐标轴不过原点的解决办法

【2000.】临床预测模型SPSS验证实战

------------------------------


    转藏 分享 献花(0

    0条评论

    发表

    请遵守用户 评论公约

    类似文章 更多