StarTrek: Combinatorial Variable Selection with False Discovery Rate Control

TL;DR

This paper introduces the StarTrek filter, a novel method for selecting hub nodes in high-dimensional networks while controlling the false discovery rate, using Gaussian multiplier bootstrap and new theoretical bounds.

Contribution

The paper develops a new inferential method for identifying network hubs with FDR control, addressing combinatorial and dependence challenges with novel Gaussian comparison bounds.

Findings

StarTrek filter effectively controls FDR in simulations.

Method successfully identifies key hub nodes in real data.

Theoretical bounds ensure accurate FDR control under dependence.

Abstract

Variable selection on the large-scale networks has been extensively studied in the literature. While most of the existing methods are limited to the local functionals especially the graph edges, this paper focuses on selecting the discrete hub structures of the networks. Specifically, we propose an inferential method, called StarTrek filter, to select the hub nodes with degrees larger than a certain thresholding level in the high dimensional graphical models and control the false discovery rate (FDR). Discovering hub nodes in the networks is challenging: there is no straightforward statistic for testing the degree of a node due to the combinatorial structures; complicated dependence in the multiple testing problem is hard to characterize and control. In methodology, the StarTrek filter overcomes this by constructing p-values based on the maximum test statistics via the Gaussian multiplier bootstrap. In theory, we show that the StarTrek filter can control the FDR by providing accurate bounds on the approximation errors of the quantile estimation and addressing the dependence structures among the maximal statistics. To this end, we establish novel Cram\'er-type comparison bounds for the high dimensional Gaussian random vectors. Comparing to the Gaussian comparison bound via the Kolmogorov distance established by \citet{chernozhukov2014anti}, our Cram\'er-type comparison bounds establish the relative difference between the distribution functions of two high dimensional Gaussian random vectors. We illustrate the validity of the StarTrek filter in a series of numerical experiments and apply it to the genotype-tissue expression dataset to discover central regulator genes.