单臂,单个率的meta分析我们经常使用到,因为缺乏对照组,我们不能使用revman进行处理,下面小编介绍,如何使用stata进行单个率的meta分析。
01输入数据
我们使用下面的数据进行计算,其中a是时间数、n是总数、r是比例。
首先在stata中输入我们的数据
Study的格式是这样的
其他的格式是这样的
下面执行代码:. generate ser=sqrt(r*(1-r)/n),这一步是进行数据转换的。
02制作森林图
执行代码metan r ser, random label(namevar=study)。
固定效应模型是这个:metan r ser, fixed label(namevar=study)。
我们执行随机效应模型,结果如下图:
产生了森林图是这样的:
数据表上增加了新的参数:
这些是中间变量,不用管它。
03绘制漏斗图
执行代码:metafunnel r ser,结果如下图:
还有亚组分析,我们加一个亚组:
执行代码:metan r ser, random label(namevar=study) by(subtype),结果如下图: