Meta Analysis in Stata17

本文作者：熊保维，中南财经政法大学经济学院

本文编辑：王玉洁

技术总编：方一卓

Stata&Python云端课程来啦！

       为了感谢大家长久以来的支持和信任，爬虫俱乐部为大家送福利啦！！！Stata&Python特惠课程双双上线腾讯课堂～爬虫俱乐部推出了Python编程培训课程、Stata基础课程、Stata进阶课程、Stata文本分析、正则表达式、网络爬虫、基本字符串课程。报名课程即可加入答疑群，对报名有任何疑问欢迎在公众号后台和腾讯课堂留言哦。我们在这篇推文的最后提供了每门课程的课程二维码，大家有需要的话可以直接扫描二维码查看课程详情并进行购买哦～

引言

导读

相信不少同学已经使用了一段时间Stata17，不论是DID的官方命令还是Stata17整体的提速，每一个感受过Stata17魅力的人都对其拍手称赞！

今天小编就带大家认识一下Stata17中进一步完善的一个板块——Meta分析。

1. Meta分析的基本信息

1.1 Meta是什么:

• Meta意为较晚出现的更为综合的事物，它的中文名称有元分析、荟萃分析、数据再分析等，意味着对他人分析结果的再分析。20世纪初以来Meta分析成为医学和心理学等实验研究的常见方法，近几十年在经济领域研究中也得到了广泛应用。

• 它的原理就是对多个同类独立定量研究的结果进行汇总分析，达到增加样本含量，提高检验效能的目的，尤其是当多个研究结果不一致时，采用Meta分析可得到更加接近真实情况的结果。

1.2 为什么要用Meta:

• 弥补叙述型定性文献综述的缺陷

• 增强文献研究的客观性和科学性

1.3 可用软件：

• Revman|  Stata | R ...

2. Meta分析的类型

1. 常规Meta分析

2. 单组率的Meta分析

3. Meta回归分析

4. 累积Meta分析

5. 网状Meta分析

6. 诊断性Meta分析

7. ……

不同分析类型有不同的分析方法和质量评价方法，详细介绍可参考：https://zhuanlan.zhihu.com/p/40336871

3. Meta分析的三种模型

当有K组研究，研究模型产生的估计系数为θ时，可考虑选用以下三种模型之一：

1. Random-Effects (RE)即每一次估计的结果是具有随机性的。这是Stata对大多数元分析的默认模式。mean of the distribution of θj = θ + uj

2. Fixed-Effects (FE)如果假设这些研究有不同的真实效应量，并且你只对提供关于这些具体研究的推论感兴趣，那么FE模型是合适的。weighted average of the K true study effects

3. Common-Effects (CE)如果研究同质性的假设对您的数据是合理的，那么可以考虑使用CE模型。common effect (θ1 = θ2 = · · · = θK = θ)

4. Meta分析的步骤

1. 选题

2. 文献检索（如何尽可能保证检索出所有相关文献的同时缩小检索文献的数目）

3. 文献筛选（根据文献纳入排除标准筛选出合适文献筛选出）

4. 数据提取（提取出需要合并的数据）

5. 质量评价（利用量表对文献质量进行讨论）

6. 异质性检验（检验各项研究是否同质）

7. 数据合并（利用软件对数据进行合并）

8. 对结果的解释（关键）

5. Stata17的Meta命令汇总

help meta 可得上述界面。详细还可help具体命令。

6. Stata17新增的Meta命令

6.1 多元变量分析 Multivariate meta-analysis

• 新增命令：meta mvregress

• 将多个研究的结果结合在一起时，其中的每一个研究可能报告了多种效应规模(multiple effect sizes)，这些效应之间很可能互相关联。之前的Meta命令得出的结果会忽略相关性；
• 先可使用meta mvregress命令执行多元Meta分析，用于解释相关性。

多元变量分析的语法规则：

多元变量分析的官方示例：

*Treatment of moderate periodontal disease
*牙周疾病治疗的数据

clear all
use https://www.stata-press.com/data/r17/periodontal   //使用网络数据
describe                                              //留意数据形式
quietly meta set y1 s1, studylabel(trial) eslabel("Mean diff.")
meta forestplot, esrefline               //单变量分析
quietly meta set y2 s2, studylabel(trial) eslabel("Mean diff.")
meta forestplot, esrefline                     //单变量分析
meta mvregress y1 y2, wcovvariables(v11 v12 v22)      //多元变量分析

多元变量分析的示例输出结果：

6.2 加尔布雷斯图 Galbraith plots

• 新增命令：meta galbraithplot

• 该命令可以生成Galbraith图用于Meta分析，这类图可以用于评估研究的异质性(assessing heterogeneity of the studies)和检测潜在的极端值(potential outliners)；

• 研究数量较多时会被用来替代森林图，用于总结Meta分析结果。

  加尔布雷斯图的语法规则：

加尔布雷斯图的命令示例1：

clear all
use https://www.stata-press.com/data/r17/bcgset  //使用网络数据
meta galbraithplot                               //制作加尔布雷斯图

命令示例1输出图片：

加尔布雷斯图的命令示例2：

更多图片样式，可自行探索~

6.3 留一元分析 Leave-one-out meta-analysis

留一元分析的语法规则：作为后缀选项

• 加在summarize后的示例:

clear all
use https://www.stata-press.com/data/r17/pupiliqset
keep in 1/10
meta query, short
meta summarize, leaveoneout sort(se)

• 加在summarize后的输出结果:

• 加在森林图后的示例:

clear all
use https://www.stata-press.com/data/r17/bcgset
meta forestplot, leaveoneout rr

• 加在森林图后的输出结果：

以上就是关于Stata17中Meta Analysis的相关内容了，如果大家对Stata17的其他内容感兴趣的话，可以查阅我们之前发布的推文来进行详细的了解。

最后，我们为大家揭秘雪球网(https://xueqiu.com/)最新所展示的沪深证券和港股关注人数增长Top10。

腾讯课堂课程二维码

 对我们的推文累计打赏超过1000元，我们即可给您开具发票，发票类别为“咨询费”。用心做事，不负您的支持！

往期推文推荐                cngdf——名义GDP与实际GDP之间的摆渡船
         合并，“纵”享新丝滑：frameappend & xframeappend
什么是全局配置项？|从零开始的Pyecharts（二）帮你拿下数据可视化|从零开始的Pyecharts 这十年，《金融研究》的编委和读者偏爱哪些研究话题和文章?
学好这一手，英语词典常在手 用spmap看中国空气质量
OpenCV库——轻松更换证件照背景色

关于我们 

   微信公众号“Stata and Python数据分析”分享实用的Stata、Python等软件的数据处理知识，欢迎转载、打赏。我们是由李春涛教授领导下的研究生及本科生组成的大数据处理和分析团队。

   武汉字符串数据科技有限公司一直为广大用户提供数据采集和分析的服务工作，如果您有这方面的需求，请发邮件到[e-mail]，或者直接联系我们的数据中台总工程司海涛先生，电话：[phone]，wechat: super4ht。海涛先生曾长期在香港大学从事研究工作，现为知名985大学的博士生，爬虫俱乐部网络爬虫技术和正则表达式的课程负责人。

此外，欢迎大家踊跃投稿，介绍一些关于Stata和Python的数据处理和分析技巧。

投稿邮箱：[e-mail]投稿要求：
1）必须原创，禁止抄袭；
2）必须准确，详细，有例子，有截图；
注意事项：
1）所有投稿都会经过本公众号运营团队成员的审核，审核通过才可录用，一经录用，会在该推文里为

作者署名，并有赏金分成。

2）邮件请注明投稿，邮件名称为“投稿+推文名称”。
3）应广大读者要求，现开通有偿问答服务，如果大家遇到有关数据处理、分析等问题，可