sample_gamma_model function

Sample gamma model(s) with parameters according to priors