1
在R命令框输入以下命令：
install.packages("rmeta") #安装rmeta程序包
library(rmeta)
data(cochrane) #加载演示数据库
cochrane #显示数据库cochrane
2
计算fixed effect model 固定效应模型，输入以下命令：
model.FE <- meta.MH(n.trt,n.ctrl,ev.trt,ev.ctrl, names=name,data=cochrane)
summary(model.FE)
3
计算random effect model 随机效应模型，输入以下命令：
model.RE <- meta.DSL(n.trt,n.ctrl,ev.trt,ev.ctrl, names=name,data=cochrane)
summary(model.RE)
4
绘制森林图，输入以下命令：
tabletext<-cbind(c("","Study",model.FE$names,NA,"Summary"),c("Deaths","(steroid)",cochrane$ev.trt,NA,NA), c("Deaths","(placebo)", cochrane$ev.ctrl, NA,NA), c("","OR",format(exp(model.FE$logOR),digits=2),NA,format(exp(model.FE$logMH),digits=2)))
#上面的命令帮我们构建森林图需要的文字描述部分，包括1、name， 2、治疗组死亡人数，3、对照组死亡人数，4、单个研究的OR值。
m<- c(NA,NA,model.FE$logOR,NA,model.FE$logMH)
l<- m-c(NA,NA,model.FE$selogOR,NA,model.FE$selogMH)*2
u<- m+c(NA,NA,model.FE$selogOR,NA,model.FE$selogMH)*2
#上面三行命令是准备森林图需要的参数包括OR值和OR值95%CI区间的上下限
forestplot(tabletext,m,l,u,zero=0,is.summary=c(TRUE,TRUE,rep(FALSE,8),TRUE), clip=c(log(0.1),log(2.5)), xlog=TRUE, col=meta.colors(box="royalblue",line="darkblue", summary="royalblue"))
#上面的命令是最后绘制森林图
