R语言是处理生信的神器，处理的图形非常高端洋气上档次。但是敲代码着实麻烦。看完这篇文章，你只要按照笔者的内容复制粘贴代码即可，也可以根据需要更换一点简单的代码，要是不愿意更换，把你的excel表头改成和笔者一样的，直接复制即可。

　　（一）使用R语言进行单量率的meta分析

　　1软件下载

　　R语言新版本R version 3.6.1(截止到2019年12月20日)，下载地址https://www.r-project.org/。使用的镜像大陆的镜像，要是发现下载不流畅，选择另一个镜像。注意，镜像选错，就下载不了!

　　Rstudio是从官网下载的新版本，下载地址https://www.rstudio.com/。

　　2准备数据库和Packages

　　将数据库格式改为.csv格式。.xls文件也能在R语言下运行，但是.csv文件更小，在R中运行效果更好。

　　打开文件可以执行代码：

　　library(readr)

　　test <- read_csv("C:/Users/Administrator/Desktop/test.csv")

　　View(test)

　　也可以使用菜单，打开菜单顺序：file-import dateset-fromtest(readr)

　　此时，截取的结果便在软件中展示出来：

image.png

　　导入数据后：

image.png

image.png

　　执行代码：strepto=read_csv("C:/Users/Administrator/Desktop/test.csv")

　　代码解读：

　　“C:/Users/Administrator/Desktop/test.csv”值得是文件的位置，“<-”相当于C语言中的“=”，是赋值作用。

　　打开使用的函数包

image.png

　　点击后，在命令框自动生成以下代码

image.png

　　注意：要是没有安装“matrix”的需要安装一下，执行以下代码即可：

　　> install.packages("matrix")

　　> head(strepto)

　　看看我们的图标情况，只能显示六行

image.png

　　下一步我们要做的就是对数据进行处理，执行以下代码：

　　>res=metabin(tn,tt,cn,ct,data =strepto,sm="OR",method="MH",studlab=paste(study),comb.random=TRUE)

　　> forest(res)

　　tn,tt,cn,ct可以根据自己的表头更换，comb.random=TRUE可以改为comb.fixed=TRUE(固定效应模型)

　　结果如图所示：

image.png

　　3异质性分析

　　首先来一个亚组分析

　　>res=metabin(tn,tt,cn,ct,data =strepto,sm="OR",method="MH",studlab=paste(study),byvar=group,comb.random=TRUE)

　　> forest(res)

image.png

　　漏斗图来一个

　　> funnel(res)

image.png

　　再来一个敏感分析

　　> metainf(res)

image.png