数据
variable point low high
金湖(协变量借用) 148.40 98.74 198.05
金湖(地理借用) 134.66 85.93 183.40
灵璧(协变量借用) -25.76 -54.69 3.18
灵璧(地理借用) -7.65 -37.12 21.83
墉桥(协变量借用) 17.01 -30.51 64.54
墉桥(地理借用) -31.63 -78.27 15.01
盱眙(协变量借用) 23.99 4.07 43.90
盱眙(地理借用) -74.25 -94.78 -53.72
代码
getwd()
setwd("F:/statistics root/r workspace")
data <- read.csv("data.csv")
pacman::p_load(forestplot)
##提取所需数据
txt <- data[,1]
hr <- data[,c(2:4)]
##绘制森林图
forestplot(txt,hr,
ci.vertices = TRUE, # 穿过box的线段两端是否带下垂线,默认FALSE。
ci.vertices.height = 0.1, #可是设置线段两端垂线的长度
graph.pos = 2, #控制图形输出在第几列
graphwidth = unit(100,"mm"), #横轴长度
clip=c(-100,200), #设置横轴范围,上一条语句必须有本语句,x轴范围会依照数据显示
xticks = c(-100,-50,0,50,100,150,200), #设置x轴显示的坐标,包括显示范围
txt_gp = fpTxtGp(ticks = gpar(cex = 1),xlab = gpar(cex = 1.2), cex = 1.2),
#cex为文本字体大小,ticks为坐标轴大小,xlab为坐标轴文字字体大小
col=fpColors(box="black",line="black",summary="black"), #图形颜色
fn.ci_norm = fpDrawCircleCI, #置信区间中间图形状
boxsize = 0.1, #置信区间中间图形大小
vertices = TRUE, #为端点加垂线
xlab="预测死亡率-真实死亡率(/10万)" )
示图