贝叶斯软件genle教程_手把手教你用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("D:R网状教程") #你可以更改为自己的工作目录，在后面环境变量增多的时候，在这一步前可以加一个rm(list=ls())清除环境变量。

然后敲入：

library(gemtc) #装载gemtc包

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

其中每一行代表研究的一个臂，study代表研究，treatment代表干预名称，mean代表连续变量的均值，std.dev代表连续变量的标准差，sampleSize代表每一个研究中每一个臂的样本量。注意，以上英文单词必须准确无误，大小写都要一致，不然等会分析无法进行。

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

把刚刚准备好的数据导入到R里，敲入：

data <- read.csv("data.csv", sep="