Stata率的meta分析.doc

用stata软件做单个样本率的meta分析

本实例采用的数据是本人的另一个贴子 用的数据是这个贴子中我能用stata做患病率的meta分析了，大家交流交流啊 - 丁香园论坛 中网友在别的文章中看到并上传的一个森林图：

这个森林图中的数据其实不太好，因为异质性太大，但数据简单，也就给大家摸拟一下。重要的是知道怎样用stata软件做单样本率meta分析就可以了。
在stata中要做meta分析，最重要的就是要知道两个变量，一个是ES也就是效应量，另一个是seES,也就效应量，所有关于率的meta,做meta的关键也就是如何去寻找这两个东东了,
我这里采用的就是直接用率做为ES,而ES的标准误其实也不难求出,大家可以看看孙振球教授《医学统计学》中的这个例子或许会有所启发:如下图

这个是基于正态近似法的公式,在样本量较,数据正态时使用,这算ES的标准误也是如用的这个公式
下面开始具体操作:
1,输入数据,数据的格式是study,率，以及样本量

第二步：generate ser=sqrt(r*(1-r)/n)

3，用随机效应模型进行分析命令如下：metan r ser, random label(namevar=study)

继续

4，输入如下命令得到漏斗图:metafunnel r ser