森林图的绘制

res=read.table('PUD_PTSD_MR.txt',header = T,sep = '\t')
res= generate_odds_ratios(res)
write.table(res,file='PUD_PTSD.txt',sep = '\t',row.names = F,quote = F)
library(forestplot)
inputFile="PUD_PTSD.txt"
outFile="PUD_PTSD_forest.pdf"
rt=read.table(inputFile,header=T,sep="\t",row.names=5,check.names=F)
gene=rownames(rt)#读取基因列
hr=sprintf("%.3f",rt$"or")#获取HR列取小数点后3位
hrLow=sprintf("%.3f",rt$"or_lci95")#获取95%置信区间取HR.95L列小数点后3位
hrHigh=sprintf("%.3f",rt$"or_uci95")#获取95%置信区间HR.95H列小数点后3位
pVal=ifelse(rt$pval<0.001, "<0.001", sprintf("%.3f", rt$pval)) #获取P值
Hazard.ratio=paste0(hr,"(",hrLow,"-",hrHigh,")")#合并为一个数据集
#输出图形
pdf(file=outFile, width = 8, height =4.5)
n=nrow(rt)
nRow=n+1
ylim=c(1,nRow)
layout(matrix(c(1,2),nc=2),width=c(4,2))
#森林图左边的基因信息
xlim = c(0,3)
par(mar=c(4,2,1.5,1.5))
plot(1,xlim=xlim,ylim=ylim,type="n",axes=F,xlab="",ylab="")#axes=F表示禁止生成坐标轴
text.cex=0.8 #放大0.8倍
text(0,n:1