GLAM_2d_no_covariates function

Fit the 2d GLAM without covariates