add_covariate function

Add a covariate effect to a parameter model