森林图-R语言代码-(含备注-附练习数据).pdfVIP

森林图R语言代码

（含备注，附练习数据）

#代码

library(forestplot)

library(readr)

#从本地导入forest.csv数据，命名成mydata

格式的文件路径径例如：mydata-mydata-

#定义表头

of.of.

#增加一行

mydata-rbind(rep(NA,3),mydata)

#查看数据

View(mydata)

forestplot(labeltext=row_name,

#设置用于文本展示的列，此处我们用数据的前三列作为文本，在图中展示

zero=1,#设置参照值，此处我们展示的是OR值，故参照值是1

xlog=T,#对数x轴

xticks=c(0.5,1.0,2.0,4.0),#x轴刻度

boxsize=0.1,#设置点估计的方形大小

lineheight=unit(10,mm),#设置图形中的行距

colgap=unit(3,mm),#设置图形中的列间距

lwd.zero=2,#设置参考线的粗细

lwd.ci=1.5,#设置区间估计线的粗细

#使用fpColors()函数定义图形元素的颜色，从左至右分别对应点估计方形，汇总

值，区间估计线，参考线

设置x轴标签

lwd.xaxis=2,#设置x轴粗细

graph.pos=3,#设置森林图的位置，此处设置为3，则出现在第三列

is.summary=c(T,F,T,F,F,T,F,F,T,F,F))#该参数接受一个逻辑向量，用于定义数据中

的每一行是否汇总值，若是，则在对应位置设置TRUE,若否，则设置为FALSE,

设置为TRUE的行则以粗体出现

library(ggplot2)

forestplot(labeltext=row_name,

zero=1,

xlog=F,

xticks=c(0.5,1.0,2.0,4.0),

boxsize=0.1,

lineheight=unit(10,mm),

colgap=unit(3,mm),

lwd.zero=2,

lwd.ci=1.5,

lwd.xaxis=2,

graph.pos=3,

坐标轴

轴

is.summary=c(T,F,T,F,F,T,F,F,T,F,F))

附件：练习数据

SubgroupsNo.ofpatientsOR(95%CI)ORLLUL

总人群15631.79（1.41,2.28）1.791.412.28

性别

男性7691.81（1.36,2.39）1.811.362.39

女性7941.79（1.18,2.79）1.791.182.79

年龄分组

老年人10561.86（1.43,2.42）1.861.432.42

青年人5071.71（0.63,5.17）