sample_glm function

Sample the conditional distribution created by a generalized linear model