Goodness-of-fit testing in high dimensional generalized linear models


Type
Article
Change log
Authors
Janková, J 
Shah, RD 
Bühlmann, P 
Samworth, RJ 
Abstract

We propose a family of tests to assess the goodness-of-fit of a high-dimensional generalized linear model. Our framework is flexible and may be used to construct an omnibus test or directed against testing specific non-linearities and interaction effects, or for testing the significance of groups of variables. The methodology is based on extracting left-over signal in the residuals from an initial fit of a generalized linear model. This can be achieved by predicting this signal from the residuals using modern flexible regression or machine learning methods such as random forests or boosted trees. Under the null hypothesis that the generalized linear model is correct, no signal is left in the residuals and our test statistic has a Gaussian limiting distribution, translating to asymptotic control of type I error. Under a local alternative, we establish a guarantee on the power of the test. We illustrate the effectiveness of the methodology on simulated and real data examples by testing goodness-of-fit in logistic regression models. Software implementing the methodology is available in the R package `GRPtests'.

Description
Keywords
Debiasing, Generalized linear models, Goodness-of-fit testing, Group testing, High dimensional data, Residual prediction
Journal Title
Journal of the Royal Statistical Society. Series B: Statistical Methodology
Conference Name
Journal ISSN
1369-7412
1467-9868
Volume Title
82
Publisher
Oxford University Press (OUP)
Rights
All rights reserved
Sponsorship
Engineering and Physical Sciences Research Council (EP/N031938/1)
Engineering and Physical Sciences Research Council (EP/P031447/1)
Engineering and Physical Sciences Research Council (EP/R014604/1)
Engineering and Physical Sciences Research Council (EP/R013381/1)