Variable Selection for BART: An Application to Gene Regulation

Thumbnail Image
Penn collection
Statistics Papers
Degree type
variable selection
nonparametric regression
Bayesian learning
machine learning
permutation testing
decision trees
gene regulatory network
Life Sciences
Medicine and Health Sciences
Physical Sciences and Mathematics
Grant number
Copyright date
Related resources
Bleich, Justin
Kapelner, Adam
George, Edward I
Jensen, Shane T

We consider the task of discovering gene regulatory networks, which are defined as sets of genes and the corresponding transcription factors which regulate their expression levels. This can be viewed as a variable selection problem, potentially with high dimensionality. Variable selection is especially challenging in high-dimensional settings, where it is difficult to detect subtle individual effects and interactions between predictors. Bayesian Additive Regression Trees [BART, Ann. Appl. Stat. 4 (2010) 266–298] provides a novel nonparametric alternative to parametric regression approaches, such as the lasso or stepwise regression, especially when the number of relevant predictors is sparse relative to the total number of available predictors and the fundamental relationships are nonlinear. We develop a principled permutation-based inferential approach for determining when the effect of a selected predictor is likely to be real. Going further, we adapt the BART procedure to incorporate informed prior information about variable importance. We present simulations demonstrating that our method compares favorably to existing parametric and nonparametric procedures in a variety of data settings. To demonstrate the potential of our approach in a biological context, we apply it to the task of inferring the gene regulatory network in yeast (Saccharomyces cerevisiae). We find that our BART-based procedure is best able to recover the subset of covariates with the largest signal compared to other variable selection methods. The methods developed in this work are readily available in the R package bart Machine.

Date Range for Data Collection (Start Date)
Date Range for Data Collection (End Date)
Digital Object Identifier
Series name and number
Publication date
Journal title
The Annals of Applied Statistics
Volume number
Issue number
Publisher DOI
Journal Issue
Recommended citation