Confounding was assessed using a hybrid method that combines previous knowledge using directed
acyclic graphs (45) and a statistical method on change in point estimates in which the potential covariate was not retained in final models, if it resulted in a change in the b-coefficient of