私は、メタ分析のためにStouffer-Liptak法を使用してP値を組み合わせる必要がある立場にあります。私のデータは、行が遺伝子で、列にP値が含まれているデータフレームで構成されています。以下のシミュレートされた図。
Gene Study 1 Study 2 Study 3
a 0.03 0.001 0.3
b 0.2 0.1 0.02
c 0.04 0.05 0.03
合計で25000の遺伝子があり、Stouffer-Liptakメソッドを使用して、25,000の遺伝子すべての結合されたP値を1行に1つずつ生成するコードのスニペットが必要です。RとBioconductorを使用してこの方法を適用する方法についてのポインタをいただければ幸いです。
乾杯、アンクル。