converted downstream to using covariate list instead of model matrix