运用stata软件做meta分析

特此说明：本文为作者参考众多文献资料，汇总而成的经验成果，期间付出很多努力。所以，为尊重劳动成果，本文只给出各项分析中用到的命令语句，具体运算结果请战友们下载附件后自行参考！公平标价，童叟无欺。如对meta分析感兴趣，想进一步研究和探讨，可参考文后的相关链接。


实验数据：本文以社会科学中最常用的相关系数为例，而不是以循证医学的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

固定效应模型的效应值合并：

metan _ES _LCI_UCI, fixed lcols(authors year) favours(abc) texts(120)

随机效应模型只需将fixed改为random

备注：如果没有_ES，_LCI，_UCI值，需首先使用命令gen生成，具体步骤如下：

gen _ES=z

gen _LCI=z-sqrt(1/n-3)*1.96

gen _UCI= z+sqrt(1/n-3)*1.96

注意空格和大小写。

亚组分析：

metan z _LCI _UCI, sortby(year) lcols(authorsyear) textsize(100) by(abc) random

回归分析（将abc变量视为调节变量）

gen _ES=z

gen_seES=1/sqrt(n-3)

metareg _ES,wsse(_seES) eform

metareg _ES abc, wsse(_seES) reml

发表偏倚（漏斗图法）

metafunnel _ES _seES, xtitle(ES) ytitle(StandardError of ES) egger

漏斗图对称性检验（egger法）：

metabias _ES _seES, egger gr

失安全数

直接参照公式手动计算更简单。

参考资料：

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