EigenPrism: Inference for High Dimensional Signal-to-Noise Ratios
Autor: | Lucas Janson, Emmanuel J. Candès, Rina Foygel Barber |
---|---|
Rok vydání: | 2016 |
Předmět: |
0301 basic medicine
Statistics and Probability Inference Multivariate normal distribution 01 natural sciences Article Confidence interval 010104 statistics & probability 03 medical and health sciences Bayes' theorem 030104 developmental biology Signal-to-noise ratio Linear regression Statistics Credible interval Statistical inference 0101 mathematics Statistics Probability and Uncertainty Mathematics |
Zdroj: | Journal of the Royal Statistical Society Series B: Statistical Methodology. 79:1037-1065 |
ISSN: | 1467-9868 1369-7412 |
DOI: | 10.1111/rssb.12203 |
Popis: | Summary Consider the following three important problems in statistical inference: constructing confidence intervals for the error of a high dimensional (p > n) regression estimator, the linear regression noise level and the genetic signal-to-noise ratio of a continuous-valued trait (related to the heritability). All three problems turn out to be closely related to the little-studied problem of performing inference on the l2-norm of the signal in high dimensional linear regression. We derive a novel procedure for this, which is asymptotically correct when the covariates are multivariate Gaussian and produces valid confidence intervals in finite samples as well. The procedure, called EigenPrism, is computationally fast and makes no assumptions on coefficient sparsity or knowledge of the noise level. We investigate the width of the EigenPrism confidence intervals, including a comparison with a Bayesian setting in which our interval is just 5% wider than the Bayes credible interval. We are then able to unify the three aforementioned problems by showing that EigenPrism with only minor modifications can make important contributions to all three. We also investigate the robustness of coverage and find that the method applies in practice and in finite samples much more widely than just the case of multivariate Gaussian covariates. Finally, we apply EigenPrism to a genetic data set to estimate the genetic signal-to-noise ratio for a number of continuous phenotypes. |
Databáze: | OpenAIRE |
Externí odkaz: |