stata 16安装与meta简单应用
一、安装 根据BV号:
1LT411R7Ws
获取资源及安装，注意安装所在文件夹不能含中文，不安装在c盘也可以。 根据BV号:
1Dj411G71Y
下载外部命令包，名称为plus的文件夹。 将plus文件夹放置在
c盘
如下路径
删除该plus文件夹m文件里的名称为
metabias
的文件（应该只有一个），如图
删除后打开stata 16，电脑需要联网，运行如下命令
安装完成后根据bv号:
1cQ4y1Q7i2
制作profile.do文件，注意复制代码时窗口需最大化，否则运行失败。成功后可见如下
二、meta分析（干预性实验） 1.森林图绘制 1.1连续性变量（什么是连续性可以百度） 1.1.1所需数据：干预组人数n1，干预组数据平均值mean1，干预组标准差sd1，对照组人数n2，对照组数据平均值mean2，对照组标准差sd2 最好制作为表格，方便复制
打开stata 16，打开数据编辑器，路径如图
将表格数据粘贴在每列第一个空格，若跳出窗口均选数据
输入完成后可以单击2次var更改名称
为方便自己认得出。更改完成后如图
回到命令界面。输入命令db metan运行，或如图
运行后
count 二分类数据，continuous连续性数据，effect/ci 效应值/置信区间，effect/se 效应值/标准误。我们选第二个 下拉菜单后单击填入
点击该窗口第三个按钮，如图，根据文献中类型选择模型，随机实验一般为random和nostandard。选择后确定。
1.2二分类数据 所需数据:a干预组发生事件数，b干预组未发生事件数（
注意！不是总人数
）
，c对照组发生事件数，d对照组未发生事件数（
同b
）。同样运行db metan命令。选择数据种类为count，下拉菜单点击abcd，顺序不要错
根据文献内容点击该窗口第二个按钮，选择模型和统计量，选择后确定。
1.3 EFFECT/CI适用于文献仅有相对危险度rr及置信区间（本小节未亲自试验，存疑） 所需数据，rr，最低置信区间lowerCI（lcl），最高置信区间upperCI（ucl），同样db metan依次选择rr，lcl，ucl。
在该窗口第四个按钮选择模型，效应值标签为rr（你设rr叫什么就是什么）。设定后运行。
1.4已知相关系数r值。本小节采取文献《中国老年人睡眠质量与抑郁、焦虑相关性的meta分析》 已知相关系数r，本文献采取z值计算。 z=0.5*ln((1+r)/(1-r）），z值方差Vz=1/（n样本量-3），z值标准误Se=Vz＾0.5。网上搜的公式，用就好了，我也不知道为什么。在office或wps表格中先计算出vz和se
type选择effect/se，下拉菜单依次选z，se。第四个按钮中标签label改为z值。点击确定
绘制后。
2发表偏倚 我用stata就为了他的发表偏倚，stata16画完begg它同时就把egger画出来了，很方便。
画发表偏倚前按1步骤画森林图
画完森林图后
2.1连续性变量
输入命令:metabias6空格_ES空格_seES, 空格graph（begg）后运行，此处写空格是为了便于讲解。stata16中做发表偏倚
metabias后必须有个6，逗号后的空格不要忽略。标点全英文
metabias6 _ES _seES, graph(begg) 2.2二分类变量 gen命令可以在数据编辑器里新建一列你所需要的数据。 gen logES=log(_ES)回车 metabias6 logES _selog(ES), graph(begg)回车 2.3EFFECT/CI 据网络所述，需要先定义三个数据。每个命令回车一次。 gen logrr=log(rr) gen loglcl=log(lcl) gen logucl=log(ucl) metabias6 logrr loglcl logucl, graph(begg) 或者不定义，我没试过 metabias6 log(rr) log(lcl) log(ucl), graph(begg) 2.4EFFECT/SE 运行命令metabias6 z se, graph(begg) 依据上述文献如图
书写格式，写begg法的z和p，egger法的t和p
3，亚组分析 如本文献，建议在表格中先分组
同样根据文献类型输入数据
同样用绘制森林图那个窗口，下方有个by variable，选择你的亚组分组依据
运行
如何统计亚组的发表偏倚？ 此例中，运行 metabias6 z se, by(year2) graph(begg) 所以说，连续性变量啊，二分类变量啊，在graph前面加 by(你分组的自定义标签)就可以了。只不过需要在excel中提前分好组。stata好像一次只能根据一个亚组分析，就是by variable下拉标签里你选两个就报错了。 三随便说说 revman也是不错的软件，绘制森林图很方便，可是不如stata编辑性强一点，可以随便加文本。revman的发表偏倚也很方便，只不过出来的是图，没有begg和egger这种具体的数字。改天我要继续学习一下revman的亚组分析，stata的楼斗图也还没学会。 感谢上述列举视频中b友的资源，愿看到这篇文章的你前程似锦，sci随便发👍。
