library(forestplot)
fp<-read.csv('/Users/Administrator/Desktop/forestplot1.csv',header=T)
forestplot(labeltext=as.matrix(fp[,1:3]),
lower=fp$OR_1,
upper=fp$OR_2,
mean=fp$OR_mean,
zero=1,
boxsize=0.4,
xticks=c(0,1,2),
lineheight = unit(7,'mm'),
colgap=unit(2,'mm'),
lwd.zero=2,
lwd.ci=0.1,
col=fpColors(box='#458B00',
summary='#8B008B',
lines = 'black',
zero = '#7AC5CD'),
xlab="OR",
lwd.xaxis =2,
txt_gp = fpTxtGp(ticks = gpar(cex = 0.85),
xlab = gpar(cex = 0.8),
cex = 0.9),
lty.ci = 1,
title = "Forestplot",
line.margin = 0.08,
graph.pos=2)
mydata <- read.csv('/Users/Administrator/Desktop/forestplot1.csv')//读入模型
mylogit <- glm(admit ~ Age+ Cr + Bleeding+RBCs+Plasma+OT, data = mydata, family = "binomial")//拟合模型
exp(cbind(OR = coef(mylogit), confint(mylogit)))//查看模型的OR值和置信区间
plot_model(mylogit)//得到结果
中间有个逻辑有问题