add_confregions function

Compute the confidence regions of covariate effects