restrict_coef function

Restrict coefficients in the model