来源:微信公众号小侃数据
大家好,本教程将介绍如何使用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=”