单臂,单个率的meta分析我们经常使用到,因为缺乏对照组,我们不能使用revman进行处理,下面小编介绍,如何使用stata进行单个率的meta分析。

  01输入数据

  我们使用下面的数据进行计算,其中a是时间数、n是总数、r是比例。

stata:如何进行单个率的meta分析

  首先在stata中输入我们的数据

stata:如何进行单个率的meta分析

  Study的格式是这样的

stata:如何进行单个率的meta分析

  其他的格式是这样的

stata:如何进行单个率的meta分析

  下面执行代码:. generate ser=sqrt(r*(1-r)/n),这一步是进行数据转换的。

  02制作森林图

  执行代码metan r ser, random label(namevar=study)。

  固定效应模型是这个:metan r ser, fixed label(namevar=study)。

  我们执行随机效应模型,结果如下图:

stata:如何进行单个率的meta分析

  产生了森林图是这样的:

stata:如何进行单个率的meta分析

  数据表上增加了新的参数:

stata:如何进行单个率的meta分析

  这些是中间变量,不用管它。

  03绘制漏斗图

  执行代码:metafunnel r ser,结果如下图:

stata:如何进行单个率的meta分析

  还有亚组分析,我们加一个亚组:

stata:如何进行单个率的meta分析

  执行代码:metan r ser, random label(namevar=study) by(subtype),结果如下图:

stata:如何进行单个率的meta分析

stata:如何进行单个率的meta分析