利用forestplot包绘制效应量亚组分析的类森林图

在临床研究中，为了验证主要结论在不同特征人群中的一致性，常需要在完成主效应分析后，再进行亚组分析，如发表在《Cardiovascular Diabetology》杂志上的一篇论文[1], 作者在得到TyG与动脉硬化发生风险的OR后，对性别、年龄、糖化、病程等6个变量进行了亚组分析，为了更加直观的展示各亚组的OR及95%CI，作者采用类森林图进行展示（图1）。类似的，图1中的P
值也可以换为该指标与待分析指标的交互作用P值（图2）[2]。

绘制此类图，常用的是forestplot包，在绘制图之前，需要将单因素或多因素回归得到的OR值以及95%置信区间整理成如图3所示的表格。为了方便，可将带有数值的6列分别命名为v1-v6。将表格命名为plot，并放置于R项目所在路径下。

1. 首先读取原始数据文件

plot<-read.csv("plot1.csv"