Real-world data often exhibit nonlinear, interactive associations between predictors and an outcome. For example, in public health studies, environmental exposures and demographic variables can have nonlinear associations with health outcomes and modify the associations of other predictors. Parametric regression models facilitate interpretable inference and hypothesis testing with Type I error rate control, but may be misspecified when associations are nonlinear. In contrast, flexible nonparametric regression models such as Gaussian processes (GPs) can capture complex nonlinearity, but are often more difficult to interpret and do not readily provide hypothesis tests with Type I error rate control. To bridge this gap, we propose a hypothesis testing framework based on additive GPs restricted to pairwise interactions, allowing tests of variable associations and interaction effects with Type I error rate control. We further establish contraction rates for additive GPs, providing theoretical support for their reduced dependence on predictor dimension when higher-order interactions are absent. Through simulations, we show that restricting the model to pairwise interactions can substantially improve testing power relative to methods based on standard GPs. Applications to environmental health studies identify nonlinear associations between metal exposures and health outcomes as well as effect modification by demographic factors.

翻译：暂无翻译