贝叶斯软件genle教程_手把手教你用R的gemtc包对生存数据进行贝叶斯网状Meta分析...

大家好，本教程将介绍如何使用R的gemtc包对生存数据(HR为效应量)进行贝叶斯网状Meta分析。

前提条件

需要下载R软件(推荐使用的R版本为3.5.3)，以及RStudio(一个R的友好交互界面)。

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

install.packages(“gemtc”)

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

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

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

然后敲入：

library(gemtc)     #装载gemtc包

关于如何将生存数据准备成可以用R的GEMTC分析的格式，以及如何计算多臂研究中参考臂的标准误，请参考如下教程：

生存数据的网状Meta分析||如何整理为可以用R的gemtc包进行分析的格式

主要分析

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

data "sample1.csv", sep=",", header=T)   #sep代表逗号为分隔符，header=T代表把首行的数据作为标题，

把数据弄成gemtc的网络格式，敲入：

network #这里

输入：

plot(network)

可以看到网状图如下：

当然了，这些大小、颜色、字体等等的参数都是可以调整的&#x