Deeply Digging the Interaction Effect in Multiple Linear Regressions Using a Fractional-Power Interaction Term

Xinhai Li,Baidu Li,Guiming Wang,Xiangjiang Zhan,Marcel Holyoak

DOI: https://doi.org/10.1016/j.mex.2020.101067

2020-01-01

MethodsX

Abstract:In multiple regression Y similar to beta(0) + beta X-1(1) + beta X-2(2) + beta X-3(1) X-2 + epsilon., the interaction term is quantified as the product of X-1 and X-2. We developed fractional-power interaction regression (FPIR), using beta X-1(M) X-2(N) as the interaction term. The rationale of FPIR is that the slopes of Y-X-1 regression along the X-2 gradient are modeled using the nonlinear function (Slope = beta(1) + beta 3MX1M-1 X-2(N)), instead of the linear function (Slope = beta(1) + beta X-3(2)) that regular regressions normally implement. The ranges of M and N are from -56 to 56 with 550 candidate values, respectively. We applied FPIR using a well-studied dataset, nest sites of the crested ibis (Nipponia nippon).We further tested FPIR by other 4692 regression models. FPIRs have lower AIC values (-302 +/- 5003.5) than regular regressions (-168.4 +/- 4561.6), and the effect size of AIC values between FPIR and regular regression is 0.07 (95% CI: 0.04-0.10). We also compared FPIR with complex models such as polynomial regression, generalized additive model, and random forest. FPIR is flexible and interpretable, using a minimum number of degrees of freedom to maximize variance explained. We have provided a new R package, interactionFPIR, to estimate the values of M and N, and suggest using FPIR whenever the interaction term is likely to be significant. Introduced fractional-power interaction regression (FPIR) as Y similar to beta(0) + beta X-1(1) + beta X-2(2) + beta X-3(1)M X-2(N) + epsilon to replace the current regression model Y similar to beta(0) + beta X-1(1) + beta X-2(2) + beta X-3(1) X-2 + epsilon; Clarified the rationale of FPIR, and compared it with regular regression model, polynomial regression, generalized additive model, and random forest using regression models for 4692 species; Provided an R package, interactionFPIR, to calculate the values of M and N, and other model parameters. (C) 2020 The Authors. Published by Elsevier B.V.