R语言广义线性模型Logistic回归模型亚组分析及森林图绘制

#Logistic回归案例 6 亚组分析森林图

library(forestplot)

rs_forest <- read.csv( s_forest.csv,header = FALSE)

# 读入数据的时候大家一定要把header设置成FALSE，保证第一行不被当作列名称。

# tiff(Figure 1.tiff,height = 1600,width = 2400,res= 300)

forestplot(labeltext = as.matrix(rs_forest[,1:3]),

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

mean = rs_forest$V4, #设置均值

lower = rs_forest$V5, #设置均值的lowlimits限

upper = rs_forest$V6, #设置均值的uplimits限

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

#该参数接受一个逻辑向量，用于定义数据中的每一行是否是汇总值，若是，则在对应位置设置为TRUE，若否，则设置为FALSE；设置为TRUE的行则以粗体出现

zero = 1, #设置参照值，此处