这是一段使用forestplot包画森林图的一段程序,已测能跑,需者自取~
tips1: 首先需要安装和加载forestplot包,如果通过常规途径无法加载,请检查一下自己是否有装Rtools
install.packages("forestplot")#安装forestplot包
library(forestplot)#加载forestplot包
tips2: 用于画森林的数据列(lower,or,upper)需要是数字形式,当无法跑出图片时,可以检查一下是不是这三列里面掺进了文本(笔者的血泪教训)
tips3: 以下代码只是forestplot的一小部分,但应该够用了(之后学习了会再加的)
senlin <- read.csv("D:\\analysis learn\\forest2.csv",header=T,sep=",")#导入数据
labeltext <- as.matrix(senlin[,1])#将数据框的第一列转换成矩阵
attach(senlin)#绑定数据框
library(forestplot)
windowsFonts(myFont1 = windowsFont("Times New Roman"))#调用Times New Roman字体
forestplot(labeltext,#图形文本部分
mean=HR,#图形HR/OR部分
lower=LowerCI,#95%CI下限
upper=UpperCI,#95%CI上限
is.summary=c(T,F,F,T,F,F,T,F,F,T,F,F,T,F,F,T,F,F,F,F),#指定突出显示和非突出显示
align=c("l"),#文字对齐方式,l,r,c,分别表示左对齐、右对齐、居中对齐
graph.pos=2,#表示图形元素在第二列(也可设置为left 或 right)
hral_lines = T,
clip=c(0,1.4),
zero=1,#设置无效线的横坐标
lwd.zero=1,#设置无效线的线条宽度
lty.zero="solid",#这个好像是设置线条类型,但是我不知道虚线用什么表示
title="森林图",
xlab="Risk ratio(95%CI)",#设置x轴的标签
lwd.xaxis=1,#设置x轴的宽度
xticks=c(0,0.2,0.4,0.6,0.8,1.0,1.2,1.4),#设置x轴的刻度
lty.ci=1,#设置置信区间的线条类型
lwd.ci=1,#设置置信区间的线条宽度
ci.vertices=T,#两端是否有竖线
ci.vertices.height=0.2,#两端竖线的高度
boxsize=0.2,#设置箱体的大小
col=fpColors(
box = "red",#箱体的颜色为黑色
lines = "orange",#置信区间端点线条的颜色
summary = "yellow",#文字的颜色
zero = "green",#无效线的颜色
text = "purple",#文本的颜色
axes = "black",#x轴的颜色
hrz_lines = "black"))#设置线条的颜色