2016-04-25 11 views
2

異なる研究のp値をメタ分析しようとしています。私はp-value3p-values1を分析したメタの4番目の列を取得しようとしていたデータフレームに異なる観測のp値をメタ分析する方法

DF1 
p-value1 p-value2 pvalue3 m 
0.1   0.2  0.3  a 
0.2   0.3  0.4  b 
0.3   0.4  0.5  c 
0.4   0.4  0.5  a 
0.6   0.7  0.9  b 
0.6   0.7  0.3  c 

を持っています。 は私が

p<–rbind(DF1$p-value1,DF1$p-value2,DF1$p-value3) 
pv–split (p,p$m) 
library(metap) 
for (i in 1:length(pv)) 
{pvalue <- sumlog(pv[[i]]$pvalue)} 

のMetAPパッケージを使用しようとしました。しかし、それは1つのp値になります。何か助けてくれてありがとう。 `作業:

+2

を試すことができますが、'(、1、sumlog DF1 [3,1])を適用するのでしょうか? – Sotos

+0

ああ、ありがとう! – HoHoHo

答えて

2

あなたは

apply(DF1[,1:3], 1, sumlog) 
関連する問題