stata 亚组分析_手把手教你用R的gemtc包对分类变量进行贝叶斯网状Meta分析

本教程将介绍如何使用R的gemtc包对分类变量进行贝叶斯网状Meta分析。

前提条件：

需要下载R软件，本教程使用的R版本为3.5.3，URL：

https://cloud.r-project.org/bin/windows/base/old/3.5.1/R-3.5.1-win.exe

以及RStudio，一个R的友好交互界面，URL：

https://download1.rstudio.org/desktop/windows/RStudio-1.2.1335.exe

因为gemtc包现在只支持调用JAGS，所以还需要下载JAGS的exe文件

https://nchc.dl.sourceforge.net/project/mcmc-jags/JAGS/4.x/Windows/JAGS-4.3.0.exe

以上文件下载后直接双击安装即可。

准备工作完成后，打开RStudio后，在console里输入:

install.packages("gemtc")

第一次安装的时候，安装速度可能比较慢，耐心等待一下。

 等安装完成后，设置我们的工作目录，后面可以把待分析的数据放入这个工作目录，敲入：

setwd("G:\\R网状教程\\分类变量")  #你可以更改为自己的工作目录，在后面环境变量增多的时候，在这一步前可以加一个rm(list=ls())清除环境变量。

然后敲入：

library(gemtc)     #装载gemtc包

 将我们的数据准备成如下格式(与stata做网状的格式相同)：

其中每一行代表研究的一个臂，study代表研究，treatment代表干预名称，responders代表每一个臂的反应人数，sampleSize代表每一个研究中每一个臂的样本量。注意，以上英文单词必须准确无误，大小写都要一致，不然等会分析无法进行。

大家可以在excel里准备要分析的数据，然后另存为csv文件(以逗号为分隔符的表格文件)，在这里，我把数据文件存为data.csv(你可以更改为自己想要的名字)。然后把这个文件放到我们的工作目录之下。

把刚刚准备好的数据导入到R里࿰