This figure is the directed acyclic graph used to identify variables to adjust for in the multivariable analysis
Hymel et al. (Fri,) studied this question.