Hierarchical correction of p-values via an ultrametric tree running Ornstein-Uhlenbeck process

Abstract

Statistical testing is classically used as an exploratory tool to search for association between a phenotype and many possible explanatory variables. This approach often leads to multiple testing under dependence. We assume a hierarchical structure between tests via an Ornstein-Uhlenbeck process on a tree. The process correlation structure is used for smoothing the p-values. We design a penalized estimation of the mean of the Ornstein-Uhlenbeck process for p-value computation. The performances of the algorithm are assessed via simulations. Its ability to discover new associations is demonstrated on a metagenomic dataset. The corresponding R package is available from https://github.com/abichat/zazou.

22 pages, 6 figures

References in corpus (6)

• In silico prediction of protein-protein interactions in human macrophages
• Tuning parameter selection in high dimensional penalized likelihood
• Detection of adaptive shifts on phylogenies using shifted stochastic processes on a tree
• Post hoc inference via joint family-wise error rate control
• Hierarchical inference for genome-wide association studies: a view on methodology with software
• Spatial clustering of array CGH features in combination with hierarchical multiple testing