R语⾔绘制森林图-forestplot ####森林图绘制 forestplot
#### A practice for forest plot
##t your working directory
workdir <- "D:/R"
twd(workdir)
datafile <- file.path(workdir,"ForestPlotData.csv")微信消息延迟
data <- read.csv(datafile, stringsAsFactors=FALSE)
data[1:5,1:5]
## 定义亚组
subgps <- c(4,5,8,9,12,13,16,17,20,21,24,25,28,29,32,33)
data$Variable[subgps] <- paste(" ",data$Variable[subgps])
## Combine the count and percent column，以括号相连
np <- ifel(!is.na(data$Count), paste(data$Count," (",data$Percent,")",p=""), NA)
## The rest of the columns in the table.
tabletext <- cbind(c("Subgroup","\n",data$Variable),
c("No. of Patients (%)","\n",np),
c("4-Yr Cum. Event Rate\n PCI","\n",data$PCI.Group),
c("4-Yr Cum. Event Rate\n Medical Therapy","\n",data$Medical.Therapy.Group),
红色男爵c("P Value","\n",data$P.Value))
##
library(forestplot)
png(filename = "Forestplot.png",width=960, height=640)怎样做炸鸡腿
forestplot(labeltext=tabletext, graph.pos=3,
mean=c(NA,NA,data$Point.Estimate),
座位的拼音lower=c(NA,NA,data$Low), upper=c(NA,NA,data$High),
title="Hazard Ratio",
##定义x轴我的心儿怦怦跳作文400字
xlab=" <---PCI Better--- ---Medical Therapy Better--->",
##设置透明条位置返聘协议
hrzl_lines=list("3" = gpar(lwd=1, col="#99999922"),
"7" = gpar(lwd=60, lineend="butt", columns=c(2:6), col="#99999922"),
"15" = gpar(lwd=60, lineend="butt", columns=c(2:6), col="#99999922"),
"23" = gpar(lwd=60, lineend="butt", columns=c(2:6), col="#99999922"),
白菜猪肉
"31" = gpar(lwd=60, lineend="butt", columns=c(2:6), col="#99999922")),好歌词
txt_gp=fpTxtGp(label=gpar(cex=1.25),
