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

Posted 2021-09-04 Data+Science+Insight

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

#Logistic回归案例 6 亚组分析森林图
library(forestplot)
rs_forest <- read.csv(\'rs_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, #设置参照值，此处