Adding covariates to a model