Meta分析合并P值
2015-08-09 20:22阅读：
Meta分析合并P值
很多关于Meta分析的国内著作或者相关文献，都是从P值合并入手，这个操作比较简单，本文将运用R和Stata软件介绍具体的操作步骤。
1.R软件中的操作
在R软件中，需要用到metap这个数据包，首先需要安装这个数据包，install.packages('metap')，然后加载这个数据包，library(metap)，我们采用R软中自带的数据，data(beckerp)[1]，R软件提供了多种合并的方式，主要包括使用logit方法(logitp:
Combine p values using logit
method)，求均数的方法(meanp
Combine p values by the mean p
method)等方法。我们用logit方法举例子，用到的命令是logitp。
具体过程：
d
ata(beckerp)
logitp(beckerp)
我们能得出结果，t=1.616539,
p=0.0584029。
2.Stata软件中的操作
在Stata软件中输完数据之后（具体的格式如图1），用命令metapp就可以完成操作。
得到的结果是chi2=18.53301，p=0
.04661109，结果与R软件中输出的结果有出入，两者用到的方法不一致，Stata软件是Fisher。
另外，我们认为上述的P值都是单侧的，但是大多数文章中报道的P值为双侧检验，我们假设单侧研究是P单，双侧研究是P双，可以通过公示换算：P单=P双/2(当结果为阳性时)，P单=1-P双/2(当结果为阴性时)。
参考文献
[1]Becker, B J. Combiningsignificance levels. In
Cooper, H and Hedges, L V, editors A handbook of researchsynthesis,
chapter 15, pages 215–230. Russell Sage, New York,
1994.
来源：http://mp.weixin.qq.com/s?__biz=MzA5MzEwNjM5MQ==&mid=205441439&idx=1&sn=3e32727dbb506c1a58cca4336c8ea074#rd
欢迎关注公众微信号：
实用Meta分析（meta-show）
QQ讨论群：
422021477
