贝叶斯软件genle教程_手把手教你用R的gemtc包对连续变量进行贝叶斯网状Meta分析…

  • Post author:
  • Post category:其他


96bdd6163d22c9b0ec02edf74b16b168.png

来源:微信公众号小侃数据

大家好,本教程将介绍如何使用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做网状的格式相同):

aa27fe8eca6804237f65a4642dc44181.png

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

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

把刚刚准备好的数据导入到R里,敲入:

data <- read.csv(“data.csv”, sep=”



版权声明:本文为weixin_33379297原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。