Stata是一款非常强大的统计和作图软件,本文中用的版本是Stata14.2。网上目前存在的教程多半是用mvmeta包来做网状Meta,但mvmeta相对繁琐,需要打不少的代码去转换数据格式,Network包本质上是mvmeta,不过对用户更加友好,只需要执行简单的代码就可以数据处理。本文将介绍用Stata的Network包实现分类变量的网状Meta分析。
1. 假设我们已经安装好了Stata,在命令栏敲入:net from http://www.mtm.uoi.gr,然后点击network_graphs,下载network包。
2. 将要处理的数据在Excel中整理成下图所示:
![d426d99254e671ee8969410826978d45.png](https://img-blog.csdnimg.cn/img_convert/d426d99254e671ee8969410826978d45.png)
其中id代表研究,t代表治疗,r是反应数(本例中是每组因副作用脱落的数目),n是样本量;注意每一行代表一个研究的一个臂(组),如果一个研究有三臂,那么就占三行。分析时需将id以数字来依次编号Study。同样,依次以数字来编号t。特别注意的是,同一研究中从上到下,治疗的编码只能增大,不能降低。一般来说,将安慰剂编码为1。
3. 点击Stata菜单栏Data->Data Editor-> Data Editor(Edit),将整理好的数据直接拷贝到里面,注意拷贝的时候选择首行为变量名;Excel中第1列Study不用拷贝;如下图所示:
![159711bc91e84191acd45e3c21522f95.png](https://img-blog.csdnimg.cn/img_convert/159711bc91e84191acd45e3c21522f95.png)
4. 返回Stata的首界面,在命令栏敲入:
network setup r n, studyvar(id) trtvar(t) format(augment) or
以上代码效应量为OR,如果是用RR作为效应量,则敲入:
network setup r n, studyvar(id) trtvar(t) format(augment) rr
结果如下图所示:
![c23d647c7e58b73c65c18ae0fa82f729.png](https://img-blog.csdnimg.cn/img_convert/c23d647c7e58b73c65c18ae0fa82f729.png)
继续点击Stata菜单栏Data->Data Editor-> Data Editor(Edit),发现数据变成如下模样: