Unscaled Bayes factors for multiple hypothesis testing in microarray experiments
Unscaled Bayes factors for multiple hypothesis testing in microarray experiments
复制标题
DOI：
10.1177/0962280212437827
复制
发表时间：
2015-12-01
影响因子：
2.3
通讯作者：
Racugno, Walter
中科院分区：
文献类型：
--
作者：
Bertolino, Francesco;Cabras, Stefano;Racugno, Walter
多重假设检验收集了一系列通常基于 p 值的技术，作为来自许多统计检验的可用证据的摘要。在假设检验中，从贝叶斯的角度来看，针对替代方案的特定假设的证据（以数据为条件）由贝叶斯因子给出。在本研究中，我们基于贝叶斯因子和 p 值进行多重假设检验，将多重假设检验视为多重模型选择问题。为了获得贝叶斯因子，我们假设通常不正确的默认先验。在这种情况下，由于先验伪常数的比率，贝叶斯因子通常是不确定的。我们证明，忽略先验伪常数会导致未缩放的贝叶斯因子，这不会使多重假设检验中的推理过程无效，因为它们是在比较方案中使用的。事实上，使用 p 值中的部分信息，我们能够近似未缩放贝叶斯因子的采样零分布，并在 Efron 的多重测试过程中使用它。模拟研究表明，在正常抽样模型下，即使样本量较小，我们的方法提供的假阳性和假阴性比例也低于其他仅基于 p 值的常见多重假设检验方法。所提出的程序在两项模拟研究中得到说明，并且在两项微阵列实验的分析中显示了其使用的优点。
Multiple hypothesis testing collects a series of techniques usually based on p-values as a summary of the available evidence from many statistical tests. In hypothesis testing, under a Bayesian perspective, the evidence for a specified hypothesis against an alternative, conditionally on data, is given by the Bayes factor. In this study, we approach multiple hypothesis testing based on both Bayes factors and p-values, regarding multiple hypothesis testing as a multiple model selection problem. To obtain the Bayes factors we assume default priors that are typically improper. In this case, the Bayes factor is usually undetermined due to the ratio of prior pseudo-constants. We show that ignoring prior pseudo-constants leads to unscaled Bayes factor which do not invalidate the inferential procedure in multiple hypothesis testing, because they are used within a comparative scheme. In fact, using partial information from the p-values, we are able to approximate the sampling null distribution of the unscaled Bayes factor and use it within Efron's multiple testing procedure. The simulation study suggests that under normal sampling model and even with small sample sizes, our approach provides false positive and false negative proportions that are less than other common multiple hypothesis testing approaches based only on p-values. The proposed procedure is illustrated in two simulation studies, and the advantages of its use are showed in the analysis of two microarray experiments.