A simple new approach to variable selection in regression, with application to genetic fine mapping

• 成果类型: Article
• 署名作者: Wang, Gao; Sarkar, Abhishek; Carbonetto, Peter; Stephens, Matthew
• 署名单位: University of Chicago
• 刊物名称: JOURNAL OF THE ROYAL STATISTICAL SOCIETY SERIES B-STATISTICAL METHODOLOGY
• ISSN/ISSBN: 1369-7412
• DOI: 10.1111/rssb.12388
• 发表日期: 2020
• 页码: 1273-1300
• 关键词: genome-wide association; False Discovery Rate; variational inference; causal variants; joint analysis; r-package; loci; optimization; expression; prediction

摘要

We introduce a simple new approach to variable selection in linear regression, with a particular focus onquantifying uncertainty in which variables should be selected. The approach is based on a new model-the 'sum of single effects' model, called 'SuSiE'-which comes from writing the sparse vector of regression coefficients as a sum of 'single-effect' vectors, each with one non-zero element. We also introduce a corresponding new fitting procedure-iterative Bayesian stepwise selection (IBSS)-which is a Bayesian analogue of stepwise selection methods. IBSS shares the computational simplicity and speed of traditional stepwise methods but, instead of selecting a single variable at each step, IBSS computes adistributionon variables that captures uncertainty in which variable to select. We provide a formal justification of this intuitive algorithm by showing that it optimizes a variational approximation to the posterior distribution under SuSiE. Further, this approximate posterior distribution naturally yields convenient novel summaries of uncertainty in variable selection, providing a credible set of variables for each selection. Our methods are particularly well suited to settings where variables are highly correlated and detectable effects are sparse, both of which are characteristics of genetic fine mapping applications. We demonstrate through numerical experiments that our methods outperform existing methods for this task, and we illustrate their application to fine mapping genetic variants influencing alternative splicing in human cell lines. We also discuss the potential and challenges for applying these methods to generic variable-selection problems.