我有一个数据:
controls cases
0.7433018893 0.5258425596
0.6592615656 0.4859544482
0.3695427609 0.3641427198
0.0508736108 0.171792228
0.6630316215 0.5950126248
0.4571880027 0.8807510188
0.0469831145 0.7860116142
我想将控件的p值(例如,0.7433018893)与其对应的/ pair / p值组合(例如,0.5258425596)。
必需的输出:
controls cases Combine_p_value
0.7433018893 0.5258425596 xxxxxxxx
0.6592615656 0.4859544482 xxxxxxxx
我正在使用R package metamp
,其代码如下:
x <- c(0.7433018893,0.5258425596)
sumlog(x)
chisq = 1.878813 with df = 4 p = 0.7580359
如何在R中自动执行此过程?如果你有一个列表pf p值,有一些方法可用,但在这里我们有配对的p值列表?任何帮助赞赏.... ,,, !!!!
Waqas。