sample_lm_gprior function

Sample the linear regression parameters assuming a g-prior