Relationship between collider bias and interactions on the log-additive scale.

Abstract

Collider bias occurs when conditioning on a common effect (collider) of two variables $X,Y$. In this article, we quantify the collider bias in the estimated association between exposure $X$ and outcome $Y$ induced by selecting on one value of a binary collider $S$ of the exposure and the outcome. In the case of logistic regression, it is known that the magnitude of the collider bias in the exposure-outcome regression coefficient is proportional to the strength of interaction ${\delta }_3$ between $X$ and $Y$ in a log-additive model for the collider: $P(S=1|X,Y)=\exp \{{\delta }_0+{\delta }_{1}X+{\delta }_{2}Y+{\delta }_{3}XY\}$. We show that this result also holds under a linear or Poisson regression model for the exposure-outcome association. We then illustrate numerically that even if a log-additive model with interactions is not the true model for the collider, the interaction term in such a model is still informative about the magnitude of collider bias. Finally, we discuss the implications of these findings for methods that attempt to adjust for collider bias, such as inverse probability weighting which is often implemented without including interactions between variables in the weighting model.