GLAM_1d_covariates function

Fit the 1d GLAM with covariates