Repository logo
 

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

Accepted version
Peer-reviewed

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)