运用stata软件做meta分析
特此说明:本文为作者参考众多文献资料,汇总而成的经验成果,期间付出很多努力。所以,为尊重劳动成果,本文只给出各项分析中用到的命令语句,具体运算结果请战友们下载附件后自行参考!公平标价,童叟无欺。如对meta分析感兴趣,想进一步研究和探讨,可参考文后的相关链接。
案例+固定效应+亚组分析+回归分析+发表偏倚+漏斗图对称性检验+参考
1、案例
实验数据:本文以社会科学中最常用的相关系数为例,而不是以循证医学的OR,RR,RD值为效应值的计算依据。
authors | Correlation | n | Z | V | SE | abc | year |
| | ||||||
Fonda | 0.500 | 40 | 0.549 | 0.027 | 0.164 | 1 | 1993 |
Newman | 0.600 | 90 | 0.693 | 0.011 | 0.107 | 2 | 1990 |
Grant | 0.400 | 25 | 0.424 | 0.045 | 0.213 | 1 | 1998 |
Granger | 0.200 | 400 | 0.203 | 0.003 | 0.050 | 1 | 2007 |
Milland | 0.700 | 60 | 0.867 | 0.018 | 0.132 | 2 | 2006 |
Finch | 0.450 | 50 | 0.485 | 0.021 | 0.146 | 1 | 2003 |
2、固定效应模型的效应值合并:
- metan _ES _LCI_UCI, fixed lcols(authors year) favours(abc) texts(120)
随机效应模型只需将fixed改为random
备注:如果没有_ES,_LCI,_UCI值,需首先使用命令gen生成,具体步骤如下:
- gen _ES=zgen _LCI=z-sqrt(1/n-3)*1.96
- gen _UCI= z+sqrt(1/n-3)*1.96
注意空格和大小写。
3、亚组分析:
- metan z _LCI _UCI, sortby(year) lcols(authorsyear) textsize(100) by(abc) random
4、回归分析(将abc变量视为调节变量)
- gen _ES=z
- gen_seES=1/sqrt(n-3)
- metareg _ES,wsse(_seES) eform
- metareg _ES abc, wsse(_seES) reml
5、发表偏倚(漏斗图法)
- metafunnel _ES _seES, xtitle(ES) ytitle(StandardError of ES) egger
6、漏斗图对称性检验(egger法):
- metabias _ES _seES, egger gr
失安全数
直接参照公式手动计算更简单。
7、参考
1,【系统评价/Meta分析理论与实践】PDF下载,参考书,电子书,学习资料https://bbs.pinggu.org/thread-3541797-1-1.html
2,【Meta分析】Meta分析导论,电子版,PDF格式,开放下载,文献综述宝典https://bbs.pinggu.org/thread-3538643-1-1.html
3,[学科前沿] Meta-Analysis 書籍文章(Sage、ORM..等)及最新範例https://bbs.pinggu.org/thread-1030470-1-1.html