add_covariates function

Add Covariates for Model Adjustment