我需要使用 Stouffer-Liptak 方法组合 P.values 进行荟萃分析。我的数据由一个数据框组成,其中行是基因,列包含 P.values。下面的模拟图。
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.values,每行一个。任何有关如何使用 R 和 Bioconductor 应用此方法的指针将不胜感激。
干杯,安库尔。